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(Received June 1 
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Hydrodynamic expansion of the solar corona is the basis for ‘‘solar corpuscular radia- 


tion.” 


The quiet day coronal temperatures of 1 to 2> 
hundred km/sec and a density of 5 to 50 particles ems at the orbit of Earth. 


< 108 °K yield a solar wind of several 
The solar 


wind draws out into space the lines of force of the general one-gauss solar field, to give an 
interplanetary field that is basically spiral in character with a density of the order of a few 


times 10~ gauss at the orbit of Earth. 


atures of 4X 10° °K or more immediately following a large flare. 
leads to a 1 to 2X 108 km/sec blast wave into planetary space. 


The enhanced corona may sometimes have temper- 


Such explosive heating 
The blast wave sweeps 


up the quiet-day solar wind ahead and kinks the quiet-day field, raising the density of the 


field to several times 10-4 gauss on some occasions. 
kink, producing a Forbush-type intensity decrease. 


swept back by the 


The galactic cosmic ray particles are 
Altogether, the 


blast wave and its magnetic field constitute the “enhanced solar corpuscular radiation.” 
It is shown that the alternative magnetic tongue model of Gold and others [1960] is 


untenable because the ‘‘tongue”’ 


of the blast wave. 


1. Introduction 


It was pointed out several years ago [Parker, 
1958a] that the solar corona is too hot to ‘be confined 
by the solar gravitational field. In place of the 
older view of the corona as an extended, but static, 
atmosphere it was shown that the corona must 
expand continually into space. A study of the 
hydrodynamic equations showed that the solar 
corona is steadly expanding with a velocity that 





is only a few kilometers per second in the lower | 


becomes supersonic at large 
distances from the sun. The result is a continuing 
outward flow of coronal gas through interplanetary 
space, which we have called the solar wind to empha- 
size its hydrodynamic nature. The observed quiet- 
day coronal temperatures of 210° °K indicate 
that the quiet-day solar wind velocity is of the order 
of 150 to 500 km/see. The density of the solar wind 
at the orbit of Earth is determined by the existing 
temperature observations to lie in the range 5 to 200 
atoms/em*®. Thus the quiet-day solar wind re- 
presents a proton flux in the range 10® to 10'°/em? 
sec at the orbit of Earth. 

Some time before these calculations were made 
it had been pointed out by Biermann [1951, 1952, 
1957] that the observed outward acceleration and 
the observed ionization and excitation of type I 
comet tails can be explained only as the result of 
continuing quiet-day solar corpuscular radiation. 
He originally estimated that the corpuscular flux 
was of the order of 10'° protons/cm? sec at the orbit 
of Earth, recently [Biermann, 1960] revising the 


1 Paper presented at 1961 Spring meeting of the URSI, Washington, D.C. 


corona but which 





would not reach the earth until many hours after the arrival 


estimate downward to 10° protons/em? sec on the 
basis of laboratory measurements of the appro- 
priate charge exchange cross sections. 

Upon completion of the solar wind calculations it 
was at once obvious that Biermann’s solar cor- 
puscular radiation was in fact the hydrodynamic 
solar wind from the steadily expanding corona. 

There are many other phenomena, such as the 
aurora, geomagnetic activity, etc., which are con- 
ventionally ascribed to solar corpuscular radiation. 
They too are the consequence of the solar wind. On 
an interplanetary scale the solar wind is hydro- 
dynamic, but on a planetary scale its low density 
renders it corpuscular in many respects. 

Since this first identification of solar corpuscular 
radiation with the expanding corona, we have ex- 
plored a number of consequences of the solar wind 
[Parker, 1958 b, ec; 1960, 1961 a, b]. It now appears 
that the solar wind is the dominant interplanetary 
dynamical force and is responsible for the inter- 
planetary magnetic field configuration, the observed 
modulation of the galactic and solar cosmic ray 
intensity, the quiet day and the enhanced geomag- 
netic activity, ete. 

Recent direct observation of the solar wind near 
Earth [Bridge, 1961] has determined its quiet-day 
velocity to lie in the range 250 to 400 km/sec, with 
a density of 10 to 20 protons/em*. The total flux is 
of the order of 0.5 10° protons/em? sec. An observ- 
ation of the quiet-day solar wind in January and 
September of 1959 had been reported earlier by 
Shklovskii et al. [1960] to suggest a flux of 0.2 10° 
protons/em? sec, but the rather large corrections 
that had to be applied to the data make it uncertain 
as to how meaningful the flux estimate really was. 
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2. The Expanding Corona 
2.1. Quiet Sun 


The simplest model of the quiet corona is one 
which assumes spherical symmetry about the center 
of the sun. Then ifr is the radial expansion velocity 
as a function of distance from the center of the sun, 
the hydrodynamic equation for stationary expan- 
sion is 

my 1 dp,GMo 


"iy NM dr. 7 r —=), (1) 


where N is the number of atoms per unit volume, 
G is the gravitational constant, and p is the hy- 
drostatic pressure. Since the coronal gases are fully 
ionized and largely hydrogen, M is the mass of the 
hydrogen atom and the hydrostatic pressure is 
p=2NkT. Conservation of mass requires that 


Nor?= Nwvoa 2 (2 ) 


where the subscript zero denotes the value at the 
reference level r=a. Choosing a=10° km, N, is 
of the order of 10’/em*. It will be sufficient for our 
purposes here to discuss the simple case of an 
isothermal corona. The more general case has been 
considered elsewhere [Parker, 1960]. 

For an isothermal corona, T= 7)~%210° °K, and 
(1) can be integrated to give 


sp 2AT ty © AE yy 
5 (0 rs. hh 





r ree o (1 (1 —*)=0 


(3) 


upon using (2) to eliminate N. The result is a one 
parameter family of curves v (r) for any given value 
of 7), with v as the parameter. The general form 
of the family is sketched in figure 1. The solution 
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The family of solutions of equation (3) for a given 
coronal temperature T.. 
The solution of physical interest is the one passing through the critical point, 
given by eq (4). The critical point is analogous to the point of sonic transition in 
a Laval nozzle. 


FiGcure 1. 
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of physical interest is the solution starting from the 
origin and passing up through the critical point to 


supersonic velocity at infinity. The coordinates of 
the critical point are (v,, r,) where 


1 
5 


~ 





Mv?=kT, GM © M/r.=4kT. (4) 
The value of v for this solution is the lower positive 
root of 


Mes_), Meé__,, GMoOM 


, 1, GMOM 
OkT "2k “ee 


akT ~ 4akT 











(5) 





The coronal expansion will automatically drift into 
this solution whatever its initial expansion rate: 
The solutions with smaller v go to r= on the 
lower branch of v (7), which can be maintained only 
by a large inward pressure from v=o. In the 
absence of such an inward pressure the expansion 
will accelerate to the solut on through the critical 
point. There are no solutions above the solution 
through the critical point which start at the origin. 
Thus, the only available solution for expansion in 
the absence of an inward pressure at r= is the 
solution through the critical point, which starts 
with the velocity % given by (5). It has been 
pointed out by Clauser [1960] that this expansion 
of the corona is analogous to the expansion of a 
gas through a Laval nozzle, with the gravitational 
field playing the role of the throat. 

The expansion velocity as a function of radial 
distance is shown in figure 2 for a number of coronal 
temperatures. The Observed coronal temperatures 
are of the order of 210° °K [Billings, 1959], but 
even low temperatures like 110° °K give solar 
wind velocities of the order of several hundred 
kilometers per second. 
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Ficure 2. The steady expansion velocity as a function of 
radial distance for a model isothermal corona. 
The two lower curves, for 1 and 2108 °K, are applicable to the quiet corona. 
The upper curve, for 4X10®,°K, gives an indication of the rate of expansion of 
the enhanced corona. 
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2.2. Active Sun 


It is observed that following a large solar flare the 
coronal temperature over the general region of the 


flare may increase to 410° °K or more [Christian- | 


sen, Yabsley, and Mills, 1949; Christiansen and 
Warburton, 1953]. The result is a hydrodynamic 
explosion of the corona outward into space. The 
expansion velocity may be in excess of 1000 km/sec, 
as is evidence from figure 2. The asymptotic form 
of the resulting blast wave from the sun may be 
treated using the similarity transformations of the 
progressive wave [Courant and Friedrichs, 1949]. 
The idealization is made that the corona has spheri- 
cal symmetry; i.e., the entire corona is heated to 
410° °K, rather than the portion over the active 
region in which the flare occurred. The 300 km/sec 
velocity of the quiet day solar wind ahead of the 
blast wave is neglected, and the temperature of the 
quiet day wind is assumed to be so small that the 
10° km/see velocity of the blast wave represents a 
very high Mach number. The resulting blast wave 
profiles are shown in figure 3. The shock transition 
at the front is of the collisionless type and involves 
a density increase by a factor of four above the 
quiet-day value. The shape of the blast wave 
behind the shock transition depends upon how hard 
the enhanced corona pushes from behind. In the 
event that the corona should push so hard that the 
energy of the blast wave increases linearly with time 
after leaving the sun, then the relatively thin, high- 
density profile identified by \=1 in figure 3 is the 
result. If the corona should not push at all once 
the blast wave is on its way, then a linear decrease 
of density behind the front is the result, identified 
by the sawtooth profile for \=3/2 in figure 3. In 
either case an enhanced coronal temperature of 
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Ficure 3. Blast wave profiles from the active corona. 


The quiet day solar wind density is given by po (r), and is essentially propor- 
tional to 1/r?. The shock transition at the head of the blast wave is at r=, 
where the density jumps by a factor of four from the quiet day value. The 
parameter \ is an inverse measure of how strongly the enhanced corona is driving 
the blast wave from the rear: If the corona does not push at all on the blast wave, 
then A=3/2 and the sawtooth profile is the result; if the corona pushes so hard 
that the energy of the blast wave increases linearly with time, then A=1 and the 
blast wave lies between r/R, =0.84 and r/Ri=1.0. 





410° °K near the sun yields a blast wave with a 
velocity of 1-210° km/sec, and a density of the 
order of 10? em at the orbit of Earth [Parker, 1961, 
a and bj. Thus the blast wave represents the 
“enhanced solar corpuscular radiation” responsible 
for the geomagnetic storm, the Forbush decrease, 
ete. 

An important point to be noted is that’ the blast 
wave consists of gas which had already left the sun 
before the occurrence of the flare. It is material 
which has been swept up as the wave advances out- 
ward into space. The material in the corona at the 
time of the flare is behind the rear of the blast wave 
profile. This must be borne in mind in the next 
section when we discuss the interplanetary magnetic 
field configuration. 


3. Interplanetary Magnetic Fields 
3.1. Quiet Sun 


The magnetic fields in interplanetary space are 
composed of the magnetic lines of force of the 
general solar field [Babcock and Babcock, 1955; 
Babcock, 1959] which are pulled outward from the 
sun by the expanding corona. The magnetic fields 
associated with active regions on the sun are usually 
considerably in excess of 1 gauss, so that they can 
to a large degree withstand the expansion of the 
local corona, remaining fixed over the active regions. 
The hydrostatic pressure in the solar corona is of 
such an order, 10~-? dynes/em?, as to suggest that 
only fields of the order of one gauss can be extended 
by the solar wind. The simplest model, then, is 
one in which a corona with spherical symmetry 
expands outward, carrying with it the dipole-like 
fields of the polar regions of the sun. Such an 
idealized picture ignores the many complications of 
the disordered one gauss fields which are undoubtedly 
carried outward from the equatorial regions. The 
idealization also ignores the many interesting com- 
plications that result from the fact that the actual 
corona expands more rapidly and densely in some 
directions than in others. However, the idealization 
of a radial solar wind with spherical symmetry in 
a dipole field serves to illustrate the general principles 
and to give a rough idea of the overall interplanetary 
magnetic configuration to be expected. When more 
comprehensive observations allow greater detail to 
be filled into the present simple theoretical models, 
the many irregularities and swirls which undoubtedly 
exist in the interplanetary field may be taken into 
account. Ignoring the irregularities now will give 
a conservative estimate of the cosmic ray effects to 
be discussed in the next chapter. 

The lines of force of a general solar dipole field 
represent the two parameter family of curves 


r—az - (¢—¢0), 6=4 


upon extension into space by the quiet day solar 
wind of velocity y. The solar angular velocity Q is 
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ORBIT OF EARTH 






300 km/s 


Figure 4. Projection onto the equatorial plane of the lines 
of force of the solar fields extended by a quiet-day radial solar 
wind of 300 km/sec. 


(00 ’ do) 


approximately 2.610~* radians/sec, and 


is the position, in spherical polar coordinates, of | 


the intersection of the line of force with the reference 
level r—a in the lower corona. The radial compon- 
ent of the field is 


B,(r, =)=B,(a, 6) (*) 


2 
’ 


the meridional component is zero, and the azimuthal 
component 1s 
Bz(r, 0=B,(a, 0)rQ/v. 


Projection of each line of force onto the solar equa- 
torial plane is an Archimedes spiral, shown in figure 
4 for y=300 km/sec. It is readily seen that the 
field is principally radial inside the orbit of Earth 
and principally azimuthal beyond. One gauss at 
the solar photosphere yields about 210~* gauss in 
the radial direction at a distance of one astronomical 
unit. 


3.2. Active Sun 


The quiet-day interplanetary field shown in 
figure 4 is deformed by the blast wave from the 
corona over a flare as shown in figure 5. Again the 
case \=1 represents the situation when the enhanced 
corona is pushing hard on the rear of the blast wave, 
and \=3/2 when the blast wave is coasting. In 
either case the field may reach 10 107° gauss at the 
front of the wave. When the corona is pushing hard 
on the rear of the blast wave the field toward the 
rear may easily be 40 10~° gauss or more. 
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It has been suggested that the “enhanced solar 
corpuscular radiation,” which we suggest is a hydro- 
dynamic blast wave from the corona, sometimes 
draws out a loop of field from the active region, 
giving a magnetic tongue which is responsible for 
many of the effects associated with the enhanced 
| corpuscular radiation (see discussion and references 


in Gold, 1960). Presumably the expanding interior 
of the magnetic tongue is partially shielded by the 
fields of the tongue from the galactic cosmic ray 
intensity, giving the Forbush-type cosmic ray 
decrease in coincidence with the geomagnetic storm 
resulting from the impact of the gas in the tongue 
against the geomagnetic field; once the tongue has 


| engulfed Earth, it affords free access of solar protons 


| ary field on a number of points. 


| 
| 


from the sun to Earth; the tongue stores the solar 
protons following the emission from the sun, ete. 
In this way it has been argued that the tongue 
explains most of the observed cosmic ray phenomena 
associated with ‘‘enhanced solar corpuscular radi- 
ation.” 

We object to the tongue model of the interplanet- 
First of all, the 
magnetic configuration of the blast wave fields shown 
in figure 5 possess the same properties as the tongue 
so far as modulating and channeling galactic cosmic 
rays and solar protons. Thus the tongue is not 
unique and seems an unnecessary complication. 
Second, and more serious, is the objection that the 
tongue cannot possibly arrive at Earth with the 
“enhanced solar corpuscular radiation.”” Whatever 
the nature of the outburst on the sun which leads to 
the enhanced radiation, the result is a blast wave 
which scoops up the interplanetary gas ahead to 
form a wave as shown in figure 3. The active 


_ coronal material is behind the blast wave, and hence 





so is the magnetic tongue which the coronal material 
is supposed to carry. Even in the extreme case 
that the enhanced corona is pushing so hard on the 
rear of the blast wave that X=1 (which requires that 
the 10° °K coronal temperatures extend all the way 
to the rear of the wave) the tongue is 0.16 a.u. 
behind the front of the blast wave, which means a 
delay of the order of 6 hr in the arrival of the tongue. 
Yet it is clear in many cases that the cosmic ray 
modulation ete., begins with the first onset of geo- 
magnetic activity. The tongue configuration is 
shown in figure 6 for the special case that \=3/2. 
The tongue extends barely beyond the orbit of 
Mercury when the blast wave reaches the orbit of 
Earth. 


4. Cosmic Ray Modulation 
4.1. Quiet Sun 


The general outward flow of nonradial magnetic 
fields from the sun tends to push back the galactic 
cosmic ray particles and thereby decrease the cosmic 
ray intensity in interplanetary space below the in- 
terstellar level. This reduction follows the general 
solar cycle and is usually referred to as the 11-year 
variation. There are many individual effects which 
probably contribute to the total reduction: the ideal- 











Ficure 5. 
4 by the freely coasting blast wave, \=3/2, and by the driven 
blast wave \=1. 


Deformation of the quiet-day field shown in figure 


The shock transition is at r=R,. The rear of the blast wave is at »= Ro. 


ized quiet-day interplanetary field shown in figure 4 
is, in effect, a giant corkscrew conveyor which 
revolves once every 27 days with the sun and tends 
to push outward the galactic particles slidirg in- 
ward along the lines of force. There is reason to 
expect disordering in the interplanetary field, some 
examples of which have been cited elsewhere [Mever, 
Parker, and Simpson, 1956; Parker, 1958c]; the 


disordering contributes to the trapping of solar 


particles and to the impediment of galactic cosmic 
ray particles diffusing inward from _ interstellar 
space. 

These effects are large and contribute signifi- 
cantly to depression of the cosmic ray intensity. 
The depression is of the order of a factor of three 
at 1 Bev (see for instance Simpson, 1960) and prob- 
ably more at lower energies. We may hope that 
one ‘day interplanetary observation may decide 
what effects predominate. 








SUN 





Figure 6. A re-entrant loop of magnetic field in the solar 
corona before an outburst is shown, followed by its extension 
into space by the blast wave \=3/2. 

The loop remains in association with the coronal gas with which it moves 
following the outburst. The blast wave ahead of the loop consists of inter- 
planetary gas (the quiet-day solar wind) swept up to form the blast wave. 


4.2. Active Sun 


The magnetic configuration carried by the blast 
waves from the active corona, and shown in figure 5 
for an idealized case, have the basic characteristic 
of a strong magnetic field (probably 10-501075 
gauss) in the wave, with an essentially radial field 
connecting into the sun from behind. The strong 


fields result from the twisting and pinching to- 
gether of the quiet-day lines of force. The re- 


sulting constriction in the lines of force represents 
a serious impediment to the passage of cosmic ray 
particles. Thus the blast waves will have a ten- 
dency to store energetic protons of solar origm 
behind them, and to push back galactic cosmic ray 
particles ahead of them, as required by observation 
(see for instance, Steljes et al., 1961). The Forbush 
decrease for the idealized model shown in figure 5 
may be as large as 40 percent, and more if any 
complications are permitted in the idealized field 
configuration. The energy dependence of the cos- 
mic ray decrease can be estimated if it is remembered 
that the blast waves are in fact not spherical but 
occupy a solid angle of the order of one steradian. 
The diffusion of cosmic rays from the sides of the 
region swept out by the blast waves gives an energy 
dependence to the fractional decrease of the cosmic 
ray intensity as shown in figure 7. The parameter 
v in the figure is a measure of the angular width 
of the blast wave. ‘The computed energy depend- 
ence is found to be not unlike the observed energy 
dependence of the Forbush decrease (see for instance 
MacDonald and Webber, 1960) viz, proportional 
to reciprocal rigidity, or flatter. 
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Figure 7. The fractional depression Au/p of the cosmic ray 
intensity as a function of proton energy x (in units! of 931 
Mev) for a freely coasting shock, \=3/2, and for the hard 
driven shock, \=1. 


The parameter » is an inverse measure of the angular width of the blast wave as 
seen from the sun, The broken lines represent reciprocal rigidity, for purpose 
of comparison, 


5. Summary 


The purpose of the paper has been to describe 
in a semiquantitative way the solar wind model 
of interplanetary dynamical processes. The model 
affords a deductive approach to the problem of 
interplanetary plasmas, fields, and cosmic ray 
variations, starting with the hydrodynamic theory 
of the expanding solar corona. The observed 
coronal temperatures and their variations lead 
through such idealized models as we have discussed 
to the kind of “‘corpuscular’’ and cosmic ray effects 
that are known from observation. As direct obser- 
vation of interplanetary conditions become more 
comprehensive with the progress of space experi- 
mentation, we look forward to filling in the many 
details which are missing from the present simple 
model. 
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Attenvation coefficients for propagation under SID conditions are deduced, from records 


of atmospheric noise, for the frequenzy range of 3.5 to 50 kilocycles per second. 
values under normal daytime circumstances. 


compared with the 


These are 
It is shown that the 


advent of SID implies little change in attenuation between about 12 and 20 kilocycles 
per second; above this range there is a markedly decreased attenuation, while below 12 
kiloeyeles per second there is an even more pronounced increase in the attenuation coef- 


ficients. 
ated with SID conditions is outlined; 


The principle of an improved recorder of the changes in atmospheric noise associ- 
this recorder would discriminate between source 


effects independent of the SID and propagation influences solely attributable to the SID. 


1. Introduction 


The behavior, during propagation under normal 
daytime conditions, of radio waves of very low 
frequency (VLF), is now well established. It is 
known that the propagation is within the wave- 
guide formed by the earth and the lower edge of 
the ionospheric D layer. Attenuation coefficients 
in the VLF range (3 to 30 ke/s) have been deduced 
using the mode theory [Wait, 1957], and have also 
been derived experimentally from observations of 
atmospherics [Chapman and Macario, 1956] and of 
longwave transmitters [Watt and Plush, 1959]. 
The agreement between experimental and theoreti- 
cal results has been good, and the increasing rec- 
ognition of such refinements as the influence of 
ground conductivity, the effect of the earth’s mag- 
netic field, and so on, is likely to lead to still better 
accord. 

One influence that considerably affects daytime 
VLF propagation is the sudden ionospheric dis- 
turbance (SID). It is now generally agreed that 
a SID is caused by an increased flux of X- rays 
emitted by the sun during a solar flare [Friedman, 
1959]. These X-rays enhance the ionization in 
the D layer. One aspect of this increased ionization 
is that an additional shallow bulge of electrons 
forms upon the underside of the normal D region. 
The exact shape of the bulge depends upon the 
spectral distribution of the X-rays in the incident 
flux [Pierce and Swift, 1960]; in all cases, however, 
the depression of the D layer will be greatest at 
the subsolar point, and there will be appreciable 
effects upon the D region within a radius of at least 
6,000 km from the subsolar point. 

__ As regards propagation at VLF the influence of a 
SID may be represented as a narrowing of the wave- 
guide and, possibly, a change of the effective con- 
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ductivity at the upper boundary of the guide. Con- 
sequences upon propagation may therefore be 
anticipated and these are indeed found. The main 
effects are sudden phase anomalies (SPA) in which a 
phase-change in the skywave from a manmade 
transmitter is detected, due to the depression of the 
D reflecting layer; sudden anomalies in field strength 
(SFA) which result from the changed phase relation- 
ship between the skywave and the groundwave, 
again from an artificial transmitter; and sudden 
enhancements of atmospherics (SEA) in which an 
increase of the noise due to distant oo is 
observed (typically at a frequency of 27 ke/s). SPA 
and SFA studies indicate that for a specific oa ‘tron 
density an average SID lowers the D layer by some 
5 km [Bracewell and Straker, 1949]. 

Although, as already mentioned, VLF attenuation 
coefficients have been widely obtained for normal 
conditions, the same information does not appear to 
have been derived under SID circumstances. The 
essential raw material, namely the changes observed 
in field strengths, has however been available for 
some time. It is the purpose of this note to attempt 
to obtain some indication of the VLF attenuation 
coefficients, during a SID, from data upon atmos- 
pherics. 


2. Normal Daytime Attenuation Coefficients 


Figure 1 depicts the relationship between the 
attenuation coefficient and frequency during normal 
daytime propagation. Two curves are shown. The 
first curve was deduced by Chapman and Macario 
[1956] from observations of atmospherics. The 
second is that derived by Watt and Plush [1959], 
largely on the basis of observations of manmade 
transmitters, and is described as being an average 
for a mixed path. There are divergencies between 
the two curves but these are not entirely unexpected. 
Even if D layer ionospheric conditions were uniform 





AFTER CHAPMAN AND MACARIO (1956) 





----AFTER WATT AND PLUSH (1959) 





TENUATION COEFFICIENT ¢ 


AT 











FREQUENCY , kc/s 
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tions. 

during daylight irrespective of solar altitude, such 
influences as ground conductivity and the direction 
of propagation with respect to the earth’s magnetic 
field would still be important. These influences 
may lead to systematic deviations from ‘average’ 
behavior. For instance, it is not unusual for the 
main sources of atmospherics received at a particular 
station to lie within a restricted geographical area; 
this implies that the attenuation coefficients, ob- 
tained from such observations of atmospherics, will 
be representative of a certain path, of its orientation 
with regard to the magnetic field and of the con- 
ductivity along its ground boundary, rather than of 
more general conditions. 


3. Basic Data and Its Analysis 


Wait [1958] has given the following expression for 
the vertical electric field E (in millivolts/meter) at a 
great circle distance d from a source radiating P kw: 


, 300 Px ins, 
pnts Ph} exp (—ad), (1) 


where a is the radius of the earth, A the height of 
the ionospheric reflecting layer in kilometers, \ is 
the wavelength in kilometers, and a the attenuation 
factor. Equation (1) assumes that only a single 
mode is dominant, but this is a fair approximation 
up to 30 ke/s and for distances in the range 1,000 
to 2,000 km, and almost completely true for d greater 
than 2,000 km; thus (1) is valid in the study of 
distant atmospherics. Equation (1), however, is 
likely to be increasingly incorrect as the frequency 
rises above 30 ke/s. Under normal conditions 
(suffix NV) (1) may be written as 


= MD Pr Roa ‘ 
cert © i iota } i i a 


If atmospheric noise is being received, then the 
advent of SID conditions (suffix S) will not cause 
P or d to change and 


+ 300 Pr sm Xe 
Es {Pa} exp (—asd). (3) 


~ hs 
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Hence hyEy/hsEs= exp (as—ay)d, (4) 
it being assumed that Ey and Es are compared at 
the same frequency. Equation (4) may be con- 
verted into the form 


(as—ay)d=2 pod 10 logo (hy Ey/h,E,) ; (5) 


in which the units for as and ay are now the familiar 
db/1,000 km. Thus if records of Ey and Eg are 
available together with estimates for d, hy, and 
hs, (as—ay) can be deduced, and, knowing ay, as 
obtained. 

Four sets of data in the published literature give 
the ratio Ey/Es at different frequencies and are 
therefore particularly suitable for the application of 
eq (5). These instances are in the papers by Gardner 
[1950], by Obayashi, Fujii, and Kidokoro [1959], 
and by Obayashi [1960]; this last paper contains two 
examples. Internal evidence within Gardner’s paper 
suggests that 1,500 km is a good estimate for the 
distance of the noise sources; the reliability of this 
estimate is enhanced by the fact that Gardner 
worked in close collaboration with the Sferics 
Locating Network of the British Meteorological 
Office. It remains to determine hy/hs before (as— 
ay) can be ealeulated. hy/hs is not far removed 
from unity and the errors in its estimation are 
considerably less than those implicit in assigning 
a value to d. 70 km is a commonly accepted value 
for hy, while hs may be taken as (65+5) km [Brace- 
well and Straker, 1949] the extremes corresponding 
to the slightest and to the most intense SID. Thus 
hy/hs=1.08. Values of Ey/Es were obtained from 
Gardner’s paper and, using eq (5), (as—ay) was 
calculated over the frequency range 3.5 to 50 ke/s. 
The results are represented by curve (a) on figure 2. 

Similar information to that of Gardner is given 
in the much later paper by Obayashi and coworkers 
[1959]. Here the frequency range is restricted to 
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10 to 50 ke/s. There is no direct evidence in the 
paper as to the distances of the sources of atmos- 
pheric noise. However, the SID concerned occurred 
between 0430 and 0500 U.T. on 16 August 1958. 
Atmospheric noise records for this day indicate no 
local sources, while sferics bearing observations 
covering the ‘period of the SID show the activity 
to be to the east “< to the southwest of Japan 
(IGY Data, Japan, 1959]. It seems plausible to 
associate these cue respectively with the 
main storm centers in the Americas and in SE Asia. 
The activity in the Americas will be dying with the 
advent of local night; that in SE Asia will be rapidly 
increasing the time being adjacent to local noon. 
In atmospheric noise measurements the nearest 
storms are the dominant factors [Horner, 1960]; 
thus the SE Asia center is likely to be more important 
than that in the Americas. Furthermore, it is the 
closer parts of the SE Asia center that will be most 
vital; this effect is enhanced by the influence of time 
since maximum thundery activity occurs in the local 
afternoon, and it is generally in the nearer parts of 
the SE Asia center to Japan, that, at the time of 
the SID, the local maximum is most closely ap- 
proached. Considering all factors, it seems reason- 
able to identify the principal sources of the atmos- 
pherics as being within the nearest landmasses to 
the southeast of Japan, and to take d as being 2,000 
km. Using hy/hs=1.08, the values of (as—ay) 
were calculated as before and are represented by 
curve (b) on figure 2. 

Obayashi [1960] has given two further examples 
of enhancements of atmospheric noise which are 
suitable for analysis. These are for 29 March 1959 
and 16 June 1959. The distance of the sources 
of noise for the first instance is estimated at 2,000 
km, and Obayashi states that the storms were 
rather closer for the second enhancement; d=1,500 
km would seem to be a reasonable mean estimate in 
this case. Using hy/hs=1. 08 as before, the results 
obtained are shown in curves (c) and (d) on figure 2. 

The agreement between the four curves in figure 
2, for the range of 8 to 35 ke/s, is surprisingly good 
considering the speculative nature of the estimates 
for d in particular. Above 35 ke/s and below 8 
ke/s, the signal strengths involved are low, and 
thus the results are more prone to errors. This 
effect may account for the spread between 35 and 
50 ke/s, but seems quite inadequate to explain the 
discrepancies below 8 ke/s. Here the divergence 
between Gardner’s and Obayashi’s work is very 
considerable. The estimate for d is probably better 
for Gardner’s observations than for Obayashi’s; 
however, this factor is insufficient to enable the 
disagreements to be reconciled, and further experi- 
mental work is obviously desirable. 

An intriguing feature of three of the curves on 
figure 2 is a maximum value and then a drop in 
(as—ay) as frequency decreases below 10 ke/s. If 
this trend is continued it is probable that at some 
lower frequency, possibly in the extremely low fre- 
quency (ELF) range (<3 ke/s), (as—ay) will become 
negative as it is at frequencies above about 15 ke/s. 
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It might be predicted that a SID will produce an 
enhancement of atmospheric noise above 15 ke/s; a 
decrease at intermediate frequencies; and an en- 
hancement again in the ELF band. This kind of 
predicted behavior, which would have analogies with 
the usual variation of ay with frequency, could 
readily be verified by an extension of experimental 
observations to lower frequencies. 


4. Discussion 


Figure 3 is an attempt to summarize the preceding 
information. Curve (a) is an average representation 
for ay derived from the two graphs of figure 1. 
Curve (b) is a composite graph for as; it was ob- 
tained by considering the four curves of figure 2, 
deducing weighted mean values for (ay—as), and 
using these, in conjunction with curve (a) of figure 3 
to derive curve (b). 

It is apparent from figure 3 that the frequency 
of 27 ke/s which has been established for many 
years [Bureau, 1937] as an indicator of enhance- 
ments in atmospherics is a good choice, since the 
minimum of ag is close to 27 ke/s and the change 
between ay and ag is considerable. It is equally 
true that the differences between ay and ag are 
most marked for the lower frequencies, and the 
decrease in atmospherics over that frequency range 
is potentially a more sensitive indicator of the effects 
of a SID than is the enhancement at frequencies 
exceeding some 15 ke/s. This is particularly so 
if the normal storm centers are sufficiently close 
to the observing station for adequate daytime 
signal strength to be always present. An interest- 
ing feature of figure 3 is that the band (15 to 20 
ke/s) occupied by the most important VLF trans- 
mitters is little affected in attenuation by the onset 
of SID conditions. 

The normal atmospherics receiver on 27 ke/s 
which is usually intended as an indicator of the 
sudden enhancements associated with a SID, can 
also give increased signals produced by source 
effects and not by the propagation changes due 
to the SID. Three main causes may be assigned 
to such increases, namely, an augmentation of the 
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Ficure 3. Composite values for attenuation coefficients under 
normal daytime and under SID conditions. 
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activity in an existing storm center; the develop- 
ment of a new storm center; and the movement 
of existing activity closer to the observing station. 
Tt seems possible, upon the basis of figure 3, to pro- 
duce a reciprocal instrument which would not 
respond to the above source effects, but would 
only register a change in propagation conditions 
such as that produced by a SID. The principle 
would be to select two frequencies f; and fz, one on 
each side of the normal minimum for ay at around 
18 ke/s, and for which the values of ay are equal; 
as an example figure 3 indicates that 9 and 35 ke/s 
could be chosen. for f; and f,. Then if E(f,) and 
E(f2) are the measured fields at f,; and fs, the ratio 
E\ (fi) /E( fo ) should remain effectively constant irre- 
spective of alterations in the strengths of the sources 
of atmospheric noise, or of changes i in the distances 
of the sources, provided that these are not such as 
greatly to modify the validity of the basic one- 
mode assumption; this implies distances always 
greater than some 1,000 km. However if prop- 
agation conditions are affected as in a SID, the 
ratio E(f,)/E(f2) will alter very considerably be- 
cause of the changes in attenuation coefficient 
from ay to as. The effect will be especially marked 
being reciprocal at the two frequencies, as is shown 
in figure 3 for the selected frequencies of 9 and 35 
ke/s. Construction of a simple recorder presenting 
the ratio E( f,)/E( jz) directly offers no great difficulty. 

It is outside the scope of this paper to examine the 
interpretation of the as values of figure 3 upon the 
waveguide mode theory. However, there is no 
reason to believe that the values are incompatible 
with the theory, indeed the opposite is rather the 
vase and the results are further evidence—although 
this is hardly necessary—of the general validity of 
the mode interpretation. For instance, it can 
readily be deduced from the work of Wait [1957] 
that a reduction in fh from 70 to 60 km, without any 








considerable associated change in ionospheric con- 
ductivity, produces two graphs very similar to those 
depicted on figure 3. In each case, there is a recipro- 
cal relative displacement of the two curves, pivoting 
around a common point. 
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The problem of communication between antennas, submerged in a conducting medium 
such as sea water, is analyzed in terms of a dipole radiating in a conducting half space sepa- 
rated by a plane boundary from a dielectric half space. The theory is discussed for both 
horizontal and vertical, electric and magnetic dipoles. 

Expressions for the Hertzian potentials of the dipole in the conducting half space can 
be reduced to integrals obtained by Sommerfeld (for a dipole at the boundary) multiplied 
by an exponential depth attenuation factor. The Hertzian potentials are used to determine 
the electric and magnetic field components. 

This analysis shows that the main path of communication between submerged antennas 
is composed of three parts as follows: (a) energy flow from the transmitting dipole directly 
to the surface of the sea, (b) creation of a wave that travels along the surface refracting 
back into the sea, (ec) energy flow normal to the surface to the receiving dipole. 


1. Introduction 


Radiation from electric and magnetic dipoles in the air above a plane earth has been con- 
sidered by many writers beginning with Sommerfeld [1909]. Sommerfeld’s original paper 
treated a vertical monopole placed at the surface of a flat earth having arbitrary dielectric 
and conductive properties. The solution of the boundary-value problem was based on the 
evaluation of Fourier-Bessel integrals which were solutions to the wave equation. Weyl 
[1919] attacked the problem by resolving the dipole radiation into a spectrum of plane waves 
and evaluating the resulting integral without resorting to Sommerfeld’s cylindrical wave 
formulation. Sommerfeld [1926] expanded his original work to take into account vertical 
and horizontal, electric and magnetic dipoles above a plane earth. 

Sommerfeld’s original paper contained an error in sign which was a subject of much later 
discussion [Norton, 1935]. Although Sommerfeld’s work is followed here, the effect of the 
error in sign does not appear because of the assumption in this paper that the transmitting 
antenna is located in a highly conducting half space. 

In general, papers published prior to 1940 on the subject of radiation of dipoles were 
concerned witn dipoles in air or on the surface of the earth. Moreover, several papers have 
been written on this subject since 1940. Study of electromagnetic radiation in the sea began 
about the turn of the century and limited experimental work was done around the end of the 
first world war. Tai [1947] analyzed an electric dipole in an infinite, conducting medium and 
found extremely high dissipation in the vicinity of the dipole because of high conductive losses. 

The first extensive theoretical treatment of electromagnetic radiation by vertical and 
horizontal, electric and magnetic dipoles immersed in a conducting half space appeared in 
a thesis by Moore [1951]. This paper describes a portion of that work. Wait [1952] [1957 
considered an insulated magnetic dipole in an infinitely conducting medium, showing that the 
fields are independent of the characteristics of the insulation for an antenna diameter much 
less than the radiation wavelength in the conducting medium. Then Wait and Campbell 
[1953] and Wait [1959] analyzed a magnetic dipole of this type in a semi-infinite medium in- 
cluding special cases of frequency, antenna depth, and separation between antennas. Their 
model was a horizontal magnetic dipole (axis parallel to the surface of the sea), submerged 
in the sea (a conducting half space). 

Analysis of the electric dipole was done by Lien and Wait [1953] in which the evaluation 
of the complex integrals was reduced to forms suitable for the numerical computation. Bafos 
and Wesley [1953, 1954] carried out a mathematical analysis of the general case of a horizontal 
electric dipole in a conducting half space. The results in this paper agree with those of Bafios 
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to the first order approximation. Kraichman [1960] experimentally verified Wait’s and Bafios’ 
results in the intermediate distance region range. Also, the exponential increase of attenuation 
with depth was experimentally verified by Saran and Held [1960]. Anderson, in a thesis [1961], 
described the three-layered problem of earth-air-ionosphere for an electric dipole submerged 
in the sea. Further work by Wait [1960a, 1960b] has recently been carried out in an effort to 
unify the various approaches to the problem. 

It is concluded from our analysis that the path of electromagnetic energy, between trans- 
mitting and receiving dipoles in the sea, is the following: (a) propagation from the transmitter 
directly to the surface, (b) propagation along the surface of the sea allowing refraction of the 
energy back into the sea, (c) leakage directly downward into the sea to the receiver. 

A fundamental assumption in this analysis is that the displacement current (in the sea) 
is negligible compared to the conduction current. Since the conductivity of the sea is about 
4 mhos/m and the frequencies of practical use in undersea communication are less than 50 
ke/s, this is not at all a severe restriction for the case considered. Application to related 
problems depends upon the appropriateness of this assumption for each problem. The only 
other restriction is that the energy traveling directly through the sea between the transmitter 
and receiver is neglected. The ratio of the magnitude of the direct wave (through the sea) 
to the over-the-surface wave is of the order of n e~“-? where R, 2, 6, and n, are respectively 
antenna separation, antenna depth, skin depth, and index of refraction of the sea with respect 
to the air. Thus for R>>z, the ratio is considerably less than unity. 

Because of the “‘over-the-surface” mode of communication, the really important direction 
for transmitting dipole radiation is directly toward the surface of the sea. Since the radiation 
pattern of the vertical dipoles has a null in this direction, the vertical dipoles are not as effective 
in launching the required wave as are the horizontal dipoles. Electric and magnetic fields have 
been calculated for both the vertical and the horizontal dipoles, but the latter are of much 
greater importance in submarine communications. 


2. Waves in Conducting Media 


Discussions of the nature of electromagnetic radiation in a conducting medium and a 
conducting half space are presented by Stratton [1941] and Sunde [1949]. However, a brief 
review of the pertinent parts of the theory is included here for the purpose of defining the 
variables employed throughout this paper, as well as stating the fundamental conditions of the 
problem of dipole radiation in a conducting half space. 

The Maxwellian curl equations for a conducting medium in which displacement current 


can be neglected are 


VXE=—jouH (1) 
VXH=cE (2) 


where the time variation e%‘ has been suppressed 
E is the electric field vector (volts/meter) 
H is the magnetic field vector (ampere-turns/meter) 
w is the angular frequency (radians/second) 
uw is the magnetic permeability (henrys/meter) 
o is the conductivity (mhos/meter). 


The diffusion equation for E is obtained by combining (1) and (2): 
VE=jwucE. (3) 
The solution of (3) for a plane wave traveling in the positive 2-direction is: 


5 RXE 


E=E,e—“t28 with H=e-*'*, 
pwd 
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where 6=~+2/wuo is the skin depth (meters) 
n is the unit vector in the direction of propagation. 

It is convenient in determining the fields of a dipole to determine first either the vector 
potential A or the Hertzian potential 1. In this paper, the Hertzian potentials are used for 
both the electric and magnetic dipoles. The electric and magnetic field vectors are functions 
of the Hertzian potential as follows: ; 


Electric dipole Magnetic dipole 
H,=cVXITI, E;=—joyV X ft (4) 
E,=VV - 1,+471, H,=VV - Oy7+4i0; (5) 


where I,, M* are the Hertzian potentials for the electric and magnetic dipoles respectively, and 


ky=V—jouo (meters™'). 


The subscript 1 denotes a conducting medium. Throughout this paper the expressions re- 
lating to the electric dipole are presented in the left column while those for the magnetic dipole 
are in the right column, side by side, for convenient comparison. The Hertzian potentials 
II and II* (the Green’s function for a dipole source) in the infinite, conducting medium are 





e7imR e7imR ” 
II;=pi —R II}=pi - R (6) 
and 
/\ Pon 
_ fh Pe Is (7) 
Pine Pi= 4r , 


where F? is the distance from the dipole to the observation point 


Tr 
Ii=lim (J]), is the electric moment (ampere-meters) 
[->o 
1-0 
‘ a.” ° . e ‘ 
NIS=lim (NIJS) is the magnetic moment (ampere-meters *) 
[->o 
Ss—0 


J is the current (amperes). 

lis the electric dipole length vector (meters). 
N is the number of turns in the magnetic loop. 
S is the loop area vector (meters *). 


The corresponding statements for the electric and magnetic fields in an infinite, noncon- 
ducting medium are the following: 
E,=VV - I1,+-43M, H,=VV - 07+4303 (9) 
where k,=27/Xo, is the wave number (meters ~) 


Xo is free-space-wavelength (meters) 
€ is free space permittivity (farads/meter). 
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The Hertzian potentials in this case are 


eth e-ikeR 
II,=p, es II,=p} —_— (10) 
Om. Sg 
See. (a1) 
P2 Atrwe, aiid 4a 


The above statement for the fields in an infinite medium must be modified for the case of 
a source in a conducting half space separated by a plane boundary from a nonconducting 
half space. For the modification, the source (transmitting dipole) is located in the conducting 
half space at coordinate position (0, 0, z,) from the boundary. Similarly, the point of obser- 
vation (receiving dipole) is located at coordinate position (r, ¢, 2,) (see fig. 1). The conducting 
half space (sea) is medium 1 (subscript 1), and the nonconducting half space (air) is medium 
2 (subscript 2). 

The Hertzian potentials in the following are normalized for convenience such that |p,|= 
|p.|=|p*t|=|ps|=1. The Hertzian potentials in both media must satisfy the wave equa- 
tion (condition I below) and the radiation condition at infinity (condition IT) below: 


(I) V7, +k701,—0 VOi+ki0{=0, z,>0 (sea) 
V7,+k30,=0 V3 +k 03=0, z,<0 (air) 
(IT) lim 11,—0 jim nt=§ , 2>0 
lim II,=0 lim m0 ee 


where R= yr?+-(2;—2,)?. 

Since the source is in the conducting half space, it is also necessary that II, satisfy the radi- 
ation condition near the source: 

e~ aR 
(IIT) lim 1, =—,—= lim II. 
R-0 f R-0 

In addition to the above three conditions, the electric and magnetic fields must satisfy the 
boundary conditions at the surface of the sea (the plane,z=0). Thus the boundary conditions 
for the components of the Hertzian potentials and their derivatives are the following:' 


























(IV) — jgua=U.2 I,= 22 
= jgUr1=T,2 —jgll}, =z, 
iq MOMs on aT 
Oz Oz 0 Oz 
OM, , M1 Wy» , Oey UY, , OI, OMI» , OU, 
he is Be * Be i Se le * Oe 


where g=o/we=|ki/k3|. In condition (IV), it is assumed that the dipole is oriented either 
vertically in the z-direction, or horizontally in z-direction. (See fig. 1.) Sommerfeld [1949] 
has shown that only a II, component of the Hertzian potential is required to describe the fields 
of a vertical dipole. However, he has shown that both II, and Il, components are required 
to describe the fields of a horizontal dipole (oriented in the z-direction). 





1 In condition (IV) the quantity —jg is an approximation for the complex index of reflection squared, n’. If the frequency or conductivity 
is such that n? A—jg, then one can substitute the exact expression n’=—jg+-«, where « is the relative permittivity of the medium, throughout this 
paper, 


550 





SEA (MEDIUM |) 


VERTICAL DIPOLE 






HORIZONTAL 
OIPOLE 


(x, y,2,)(r,$ ,2,) 
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Figure 1. Distances and coordinate system. 
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3. Hertzian Potential Integrals 


The Hertzian potential can be obtained in integral form for the four basic dipole con- 
figurations: vertical and horizontal, electric and magnetic dipoles. The method used here 
for obtaining the potential integrals was first used by Sommerfeld [1909, 1926]. The method 
is straightforward and available (Stratton, 1941; Sommerfeld 1949], so it will not be presented 
here. A general Hertzian potential for each of the four dipole configurations is obtained 
satisfying conditions (I) through (III). The general solution is then substituted into condition 
(IV) allowing the boundary conditions to determine the arbitrary coefficients of the general 
solution. 

The Hertzian potentials for the observation point in the sea are presented in integral 
form for the four basic dipole pepe The superscripts v and h denote vertical and 
horizontal dipole respectively. For z,>0, 2,>0, 














TN, =A+In 2u=A+In (12) 
and 
W,=A+In Wi=A+Ia (13) 
II?,= cos oe ==]? (14) 
where 
eI Ry 9 ik Ry 
— R, a R, 
Ry=V r+ (2,—2:)? (meters) 
R.=VP4+(z,+2,)? (meters) 
and where 
%0 1 —F(2,4+24) 
In=2 | Fae (15a) 
ro p-F(2, +2), . 
In=2 7 Ps (15b) 
e FCrt2) Jo (ér)édé 
la= -— 1 0 
= 2(— jg i (F+@) (F—jgQ) ta 


= ?—k? (meters) 





=y—k2 — (meters™!). 
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The function # has two parts: exp (—jk,2,)/R,; representing the primary source at the 
position (0, 0, 2,), and exp (—jk,R,)/R, representing a secondary source at the image position 
(0, 0, —z,). In the case of the horizontal electric and the vertical magnetic antennas, the 
secondary source represents the image of the primary source. However, in the case of the 
vertical electric and the horizontal magnetic antennas, the secondary source represents an 
inverse image—that is, an image dipole of the opposite polarity. The primary source radiates 
over the direct path FR, and the secondary source radiates over the reflected or image path FR, 
(see fig. 1). Both paths FR, and F, are entirely within the sea. Consequently, each of the 
Hertzian potentials in (12) and (13) is composed of three components: (a) a primary source 
function, (b) a secondary source function, and (c) an integral. It is shown Jater that the integral 
in every case represents the major contribution to the Hertzian potentials if r>>(z,+2,). 
Therefore it is tacitly assumed that #, composed of the primary and secondary sources, can be 
neglected. Consequently, (12) and (13) can be approximated for 2,>0, 2,>0, r>>(2z,+2,), 
as follows: 


I, =I Tif) =I (16) 
Ih ~1o1 tt ~La. (17) 


I, and J,,, and J, are interdependent. It can be shown that ? for z,> 0, z,> 0, 


ol. : , 
— =Int+jgla (18a) 

8 Sal | 1 far, b’ 
s[ eR; at) ta |= Ie dz, (189) 


The Hertzian potentials for an observation point located in medium 2 are obtained in the 
same manner as were the Hertzian potentials for a point in medium 1. In the case of the 
potentials in medium 2, the function ¥ is not involved, so the potentials can be written di- 
rectly for z,> 0, z-< 0, 





II. =—jgla2 Hi? =], (19) 
and 
r= — JG] 2 73 =—jglaz (20) 
I1?,=— 99 cosy Me Ii} =cos ¢ % e (21) 
where 
°o —F2, + Gz,, R 
lag=2 | eee (22a) 
eo& —F2,+ Gz, . 
In=2 | > “ei gr)édé (22b) 
0 TU 
La) e7F +Gz, J 
? . t J (ér) )édé . 
2=2(—79g— oe 22¢ 
Les \ I9 v{: ( '+G)(G—jgG) —_ ) 


as in (18), the integrals J42, Jy, and J.2 are also interdependent. For z,> 0, 2-<. 0, 


an a (23a ) 
Ola OL ,2 9% 
ta i Dz, 1 d2z,J — 


2 (18) requires that the term (—jg—1) in (15c) be kept in this form. However, for any practical calculation (—jg—1) + —jgsince g>> 1 for the 
conductivity and frequencies of interest in this paper. 
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Now it is possible to write the electric and magnetic fields in medium 1 (sea) as a function 
of integrals 7,,; and J,,; and the fields in medium 2 (air) as functions of J,2 and Jy. Since it 
is the purpose of this paper to determine the fields in the sea and at the surface of the sea, 
attention is focused on the integrals J,; and J,;. If the fields in the sea are determined, the 
boundary conditions can be used conveniently to determine the fields in air at the boundary 
(z, = 0—). Consequently, all fields to be discussed in this paper can be expressed in terms 
of integrals J,, and J,,. It is convenient, however, to replace J,; and J,, with two new integrals 
J, and J, with a corresponding change of variables as follows: 





; =e Ko (pp)ydy 

Za ae Pe ed lew) yay ‘ 
La a af $ L—jgM (24) 

c are as eS ,(pW)vdy or 

= 11,2 |, bi L+M (25) 


where 


M=- G=yV-1 
w 

y=" k 
WwW 
(a3) 7 


=—=—_ ~=—-— 
wie Xo/2r 


w 
fat 
c Xo/27 
e=2,+2, (meters) 
c=speed of light in free space (meters/second) 
\o=free-space wave length (meters). 
The dimensionless distance factors p and ¢ are used to express distances 7 and z in terms of 
(free-space-wavelength)/27. Therefore the electric and magnetic field components in the sea, 
as functions of the integrals 7, and J, are tabulated as follows: 








Electric dipole Magnetic dipole 
Vertical dipole 
, we OT a af a 26 
, Ol, . 7 a ‘ is 
Ei y= 2 —jgla) = == ( Yom igh) (27) 
w dlp 0 OF . 
nme op Ee=j a" 3, (28) 
E,=H,=H,=0 H,=E,=E,=0 (29) 
Horizontal dipole 
2 1 OF, , of " Jae 1 Ole, We 
H,=3 o sin ¢g C dpot *) E,=—j a MSIN g saert se) (30) 
, oF ‘ . Ola , Wa 
H,=3 o COS o(s st E.=— ja is pe ae por or (31) 
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2 3 
H=—% asin g so E.=j% usin ¢ oi (32) 

3 . (Oe Z Ps as 
E,=-3 COS OIG (s3 +1») H,=—% COs ¢ 7G (5 +1.) (33) 
gi atti 8 Pp | 
E.=3 sin ¢ jg & oS +I, ) H.=3 sin ¢ & i +l.) (34) 

b ws? o7/, re w a Ol, = 

ayn cos ¢ Dede A= ai C08 Y Dpd¢ (35) 


As already stated, the electric and magnetic fields in the air at the boundary can be de- 
termined from the fields in the sea and the boundary conditions. The fields are related through 
the boundary conditions as follows, for z,=0: 


Ayn=An =P, ¢, 2 (36a) 
E2=En 1=r,¢ (36b) 
E2=—jgEn. (36¢e) 


The boundary conditions (36) state that the magnetic fields and the tangential electric fields 
in air are equal, immediately adjacent to the boundary, to those in the sea, and that the ver- 
tical electric fields in the air are related to those in the sea by the proportionality constant —7q. 


4. Evaluation of the Integrals 


In the previous section the electric and magnetic fields in the sea and at the boundary 
were written as functions of the integrals J, (24) and J, (25). Apparently these integrals can- 
not be evaluated in closed form. Asymptotic expansions for 7, and J, (as well as J,,, (15a), 
In, (15b), Lag (22a), Zn (22b)) and their appropriate derivatives can be determined by the 
method of critical points [Fredericks, 1953; Baiios, 1953, 1954; Erdelyi, 1956]. 

However, it is possible to reduce the integrals 7, and J, to those evaluated by Sommerfeld 
(1909, 1926, 1949] for the case in which the source and point of observation both lie on the 
boundary. Sommerfeld evaluated these integrals by transforming them into integrals along 
a contour in the complex plane. It is necessary to investigate the contours used to show the 
approximation possible for the case of the source and point of observation in the sea. 

An integrand map is shown in figure 2. There are four branch points and two poles. 
The branch points are associated with both integrals, but J, does not have a pole. The branch 
points arise because of the presence of the square roots in Z and M. They occur where 


L =VP+jg=0 (37a) 
M= Vy—1 =0. (37b) 
B.P. 1 and B.P. 3 are associated with Z, and B.P. 2 and B.P. 4 are associated with M. 
Coordinates of B.P. 1 and B.P. 3 are given by 
B.P.1 00 9 y=¥g/2(—14)) (38a) 
B.P.3 y=yg/2(1—9). (38b) 


It can be seen that their distance from the origin is ¥g. For low frequencies and ‘high con- 
ductivity, as with sea water, this is a rather large value (generally > 10° for communication 
in the sea). 

In the form shown in (37b) it would appear that the branch points 2 and 4 would be on 


the real axis. It is obvious from the limits of the integrals that the path of integration must 
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*% y PLANE 


x aN x/2 | 1/2g 
| ' Figure 2. Integrand map. 





POLE 2 





ae 


lie along the real axis and this could cause certain complications. The reason this branch 
point appears to lie on the real axis is that medium 2 has been assumed to have zero conduc- 
tivity. As a matter of fact, any practical medium, even air, has some conductivity and there- 
fore the effective (complex) permittivity must contain a small imaginary component. 
“ > © 

Thus, the complex permittivity, €, is 

€=€ (1—Jx) (39) 
where x is a very small value associated with the conductivity of air. Utilizing the complex 


permittivity in (37b) we find the branch points are located by 


B.P.2 yx—1+jx/2 


(40) 
B.P. 4 y~1—jx/2. 
The poles occur where (ZL—jgM) is zero. This occurs when 
Ka eke / 
v 1+g° (1 g ee 


For the case of large g (generally > 10°) the location of the poles is given approximately by 
Pole 1 y~—1+ 7/29 
(42) 
Pole 2 y~1— 7/29. 


The integrals as originally stated have a lower limit of zero and an infinite upper limit. 
To aid in evaluating them by contour methods it is convenient to write these integrals in 
forms such that the path of integration lies along the entire real axis. The principal conse- 
quence is the conversion of the Bessel functions of the first kind into Hankel functions of the 
first kind. Thus the integrals become 


7 Lt 
l= | aap HR ova (43) 
20 et : ‘ 
=| fap HP owe. (44) 
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It is therefore possible to close a contour by swinging a semicircle, whose radius is un- 
bounded, from the positive real axis to the negative real axis through the upper half plane. 
Since there are two branch points in the upper half plane, one must form the contour such 
It is, of course, possible to choose con- 





that the two associated branch cuts are not crossed. 
venient branch cuts, within limits. 

Branch line 2 lies along the contour from y=—1+jx/2 to yj for which M is pure 
imaginary. It extends, in essence along the horizontal axis from —1 to 0 and then up along 
the imaginary axis. The contour of integration is shown in figure 3, where it can be seen 
that branch line 2 is the first branch cut encountered in swinging the infinite semicircle from 
positive real axis to negative real axis. The choice of values for the square root in M on the 
two Riemann surfaces separated by branch line 2 is such that for z,>0, the exponent in the 
factor eS has a positive real part in the first quadrant of the y-plane and a negative real 
part in the second quadrant. (This aids in obtaining convergence of the integrals in the non- 
conducting medium (22) along branch line 1.) One observes as Sommerfeld did that in this 
case (high conductivity) Pole 1 and B.P. 2 are very close together so that the contours for 
J, must pass around the branch point and the pole with the same loop. This problem has been 
treated in great detail by numerous authors and will not be elaborated here. 

Branch line 1 extends from y=+/g/2 (—1+ 7) to y>jo. Along this line Z has a large 
imaginary component and a small real component. The choice of Riemann surfaces is such 
that the real part of the exponent in the factor e~“ is negative on the left side of the branch 
line, because this aids in exponential reduction of the integrand along branch line 2. 

It can be shown readily that the integral along the infinite semicircle is zero. Following 
the terminology used by other writers in the past, the integral along branch line 2 is called 
P+Qz2 (except for J, where it is just Q.). The integral along branch line 1 is called Q,. Thus 
for each of the integrals we may state that the value is given by 


I=P+@:+0; (45) 

where P is, of course, zero for Jp. 
For a highly conducting medium, except very close to the source dipole, the contribution 
to the integral along branch line 1 is negligible. One can validate this statement by examining 
the asymptotic expansion of the zero order Hankel function of the first kind (Stratton, 1940] 


which is 


¢ 


| 9 
HY (oY) ~a) et". (46) 


™p 





It is seen that the Hankel function vanishes exponentially as y increases along branch line 1. 
Since the smallest imaginary part of y associated with branch line 1 is /g/2, the contribution 
Q, is justifiably neglected. Consequently, it can be assumed that the significant contribution 
to the integral comes from the integration along branch line 2; that is 


I~P+Qo. (47) 
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Figure 3. Contour of integration. 














The most important approximation made in this treatment is the one allowing the removal 
of the exponential e-* from the integrand. To see that this is possible, it is necessary to express 
J, in the form of either (43) or (44), as the sum of two integrals for integration along branch 
line 2. The limits of the first integral are B.P. 2 and a point y on branch line 2; and the limits of 
the second integral are yy and jo. The point yo is chosen such that 1<<|W)?<<y. 
For example, using g=10", one might choose yo=j10*. Consequently, one can approximate 
L= P+ jg~ vjg over the first integral (B.P. 2, Yo) with negligible error introduced in the 
evaluation of the first integral. Over the second integral (Y,jo), the inequality 





e~ vigs <le-™| (48) 


permits the following inequality to be written: 


ee ( emt 2 (1) = 9 
(a E o (op) yd S | Iran H 0 (pb) ydy B (49) 
BLL. | B.L.2 | 
(Yo weet (Wo<v<j «0 ) | 


— 7g for I, 
where a= i 
1 for I, 


Using the approximation 1/=1y—1 ~y, over (Yo, j ©), and (46) with 10-2<<p< 102, |p| > |Wp| > 10°, 
the following upper bound is written for the inequality (49): 


TD (I? e-Ugieb —IADq Wy 3/2qy| | | er Vistt+ior ee apes 
Bs y= | : : a Sar Sa 1" v < | ae » ected (—Jp*¥y+ 2p +2) (50) 
| VamoJy, Jg+ la )y Vi 59g"? 


It is shown later that a contribution whose magnitude is less than the upper bound (50) is 
negligible. Consequently, for the following approximation in J: 

L=vjg (51) 
one observes that the error is negligible in the first integral and that the contribution of the 
second integral, where the approximation is relatively poor, is also negligible. Hence the 
approximation is justified. Using (51), the integrals 7, and J, for the conductivity and fre- 
quency range considered in this treatment are written as follows: 


~1-» vase (HS (ody Re’ 

=P , = Qe‘ <=) g cn ft 52 
I,~P+Q: ~2¢ vel2 { im" 7 (52) 
B.L 


9 


and 


omy @ £00) 
T= Q, Qe —P vals [Beene (53) 


B. L.2 





This approximation is important because it places outside the integral the factor associated 
with the distances of the dipole and of the point of observation from the boundary. Thus 
the remaining integral is simply that for observation at the surface of fields due to a source 
also located at the surface. This was the problem Sommerfeld [1909] treated for the vertical 
dipole, and the horizontal dipole problem has also been treated since then by many writers. 
Thus the problem for fields due to a dipole in a conducting half space has been reduced to the known 
problem for fields due to a dipole at an interface. 

Both convergent and asymptotic series were developed by Sommerfeld and others for 
this integral. The expression for 7, was found [Moore, 1951] for 0<p<1, to be 
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“ ee Vint lp,1 jp wr a 2 ae 
oie aillae”” in ae 213g Utse+5 oa 1+ Jot |... 


1 = 
45 EHP (0) +5 Go HO Dt} 6s 
“VIG 

This may be expressed as 
j eo te—vi8k 
I, = j2 a. r. (55) 





For distances such that 10°°<p<1, conductivity o~4 mhos/m, and frequency {<50 ke/s, 
then 7~1. A more convenient form for J, given by Sommerfeld [1949] is the following for 


>t: 

.. @ FP viot an : 

[,=j2 —— "[i+d N 2rn es -—nenv | erin | (56) 
— gp 2 0 


An expression for J, similar to that for 7, (54) was also obtained by Moore. However, J, 
can be expressed in closed form [Wait, 1952] as follows: 
e—e— viot 


I,=2 Go—De [1+jo—(1+-Vjgp)e~ vee], (57) 





For the practical purpose of submarine communications; that is, ¢~4 mhos/m, f<.50 ke/s, 
it is sufficient to approximate J, (54) or (56) and J, (57) as follows: 


e—je— viot 








IL -=2 —» 10°°*<p<10?, p>>t (58) 
and 
ene Vi0k . . 
L=—j2— 4 (1+je), 10-°*<pS1, p>>b (59a) 
gp 
or 
e—le— vias 
1y=2 (1-2) 1<p<10°, p>>¢. (59b) 
- 


Now that the integrals 7, and J, have been evaluated, the derivatives of J, and J, involved 
in the field components ((26) through (36)) must be determined. It can be shown [Moore 1951] 
that the power series (54) can be differentiated, allowing the necessary derivatives to the same 
order as of (58) to exist. J, (57) can be differentiated; consequently, the necessary derivatives 
of (59) also exist. Returning to the statement that # can be neglected in (16) and (17), one 
can use (16), (17), (24), (25), (58), and (59) to verify this approximation. Similarly, the 
approximation (51) can be justified by comparing the magnitude of the upper bound (50) to 
(58) or (59). For the case g=10", yo=710* and 10° *<p<10, the error is less than 1 percent. 


5. Fields in a Conducting Half Space 


After performing the indicated differentiations, substituting numerical values for per- 
mittivity, permeability, and velocity of propagation in free space, and inserting dipole moments, 
the resulting far field expressions (26) through (35) for p>1, are 


7/2 2/6 
Ev ~—¥j1.76X10-# 2 ee r (1-4) (60) 
3 —2/8 " 
yy ~—1.57<19-8 SE T (1-2 (61) 
Je 2/8 6 
E%™ ~1.96X10-* oe (1) (62) 
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—_ ul - 
Hi ~V—j1.76X10-” = ani (i— d 
o} a p 


bo 
wee” 


ee 3 o "Tle ** +T jT). 
H'™ ~=—¥7j1.05X107? aw” a1 aa |? 9 + sing 


Ht = — 7 5.26X107*4 ¥ a T 7; 2T — COS ¢ 
Er =w—5. 92x19-7 Te ae Mir (2T— )—% J eos ¢ 
EM = j1.18X10-* 2 = [es sing 

op? 2 2» 


2/8 . ; 
Em ~— yj6.64x10-” 2 SNe 7 (1-24) sin g 
o'/*p . Pp Pp 





; pi 4 wT § Ne-2* : 
E ie —— 9 1 33X10 - aoe za ie (1-4) cos ¢ 


H™ = 1. 1719-6 2! SNe 
p” 


seh x (1-2) cos 
bg 3 S] a A \ i 
H'™ ~ —5.83X10- a a-i4 ) sin g. 
p | + 
Similarly, the resulting near field expressions (p<1) are 
<10- oes all 


Ex ~\—j1.76) a 0) 


H% ~j1.57 107” © file” tex (1+ jp) 
op” 


a : 
Em ~—5.90X 107 0 N (1+i0—% ) 





— , w/*7S Ne *6 , 2) 
Hm ~— ¥—75.26X at i a a ( i ne 
V¥—j 5.26 X 10 oiigs 1+ Jp = 


H** x V¥—j1.05X we ee : (1+ Jp) sing 


o, wo *Ile— 7/8 
H™ x — /—75.26X10- a OT a a (1+je—p?) cos ¢ 


E* =5.92X10-” ae (1+-Jje—p?) cos 9 
Et ~1.18X107 a (1+ Jp) sin ¢ 


- 27S Ne? ’ 
EJ" =j 6.64X10-” ee (1+jp—p?) sin ¢ 


> hes a al _o9 wT. S NeW */ 
Em ~—yj 1.33X10-? 29 wo" 3 = ge) cos ¢ 
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(63) 


(64) 


(66) 


(67) 


(68) 


(69) 


(70) 


(71) 


(74) 


(75) 


(76) 


(78) 


(79) 


(80) 


(81) 











sei 2/ 


mM =147>40- res cos ¢ (82) 





HH,’ ~5.84X10 (1+jo—p?) sin ¢. (83) 


og wTS) 2 ibd 

p® 
Here the superscripts ve, vm, he, and hm refer respectively to vertical-electric, vertical-mag- 
netic, horizontal-electric, and horizontal-magnetic dipoles. TJ is determined from (54) and 
(55). 

The fields are all expressed in terms of p, the radial distance in units of (free-space-wave- 
length)/27. Since p is directly proportional to both w and r, the frequency dependence of the 
result is not what it would be if p were independent of frequency. For example, at values of 
10-?<p<1, several of the fields are independent of frequency except for the depth attenuation 
factor: E,”", E.”, ,", and H™. 

The z-components of the fields in the conducting halfspace (sea) are all small (zero to a 
first-order approximation for the vertical electric dipole) compared with the horizontal compo- 
nents, and therefore they have not been included. On the other hand, the z-component of the 
electric field in the air, for the horizontal dipoles, is the predominant component. The ex- 
pressions for the electric field vertical component in the nonconducting region at the surface 
can be obtained by applying the boundary condition (36c) for p>1; they are 


EX = — \-j 1.98107?! pal. fatal — T ( 1— 1) COS ¢ (84) 
p 
3 Tp —24/8 % 
E hm ~2 22x 19-21 2S Ne 3 ” 1—7) sin ¢. (85) 
p p 


Similarly, the results for p<1 are 


EX ~—yj 14sxi190°> er (1+ Jp) cos ¢ D (86) 
a’*p” 
Ets = — j2.22X10-** = — Sa Lead ee sin ¢. (87) 


These fields are given only at the surface as no effort has been made to evaluate the necessary 
field integrals at any distance away from the surface in the nonconducting medium. The 
horizontal electric antenna is the most practical for use in submarine communications. As an 
illustration of the manner in which the electric and magnetic fields vary as a function of distance 
from the source; the fields for the horizontal electric antenna have been plotted as a function 
of p in figures 4, 5, and 6, for typical values of frequency, conductivity, and practical antenna 
current, length, and depth. 


6. Discussion and Conclusions 


The form taken by the field equations indicates the mode of propagation that prevails as 
long as the distance from the boundary to both point of observation and source dipole is small 
compared with the horizontal separation. Each of the field expressions may be considered to 
consist of 3 parts: a multiplying factor which includes the dipole strength and parameters such 
as frequency and conductivity of the medium; an exponential attenuation factor whose expo- 
nent is the sum of the distance from the dipole to the surface and from the surface to the point 
of observation, expressed in units of skin depth; and a factor associated with variation in the 
radial (horizontal) direction. The latter factor is the same as that for surface fields of surface 
dipoles. 

Thus it appears that propagation occurs as indicated in figure 7. The wave starts at the 
dipole, proceeds by the shortest path to the surface (the path of minimum attenuation), is 
refracted at the surface, travels along the surface as a wave in air, and comes to the point of 
observation by the path of least attenuation (straight down). 
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Figure 4+. Variation of the maximum absolute value of the | Fiaure 5. Variation of the maximum absolute value of the 
magnetic field components in the sea for a submerged hori- electric field components in the sea for a submerged horizontal 
zontal electric dipoleywith radial distance in units of free- dipole with radial distance in units of free-space-wavelength/2r. 
space-wavelength/2r. 





It is possible to consider that an equivalent source on the surface sets up the wave in the 
nonconducting medium. This wave then travels as would be expected for such a source. 
Since the sea is not a perfect conductor, there is a slight tilt to the wave front, as indicated 
by the presence of both radial and axial components of the electric field; the radial (hori- 
zontal) component is much smaller than the axial (vertical) component. The Poynting vec- 
tor for such a wave is nearly horizontal, but tilted slightly toward the conducting medium. 
The component of the Poynting vector associated with that part of the wave traveling into the 
conducting medium is small but finite. This wave which travels straight down is that which 
is seen by the observer in the conducting medium. 

A similar situation prevails for waves in an imperfectly conducting waveguide. It is 
customary to calculate the losses in the wall of the waveguide by calculating the power flowing 
into the wall per unit area. In the waveguide this represents a leak of power from the desired 
wave. In this situation, where the observer is actually located in the conductor, the desired 
signal is the same as the “leak” for the waveguide. 

Because of this mode of propagation, it is possible for radiation starting with the dipole 
in a conducting half space to be observed at much greater distances than would be possible if 
the wave were generated in an infinite homogeneous conducting medium. Thus, the exponen- 
tial attenuation of 55 db per wavelength applies in the case of the conducting half space only 
to that part of the path from the dipole to the surface and from the surface to the point of 
observation. The rest of the path undergoes the same sort of attenuation that a wave gener- 
ated by a dipole in air would encounter. 
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horizontal electric dipole with radial distance in units of free- 
space-wavelength/2r. 
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Although this development has applied to the case of a conducting half space and the 
expressions given are those associated with the approximation of the earth’s surface by a 
plane, the fact that it is possible to separate out the horizontal attenuation of the wave in air 
means that it is also possible to apply the transmission loss that would be developed using 
spherical, rather than plane, earth geometry. In fact, it is also possible to use the attenuation 
calculated when the ionosphere is taken into account. Thus the results obtained here are 
considerably more general, within the limits of the approximations used, than would be indicated 
by the assumption of a plane earth. 

It is interesting to note the magnitude of the attenuation of the wave in going through 
the conducting medium for the special case of sea water. This is shown for several depths as 
a function of frequency in figure 8. The depth indication is in meters. It is obvious that if 
the attenuation is to be kept low, the frequency must be kept low also. 

It has been shown here that the fields created by a dipole in a conducting medium, which 
is near a plane boundary with a nonconducting medium, may be expressed in a form that may 
be considered as propagation from the dipole to the interface, travel through the nonconducting 
medium along the boundary, and propagation back into the conducting medium to the point 
of observation. Expressions have been presented for the fields of both vertical and horizon- 
tal, electric and magnetic dipoles. The fields of the horizontal dipoles are considerably stronger 
than those of the vertical dipoles, as would be expected because of the vertical nulls in the 
radiation from the vertical dipoles. 
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Measured radio noise values are compared with the corresponding International Radio 
Consultative Committee [C.C.I.R., 1957] predicted values at four noise measuring stations. 


Five frequencies between 0.013 and 10.0 Me/s are considered. 


The stations selected for this 


study include Balboa, Panama, near two major radio noise centers, and Byrd Station, 


Antarctic, remote from atmospheric radio noise sources. 


It is found that the predicted and 


measured noise levels are in good agreement except at some places and times, where large 


discrepancies occur. 


are based on extrapolations of data measured at other stations. 


ments are discussed. 


1. Introduction 


Because the success or failure of communications 
by radio waves depends upon the signal-to-noise 
ratio at the receiver, it has become necessary to 
obtain some knowledge of the noise levels to be 
expected. Atmospherically-generated radio noise 
levels vary with time of day, season, and geographi- 
ral location. Prediction of tbe variations on a 
worldwide basis have been published in Radio 
Propagation Unit (R.P.U.) Technical Report No. 5 
[1945, 1947] and National Bureau of Standards 
Circular 462 [1948] in terms of the minimum field 
strengths required to maintain intelligible voice 
communications ninety percent of the time in the 
presence of noise. Subsequently, the prediction 
curves were revised in NBS Circular 557 [Crichlow, 
1955] to show expected median levels of radio noise 
power during four-hour time blocks for each season. 
As more recent data accumulated it became apparent 
that large discrepancies existed between the measured 
and predicted values at some locations and times, so 
that an entirely new set of predictions was prepared 
and published in International Radio Consultative 
Committee (C.C.I.R.) Report No. 65 [1957]. The 
C.C.I.R. predictions closely follow those in Circular 
557, but are based on a wider selection of compatible 
data, and should thus be more useful than the 
previous publications. 

Cooperative measurements of received noise power 
at 16 widely separated geographical locations (fig. 1) 
are currently being conducted on a routine basis 
under the direction of the National Bureau of Stand- 
ards Central Radio Propagation Laboratory. Since 
the issuance of the C.C.I.R. predictions, several 
additional noise measuring stations have been estab- 
lished and approximately 5 years of data have been 
collected. Thus it is now possible to check the ac- 
curacy of the predictions during a time of declining 
sunspot activity (1960), and at geographical loca- 

1 This work was supported by the U.S. Air Force Cambridge Research Labora- 
tories under Contract No. AF19(604)-4092. 
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FieurE 1. Map showing location of noise measuring stations 
(reproduced from NBS Tech. Note 18). 

tions where the noise level could only be interpolated 
when the predictions were established. 

In this paper comparisons are made of the pre- 
dicted and measured noise values at the four sta- 
tions listed in table 1 for the frequencies (in Me/s) 
shown. The letter X indicates the frequencies 
analyzed at each station. These particular stations 
were selected because their noiso grades range from 
high to very low. That is, the noise grade at Balboa 
is high because it is near two major noise centers in 
Central and South America, while that at Byrd Sta- 


TABLE 1 
Frequency (Me/s) 





Station and location__.-......--- 10.013 lo.051 0.246 (0.495 (0.545 |2.5 [5.0 |10.0 
Balboa, Panama, 9° N, 80° W___|_..-_.| X xX xX 7 Poe ae 
Byrd Sta., Antarctica, 80° S, | | | | 
Uo EOE ee ee eae ox ba aS 
Enkoping, Sweden, 60° N, 17° E_|_- x a rs x | x 
of ae xX xX} X/|xX 


Kekaha, Hawaii, 22° N, 160° W- -| x 
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tion is very low because it is remote from all at- 
mospheric noise sources. The noise grades at Kekaha 
and Enkoping are between these two extremes; that 
at Kekaha being greater than at Enkoping. 


2. Noise Level Predictions 


The C.C.I.R. predictions are presented as world- 
wide noise level contour maps for 1.0 Me/s in 4-hr 
time blocks for each season. Curves are given to 
determine the noise levels at other frequencies be- 
tween 10 ke/s and about 30 Me/s. The magnitude 
of the noise is in terms of a basic parameter repre- 
senting the median value of atmospheric noise power 
passed by a narrow bandwidth, expressed for a 4- 
hour time block. The seasons of the year for the 
present purpose are defined as in table 2 

The basic parameter, F, [Crichlow, 1955], used to 
describe the noise is defined as the noise power avail- 
able from a short, vertical, grounded, loss-free 
antenna in decibels above kTB, where kTB is the 
terminal noise power in watts, available from a pas- 
sive resistance (dummy antenna). Here, 

k =Boltzmann's constant= 1.38 10~** joules per 

degree Kelvin. 

T=absolute room temperature=288° K. 

B=effective noise bandwidth in cycles per second 


(c/s). 


TABLE 2. Seasons of the Year (in accordance with C.C.I.R. 
Report No. 65 and NBS Technical Note 18) 


Season 


Months 
Northern Southern 
Hemisphere | Hemisphere 
December, January, February A or Con Winter. Summer. 
Re, OT. REO. nc eno oe ec ead eS ao Autumn. 
June, July, August_-_-- TEL Winter. 


Se ptember, October, November- Autumn... Spring. 


2.1. Prediction Curves for Each Station 


The expected atmospheric noise level at each of the 
four noise measuring stations listed in table 1 are 
shown in figures 2 to 5 by solid lines for each season 
of the year ‘and each frequency. The corresponding 
measured noise values are shown in the figures by 
dashed lines, and the estimated manmade noise levels 
at the stations are depicted as broken lines. (All 
times referred to in this paper are local times.) 
More about the measured values is given in the next 
section. 

The expected atmospheric noise levels shown in 
figures 2 to 5 were prepared from the noise grade maps 
given in the C.C.I.R. report. Because the noise 
maps are for 1.0 Me/s only, the predicted values 
of F, for the frequencies considered here were ex- 
trapolated using the curves of noise power as a func- 
tion of frequency given in figures 21 and 22 of that 
report. In a few cases, notably at Marie Byrd 
Station and at Kekaha for some time blocks, the 
interpolation between noise isopleths on the maps is 
rather broad, but the chosen values should be correct 
within one or 2 db. 














The estimated manmade noise levels shown in 
figures 2 to 5 were also obtained from the C.C.I.R. 
report. The values represent what should be ex- 
pected at receiving sites located several miles from 
populated areas and at least a quarter of a mile 
from elevated power lines and electrical machinery. 


3. Noise Level Measurements 
3.1. Seasonal Time-Block Measured Values 


The noise power is recorded for 15 min per hr for 
each frequency, and each sampling is assumed to 
represent the noise level for the full hour. From 
these data the month-hour medians, F;,,,, are ob- 
tained in terms of decibels above kTB and are com- 
piled in the NBS Technical Note 18 Series [Crichlow, 
1959 a and b}. Also given are the seasonal 4-hr time 
block values, averaged over each 4-hr period for the 
3 months of the season in accordance with C.C.I.R. 
recommendations. The tabulations include the up- 
per and lower decile values of F, to give an indica- 
tion of the extent of the day to day variations from 
the median, F,,,. The appropriate seasonal time- 
block noise power averages as reported by Crichlow 
[1959b] are plotted as dashed lines in figures 2 to 5. 

Since the noise-measuring equipment approxi- 
mates, insofar as practic able, the idealization of the 
system for which the predictions were generated, 
and corrections have been made in the measured 
data for antenna losses and the like, the measured 
and predicted values of noise power given in figures 
2 to 5 can be compared directly. 


3.2. Comparison of Measured and Predicted Sea- 
sonal Time-Block Values 


As is evident from the curves of figures 2 to 5, a 
great number of discrepancies exist. between the 
predicted and measured values of noise power. In 
several cases there is very good agreement, but in 
others the predicted noise level is as much as 40 db 
lower or 24 db higher than that actually measured. 
However, it should be noted that in most of the 
disagreements the predicted noise level is on the low 
side, and a predicted value 24 db too high occurs 
only once (winter 1600-2000 hour time block at 
Byrd Station on 0.051 Me/s) in the data analyzed 
here. Furthermore, most of the too-high predicted 
values occur at Byrd Station and Enkoping, where 
the prediction curves are based on scanty data or 
extrapolation techniques. 

Measured values that are much higher than the 
corresponding predicted values can be attributed to 
a number of causes. One of the most important 
causes is manmade noise near the receiving station. 
This can arise from automobile ignition systems, 
power lines, and electrical machinery, which is prop- 
agated primarily over power lines or by ground 
waves. Since manmade noise is random it is pres- 
ent on all the measured frequencies. Another im- 
portant contribution to the measured noise level 
might arise from local station interference, and this 
would show up on only one frequency. Atmospher- 
ics generated by local thunderstorm activity would 
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Predicted and measured noise power at Balboa, Panama. 


Prediction curves scaled from C.C.I.R. Report No. 65. 
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Figure 3. Predicted and measured noise power at Byrd Station, Antarctica. 


Prediction curves scaled from C.C.I.R. Report No. 65. 


greatly increase the received noise power, but this | made, cosmic, and atmospheric; which happened to 
has been taken into account in the prediction curves | be incident on the antenna at the time of measure- 
from a knowledge of worldwide thunderstorm activ- | ment. Data which obviously contained manmade 
ity. Thus, unless the locale of a receiving station | noise, interfering signals or equipment trouble have 
suffered a higher-than-usual number of thunder- | been omitted from the compilations, but these effects 
storms spaced over a 3-month period, local atmos- | have not been eliminated in all cases [Crichlow, 1959 
pheric noise would not be the cause of a measured | b]. Thus, the high noise level measured at Balboa on 
noise value much higher than that predicted in a | 0.495 Mc/s in the spring and summer (fig. 2) probably 
seasonal time block. contains a strong contribution from a local interfering 

The prediction curves were issued for atmospheric | signal. Were it manmade noise or local atmospherics, 


noise levels primarily, while the noise data [Crichlow, | the corresponding data on 2.5, 5.0, and 10.0 Me/s 





1959 a and b] include radio noise of all types; man- | would show its effects. The higher-than-predicted 
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FIGURE 4. 


Predicted and measured noise power at Enkoping, Sweden. 


Prediction curves scaled from C.C.I.R. Report No. 65, 


noise values measured at Balboa on 5.0 Mc/s during 
the night hours (2000-0400) of spring, summer, and 
autumn also apparently contain an interfering signal. 
In this case, however, the signal appears to have been 
ionospherically propagated from a distant trans- 
mitter, because it disappears during the daytime 
when ionospheric absorption is effective. 

The data from Byrd Station (fig. 3) are of particu- 
lar interest because this station is located in a very 
low noise region where no previous noise measur- 
ments had ever been made. The C.C.I.R. report 


602247—61——-3 





suggests caution in the use of predicted values which 
fall below the estimated manmade noise levels for a 
quiet receiving site, and this is just the case for the 
frequencies of 0.545, 2.5, and 5.0 Me/s at Byrd 
Station. Since thunderstorm activity is practically 
nonexistant at Byrd Station, the atmospheric radio 
noise must be propagated in from distant sources, 
such as the major noise centers in equatorial Africa, 
South America, and the East Indies. Other contri- 
butions to the measured noise level would, of course, 
arise from local manmade sources and interfering 
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Figure 5. 


Predicted and measured noise power at Kekaha, Hawaii. 


Prediction curves scaled from C.C.I.R. Report No. 65. 


signals. In addition, as pointed out by Crichlow 
[1959 a], the methods used for determining antenna 
losses at Byrd have not been completely accurate, so 
this may contribute to the observed discrepancies 
between the predicted and measured noise values in 
figure 3. At times during winter a large amount of 
precipitation noise can be generated at the receiving 
antenna by snow being blown against the antenna. 





Two additional noise sources may be of importance 
in polar regions, such as at Byrd Station. The first 
consists of aurorally generated noise which at times 
can be as much as 15 db above atmospheric noise on 
18 Me/s [Egan and Peterson, 1960] and significantly 
greater than atmospheric noise on 5 ke/s (Elis, 1959]. 
The second, of importance on frequencies above 
about 2 Me/s, is cosmic noise which can penetrate 
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the ionosphere near the poles where the critical fre- 
quency is low in winter, suffer ground reflection, and 
subsequently be propagated below the ionosphere by 
conventional means (Ellis, 1958]. 

To some extent the various types of noise sources 
which contribute to the measured values can be 
identified. For example, the presence of ionospheri- 
cally propagated atmospheric noise in the Byrd Sta- 
tion data is indicated by the following two facts: 

(1) A small seasonal variation can be seen (fig. 3) 
in the measured data on 0.051 and 0.545 Mc/s, where 
the minimum occurs during winter. This minimum 
probably corresponds to the absence of the D and FE 
reflecting layers during the sunless Antarctic winter. 

(2) The noise levels are generally higher in winter 
on the 5.0 and 10.0 Me/s frequencies, complementing 
the behavior of the low frequencies in fact (1). 
That is, propagation of the high frequencies is ex- 
pected to be better in winter because of the absence 
of the D region which is an absorber for these 
frequencies. 

In considering the data for 0.051 and 0.545 Me/s 
alone, one might be tempted to infer a predominance 
of manmade noise, where the seasonal minimum 
coincides with the relative inactivity of the station 
during winter. However, if this were true the noise 
level of the high frequencies would also be less during 
winter. On the other hand, that there is substantial 
manmade noise present is evidenced by an almost 
complete lack of seasonal and diurnal variation in 
the 2.5 Me/s data. It should be noted that the 
measured values for this frequency fall about 20 db 
below the estimated level of manmade noise. 

The measured noise values at Enkoping (fig. 4) 
on 0.545 and 2.5 Me/s are substantially higher than 
predicted for all seasons during the daytime, which 
would indicate either a large contribution from man- 
made noise or interference, or an error of 20 to 30 db 
in the prediction. On 0.051 Me/s the predicted 
values for the winter afternoon and evening time 
blocks are 21 and 18 db too high, respectively, and 
for the same time blocks durmg summer the predicted 
values are 16 and 5 db too low. 

In figure 5 it can be seen that generally good 
agreement obtains between the predicted and meas- 
ured noise values at Kekaha, Hawaii. During the 
day on 2.5 Me/s the predicted noise level is on 
the order of 20 db below the measured level, but 
both are below the estimated manmade level and 
must be viewed with reservation. The higher-than- 
expected values measured on 5.0 Me/s at Kekaha 
probably contain a large amount of station interfer- 
ence, as suggested by Crichlow [1959 a]. 


3.3. Month-Hour Median Measured Values 


As an illustration of the diurnal noise patterns 
occurring seasonally in various parts of the world, 
month-hour medians of Fa measured on 5.0 Me/s 
are plotted in figure 6 for selected months. The 
upper decile curves shown in this figure represent 
the noise value exceeded 10 percent of the time, 
while the lower decile value was exceeded 90 percent 
of the time. 





Comparison of the month-hour median curves in 
figure 6 with the corresponding time-block curves 
in figures 2 to 5 reveal variations in the median of 
only 2 or 3 db within a 4-hr time block, except during 
sunrise and sunset periods. Month-hour median 
plots for the other frequencies compare similarly. 
With this in mind, the diurnal features of the different 
frequencies can be discussed on the basis of the.curves 
in figures 2 to 5. 

As expected, the diurnal variation of the fre- 
quencies at 0.495 Mc/s and higher show the effect of 
daytime ionospheric absorption, while the 'lower 
frequencies at 0.013 and 0.051 Me/s exhibit little daily 
variation. The situation at Byrd Station, and to a 
lesser extent at Enkoping, is somewhat different. 
Smaller diurnal variations occur during summer and 
winter on the higher frequencies because the long 
periods of daylight during summer allows absorp- 
tion 24 hr a day, which is completely absent during 
the sunless winter. The days of the transition 
seasons, fall and spring, have hours of daylight as 
well as darkness, so the diurnal variation at Byrd 
Station during these months is greater. 

At the equator the times of sunrise and sunset do 
not change with season, so that the spread of the 
measured noise value in the 0400-0800 and 1600-— 
2000-hr time blocks should be about the same all 
year. At Kekaha (22° N) the sun rises at about 
0530 in summer and at about 0630 in winter, so this 
difference is not too important. However, at 
Enkoping (60° N) the sun rises at 0300 in summer 
(June) and not until 0900 in winter (December). 
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Thus, during summer at Enkoping large variations 
occur in the 0000-0400 hour time block (fig. 6), but 
during winter the variations due to sunrise condi- 
tions occur in the 0800-1200-hr time block. From 
this it is evident that the predicted values for time 
blocks containing sunrise and sunset can give only a 
rough indication of the median noise to be expected 
over the entire periods, and are less representative 
for individual hours than those for other time blocks. 


3.4. Measured Variations in Medians 


Within a 4-hr time block large day-to-day varia- 
tions from the median are observed in the measured 
noise power. The noise power variations are ex- 
pressed in terms of the ratio (in db) of upper decile 
to median, D,, and lower decile to median, D,, and 
are tabulated by Crichlow [1959 b] along with the 
monthly median values. An attempt was made in 
the C.C.I.R. report to predict the magnitude of D, 
and D, in each 4-hr time block as a function of 
frequency, but without regard to season or geo- 
graphical location. 

Thus, to compare the measured and predicted 
values of D,, and D, the measured data were averaged 
over a l-yr period (September 1959 to August 1960) 
in each 4-hr time block. For illustration, the average 
variations measured in the 1200—-1600-hr time block 
at each of the 4 stations on the different frequencies 
are plotted in figure 7 as a function of frequency. 
The predicted values of D, and D, as a function of 
frequency for the same time block are shown as solid 
curves. 

In figure 7 it can be seen that the measured varia- 
tions are as much as 13 db greater than predicted 
(as at Balboa on 2.5 Me/s) or 13 db Jess than pre- 
dicted (as at Byrd Station on 0.545 Me/s). At 
Kekaha, the predicted and measured values of D,, and 
D, are generally in good agreement. Differences 
between the predicted and measured values of D, 
and DP, for individual months are greater than 13 db 
in some cases. Generally, it can be said that the 
measured variations at high latitude stations are less 
than predicted while those at stations near the equa- 
tor are greater than predicted. 


3.5. Noise Power as Function of Frequency 


Figure 8 shows the decrease in noise power with 
increasing frequency. The plotted points represent 
the measured noise power in the 0000-0400 time 
block averaged over 1 yr. The solid curve for each 
station was obtained from figure 22 in the C.C.I.R. 
report and indicates the expected variation of noise 
power with frequency. 

The measured noise power variation with fre- 
quency as shown in figure 8 corresponds rather well 
with that predicted at all stations. The largest 
discrepancy occurs at Balboa on 0.013 Me/s, where 
the measured value is about 17 db higher than that 
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predicted. 
in only 10 percent of the measured points in figure 8. 
However, during individual months greater de- 
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partures appear, which should be taken into account 
in any future prediction curves. The good agree- 
ment between predicted and measured values in 
figure 8 indicates that the measured data considered 
here primarily represent noise of atmospheric origin. 


3.6. Ratio of Measured to Predicted Values 


The percentage distribution of the deviations of 
predicted from measured values are shown in figures 
9a and 9b. The ratios of measured to predicted 
values are divided into the class intervals +0 to 10 
db, +10 to 20 db, and +20 db or greater. The per- 
centage distributions are shown in two ways. In 
figure 9a, the observed deviations for 1 yr counted 
on all frequencies at each station are given to illus- 
trate the reliability of the predictions at different 
locations. In figure 9b, the observed deviations for 
the same year counted at all stations on each fre- 
quency are given to illustrate the reliability of the 
predictions for different frequencies. 

From the percentage distributions of the ratios of 
observed to predicted noise levels given in figure 9a, 
it appears that the predictions are most reliable at 
Kekaha and Balboa, where over 80 percent of the 
observed deviations are 10 db or less. At Enkoping 
and Byrd Station less than 70 percent of the observed 
deviations are 10 db or less. At Byrd Station, 18.4 
percent of the observed deviations were greater than 
+20 db, and nearly all of the observations were on 
the plus side. This strongly suggests local inter- 
ference, which also adversely affects the percentage 
distribution of 0.545 Me/s given in figure 9b. 

The percentage distributions given in figure 9b 
indicate that the predictions are most reliable for 
0.013 and 10.0 Me/s, where 100 and 90 percent of the 
observed deviations were 10 db or less, respectively. 
The large percentage of deviations greater than 20 
db (33%) observed on 0.545 Me/s is due to the Byrd 
Station data. With those data excluded, the per- 
centage deviation of the three class intervals for 
0.545 Me/s become 63, 21, and 16 percent, respec- 
tively. 

The distributions for 2.5 and 5.0 Me/s ean be 
compared with a similar reliability test of the revised 
R.P.U. predictions as reported in Circular 462 (NBS, 
1948). The comparison is not strict, however, due 
to the difference in sampling used here. The present 
distributions are based on the measurements of four 
stations, while those in Circular 462 were based on 17 
siations. Thus, it is not unexpected that the relia- 
bility indicated in figure 9b on both frequencies is 
less than in Circular 462, where the percentage of 
observations for the three 10 db class intervals (0 to 
30 db) were: 75, 20, and 5 percent for 2.5 Mc/s; and 
90, 8, and 2 percent for 5.0 Me/s. 


4. Summary 


Although good agreement is generally found 
between the predicted and measured noise power 
values, large differences still exist at some places and 
times, At places such as Kekaha, where the predic- 
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tions are based on previous noise measurements, the 
differences are small, but at locations where the pre- 
dictions are based on interpolated or extrapolated 
data (as at Byrd Station), large differences occur. 
The differences between predicted and measured 
values range from +40 to —24 db. 

In time blocks containing sunrise or sunset, the 
predicted values give only an indication of the median 
noise power to be expected over the entire periods 
and are less representative of individual hours than 
those for other time blocks. 

The measured variations from the median, aver- 
aged over a 1 year period, differ from the predicted 
variations by as much as +13 db. In individual 
months larger discrepancies than this occur. 

The noise power data considered in this paper 
follow the expected frequency dependence very well. 
Discrepancies greater than 10 db occur in only 10 
percent of the cases, and the greatest departure is 17 
db at Balboa on 0.013 Me/s. 

Of four stations investigated, it appears that the 
C.C.I.R. predictions are most reliable for Kekaha 
and Balboa, where over 80 percent of the ratio of 
measured to predicted values were 10 db or less, while 
less than 70 percent of the observations at Enkoping 
and Byrd Station were in this class. In comparing the 
predicted to measured values as a function of 
frequency, it was found that the ratio was 10 db or 
less in 100 and 90 percent of the observations for 
0.013 and 10.0 Me/s, respectively. 
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The author thanks R. Penndorf, S. C. Coroniti, 
and G. E. Hill for their critical reading of the manu- 
script and helpful suggestions, and Mrs. Carolyn 
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diction curves for each station. 


Nore: Another interpretation of the discrepancies with 
the C.C.I.R. predictions will appear in a later issue of this 
journal under the authorship of W. Q. Crichlow and R. T. 
Disney of the Central Radio Propagation Laboratory, 
National Bureau of Standards. (Editor) 
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Long range navigational aids could be devised using VLF transmissions. 


gation of such waves is controlled by the lowest 
in phase velocity which accompany the diurnal 
ground monitor stations. 
consistency in fixing of about 1 nautical mile at 


1. Introduction 


The very low frequencies are those between 3 and 
30 ke/s. This paper will not be concerned with any 
specific navigational aids but will discuss those factors 
of VLF propagation which will dominate the perfor- 
mance of any aids operating in this frequency band. 
It must be stressed that this description of the propa- 

ation should not be extrapolated into adjacent 
Sends, and in particular would not apply without 
some considerable modification to the frequency 
band 70 to 120 ke/s. 

In the VLF waveband the absorption at the iono- 
sphere on reflection at large angles of incidence is 
small at all times of day and consequently propaga- 
tion is possible even in daylight at 16 ke/s under cir- 
cumstances in which there would be strong iono- 
spheric absorption a few decades higher in frequency. 
For this reason propagation losses will not be dis- 
cussed in detail and attention will be concentrated 
on the phase velocity of propagation of the signals 
because it is this velocity which determines the 
measured range from the transmitter in any phase 
comparison navigational aid. This velocity deter- 
mines the wavelength which is the unit in which all 
measurements of distance are made. In general it 
is easier to observe changes in this velocity rather 
than determine its absolute value. 

In discussion of VLF propagation we are concerned 
only with the lowest part of the ionosphere—the D- 
layer. The apparent height at nighttime is about 92 
km while that at daytime is generally in the region 
of 72 to 80 km. For waves polarized with their elec- 
tric vector normal to the ionosphere (vertically polar- 
ized), and at large angles of incidence, the signal is 
reflected without change in polarization, and the ion- 
osphere acts as a sharply bounded conducting layer. 
Most VLF transmitting aerials radiate vertically po- 
larized signals, and consequently there is little re- 
striction in practice through pcg attention in 
this paper to vertically polarized rays 

1 This is a revised version of a paper presented at the Avionics Panel Meeting 


of the Advisory Group for Aeronautical Research and Development (AGARD), 
NATO, Istanbul, Turkey, October 1960. 
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The propa- 
level of the ionosphere. The small changes 
ionospheric changes have been studied using 


Preliminary measurements made in an aircraft show internal 


ranges of 5,000 to 6,000 miles. 


2. Mechanisms of Propagation 


The VLF band is of interest for navigational aids 
because of the long distance propagation which could 
permit the development of an aid of worldwide use- 
fulness. The long-range propagation beyond the line 
of sight is achieved through repeated reflection of the 
transmitted wave between ionosphere and earth. 
This is illustrated in figure 1 which shows how, al- 
though the transmitter 7’ is not directly visible from 
the receiver Ff, signals can be received having made 
one or two hops by way of the ionosphere which at 
these frequencies acts as a reflecting layer. Signals 
could reach 2 by a path involving many reflections— 
such as that shown by a dotted line, however because 
the reflection coefficient is small for small angles of 
incidence such signals would be rapidly attenuated, 
and the number of ray paths capable of transmitting 
an appreciable signal level to 7 is limited. It is usual 
to consider the propagation between earth and iono- 
sphere in terms of the several ray paths for propaga- 
tion up to 1,200 or 1,500 km. By taking into account 
the reflection coefficients of earth and ionosphere, it 
is possible by ray theory to explain [Wait and 
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Figure 1. Jonization contour for lower ionosphere (D region). 








Murphy, 1957] the main features of this medium- 
range propagation including the maxima and 
minima in signal corresponding to the Hollings- 
worth pattern of higher frequencies. These maxi- 
ma and minima appear in figure 2 which shows 
a plot of field strength of GBR Rugby on 16 ke/s 
against great circle ‘distance from the transmitter. 
This measurement was made in an airct raft flying at 
40,000 ft from Farnborough, to Malta in daylight. 
This shows a minimum at about 230 nautical miles 
which arises through interference between the ground 
ray and first hop skywave of figure 1, and a second 
minimum at 740 nautical miles arising through inter- 
ference between first and second hop skywaves. 

At distances greater than about 1,200 km it is 
usual to describe propagation in terms of waveguide 
modes [Wait, 1957a; Budden, 1957]. The earth and 
ionosphere are then considered as the conducting 
walls of a waveguide. This waveguide differs from 
those in general use because these conducting walls 
are relatively far apart (about 214\ at 10 ke/s and 
S\ at 30 ke/s) and consequently several modes can 
be transmitted. If both earth and ionosphere were 
perfect conductors the distributions would be true 
sine waves. Thus we have the first mode which has 
a half sine wave distribution of voltage vector 
between ionosphere and earth, and the higher order 
modes which have distributions following several 
half sine waves. However, there are resistive losses 
in both earth and ionosphere and consequently the 
waves penetrate the waveguide walls. This results 
in a distortion from the true sine wave and some 
bending of the wavefront from the vertical. The 
losses in the ground and ionosphere are most marked 
for the higher order modes and consequently these 
modes are reduced to low amplitudes before they 
have traveled far. As a result of this the first mode 
usually dominates at distances of 2,000 km or more 
although for some ionospheric conditions the zero 
order mode may be the most important at the 3 ke/s 
end of the frequency band. 


3. Attenuation of Signals 


This is a facet of the subject of great complexity. 
Towards the lower end of the VLF band there is 
marked absorption while above 10 ke/s the attenua- 
tion is of the order of 2 db/1,000 km. In general the 
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Ficure 2. JLonospheric reflection of VLF waves. 








losses for transmission paths over sea are less than 
those for a corresponding overland path by a factor 
of as much as 2:1. A detailed theoretical study 
[Wait, 1957b] has been presented which expresses 
the frequency variation of attenuation for each mode 
in terms of the height of the ionosphere, and the 
conductivities of ground and ionosphere. Inasmuch 
as it is difficult to determine the properties of the 
lowest layers of the ionosphere except by VLF 
studies, it is not easy to confirm these theoretical 
predictions in detail. However, there is no doubt 
that they are a useful guide in the interpretation of 
the experimental phenomena. Figure 3 shows the 
measured variation of signal with distance for an 
over sea path from San “Diego— Hawaii compared 
with that which would be predicted theoretically 
[Wait, 1957a] on the basis of chosen parameters for 
the ionosphere. The similarity between the forms 
of the two curves is pleasing, and the general agree- 
ment is as good as could be expected be: aring in mind 
the somewhat arbitr ary assumptions made about the 
ionosphere. 


4. Velocity 


In any navigational system at present contem- 
plated the observations will take the form of phase 
comparisons made either between a number of re- 
ceived signals, or between received signals and 
high stability oscillator carried by the aircraft. For 
this reason we are interested in the phase velocity of 
the total wave. At distances from the transmitter 
for which the waveguide mode theory is applicable 
it is to be expected that the phase velocity would be 
somewhat greater than unity because of the finite 
spacing of eround and ionosphere, on the other hand, 
the effect of losses in the waveguide walls would be to 
slow down the wave. Some measured values have 
been published. 

In the region near to the transmitter the effective 
phase velocity at any range can be deduced from the 
variation with radial distance of the resultant phase 
of the ground and sky rays, allowance being made for 
their differing ray paths and the changes of phase 
and amplitude on reflection. 
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FiGurE 3. Variation of signal amplitude with distance. 
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5. Diurnal Changes in Phase Over a Given 
Path 


If the phase of a signal received at a distant point 
is compared with that of a standard oscillator, then 
changes in velocity of propagation will become ap- 
parent as variations in relative phase. It has been 
found that there are diurnal variations in the phase 
of the received signal associated with changes from 
daylight to darkness. These are well illustrated in 
figure 4 in which J. A. Pierce [1957] depicts his 
“phase trapezium.” The far left of the figure refers 
to times at which the whole 5,200-km path between 
the transmitter at Rugby and the receiver at Cam- 
bridge, Mass., is in darkness. At a short time before 
sunrise at Rugby a steady phase change commences, 
ending as the sun rises at Cambridge. During the 
day the phase remains constant until sunset at 
Rugby when the reverse process’ begins until the 
steady night value is reached at Cambridge sunset. 
The phase change is in a sense corresponding to a 
daytime increase in phase velocity. It is explained 
as the result of a reduction of the apparent height of 
the ionosphere during the day. This would be ac- 
counted for on the ray theory by saying that with the 
lower ionosphere the sky rays had a shorter path, or 
on the basis of waveguide modes, that the phase 
velocity had been increased by bringing the wave- 
guide walls closer together. 

It may be noted that the sharp onset of these phase 
changes occurs not at ground sunrise but a little 
earlier, at the time at which the sun viewed from a 
height of about 90 km above the receiver would be 
seen to rise above the ozonosphere. 

The phase change of figure 4 has been found to be 
predictable from day to day within a few micro- 
seconds—an uncertainty in prediction which would 
introduce a positional error of about 1 mile in a 
navigational aid on a transatlantic flight. A series 
of measurements of the diurnal phase changes of the 
signals received from GBR have been made at Malta. 
These have shown variations which were basically 
repeatable from day to day but were more complex 
than those reported by Pierce. The changes for one 
day is shown in figure 5 from which it will be noted 
that instead of the neat trapezium there is a slight 
reverse drift before dawn followed by an overswing of 
phase before the steady day time value is reached. 
It has been suggested that these phase swings arise 
through the beating of signals reflected from two 
layers or that they represent the variations in phase 
of the sum of several modes, each mode being diff- 
erently affected by the changing ionosphere. It 
would be expected that such overswing would be 
most apparent on a path which lay roughly parallel 
to the dawn or dusk line, because for this dawn or 
dusk would occur almost simultaneously along the 
whole path, all parts then being in the same inter- 
mediate state between day and night. For other 
paths where at any time parts are in different inter- 
mediate states, the overswings are averaged out. 
Even if the path follows the dawn of dusk line the 
effect of overswings would not be expected to appear 
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Fiaure 5. The diurnal variation of the time of arrival of the 
GBR signal at Cambridge, Mass. (after Pierce, Proc. IRE, 
June 1957). 


at distances so great that only a single mode is re- 
ceived at an appreciable amplitude. At present it 
is not possible to check these theories because suffi- 
cient data are not available. 

The magnitude of the diurnal phase change varies 
with distance from the transmitter. At distant 
receivers the ground ray will be negligible, and in 
order to compare signals near to the transmitter with 
those received at a distance it will be most convenient 
to neglect the ground ray and consider only the 
resultant sky ray. 
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At a point near to the transmitter the sky ray will 
be reflected at near normal incidence, and so when 
the reflecting layer changes by 18 km in apparent 
height from day to night the sky-ray path will change 
by 36 km—or a time of say 120 usec. As the point 
of reception is moved further away from the trans- 
mitter the change in sky-ray path length will decrease 
rapidly at first and then more slowly as glancing in- 
cidence is approached. For greater distances the 
waveguide mode theory would predict a linear in- 
crease in phase change with distance. The initial 
rapid fall and linear rise would combine to give a 
variation of phase change with distance of the form 
of figure 6. Various measured values of phase change 
have been plotted and show a reasonable agreement 
with the expected shape of the curve. If this curve 
can be established by further measurement at other 
ranges up to say 10,000 km, then the diurnal varia- 
tions of phase for a VLF navagational aid can be 
predicted with accuracy. 

Such programs of monitoring diurnal changes 
have been planned during the summer of 1961. It 
is hoped to monitor simultaneously over 10 to 20 
paths of a variety of distances. Some of these paths 
will lie entirely in the middle latitudes while others 
will be entirely tropical and others will cross the 
Arctic icecap. 


6. Effects of Ionospheric Disturbances 


VLF propagation can be affected by a magnetic 
storm or by a “sudden ionospheric disturbance,”’ 
(SID). The effect of the magnetic storm is to give 
an increase in the fluctuations of phase during the 
night but no serious decrease in amplitude. Fluc- 
tuations of up to 5 usec at 16 ke/s lasting 2 to 3 min 
have been reported [Pierce, 1957]. Following the 
main part of the storm the ionosphere remains 
disturbed and the diurnal effects are Jess than normal, 
recovering over a period (up to a few weeks in some 
eases). A recording of phase against time during a 
magnetic storm is shown in figure 7. This was over 
the path Rugby—Malta (2,210 km) for a frequency 
of 16 ke/s. It will be noted that the steady day phase 
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Ficure 6. Diurnal variation in phase measured at Malta. 





is little different from that at night, but that the 
phase swings have the amplitude expected for a 
normal day. The features of this record are not 
yet explained although it is expected that the pres- 
ence of the overswings at a time when the diurnal 
changes were small will be significant in the explana- 
tion of the overswings on normal days. It should 
be noted that although the effect of the magnetic 
storm in this case was to decrease the diurnal effect, 
there are circumstances in which the effect may be 
the reverse. 

During the solar event of November 1960 the 
diurnal changes at Malta showed even less change 
from steady night to steady day, but the overswings 
were of full amplitude. A simultaneous phase 
record of GBR taken at Nairobi showed a trapezium 
of full amplitude. This diurnal change recorded at 
Nairobi is shown plotted as point 8 of figure 6. It 
will be noticed that it lies well on the curve, and is 
unaffected by the magnetic disturbance. 

In order to assess the effects of magnetic storms 
on the accuracy of a VLF navigational aid, it will be 
necessary to obtain a wider range of test data by 
simultaneous and continuous recording on a world- 
wide basis. This study must determine how wide- 
spread and geographically uniform are the effects of 
a storm, and in particular, how the curve of figure 6 
is altered. It may be that the effects of magnetic 
storms will best be countered by the use of an extra 
received signal which would help to establish and 
allow for changes in scale. 

Compared with a magnetic storm an SID is of 
much less consequence. It is associated with a 
solar flare and effects only the sunlit hemisphere. 
As far as VLF propagation is concerned it alters 
the phase change on reflection at the ionosphere so 
as to give an advance in phase. The effect reaches 
its Maximum in a few minutes and then decreases 
more slowly to regain a normal level within an 
hour. In general the amplitude of these changes 
for a given path is much less than the corresponding 
diurnal variations, but further data are required, 
and these could be expected as a byproduct of any 
monitoring such as that suggested for the study of 
magnetic effects. 


















































7 |_| | | | |_| 
| | | 
1S) | | 
® 100 | = = 
=i DIURNAL CHANGE IN TIME OF ARRIVAL OF TOTAL 
SKY WAVE v DISTANCE FROM TRANSMITTER. 
- 80 4 4 os 4 j 
wi | 
© = Fond 
Zz } 
aI 60 st 
3 \ 2 
| 5 2} ne — 
| 3 as all | 
| i 
° 1 2 3 a 5 6 7 8 





DISTANCE FROM TRANSMITTER, I0> km 


Diurnal change in time of propagation versus 


| 
| Figure 7. 
distance. 


578 


na 


re 
be 


an 
as: 
us 
sti 
an 


no 





7. Airborne Measurements 


In order to assess the potentialities of a VLF 
navigational aid a series of flights have been made by 
the Royal Aircraft Establishment, Farnborough. 
For this work a high stability oscillator was fitted 
in a Comet aircraft together with receivers which 
could compare the phase of this oscillator with the 
signals received from the transmitters GBR (Rugby, 
16.0 ke/s) and GBZ (Criggion, 15.2 ke/s). 

Suppose that the aircraft is directly above one 
of the transmitters and the phase of oscillator and 
transmitter compared, then if the aircraft moves 
away the phase of the transmitter will recede because 
of the time taken for the signal to pass to the air- 
craft. The recession of transmitter phase is a 
measure of great circle distance from the transmitter. 
If the aircraft is on a straight course at a steady 
speed, a plot of phase change against time will be 
a smooth curve. In figure 8 is shown a plot of such 
phase changes measured for GBR and GBZ. The 
slight scatter of the measured points about the 
smooth curve is equivalent to an uncertainty in 
great circle range of 1 nautical mile. Using these 
two transmitters this consistency has been obtained 
at distances up to 1,200 nautical miles from the 
transmitter. 

Further airborne experiments have been carried 
out comparing the phases of the signals from GBR 
and NBA (Balboa, 18.0 ke/s). This has shown that 
it is possible to track the changes in position of the 
aircraft using changes in ranges from the two stations. 
Much of the flying has been over the Mediterranean, 
and this has shown that the fixes obtained at ranges 
between 5,000 and 6,000 miles from Balboa on any 
given flight are internally consistent to within a 
nautical mile. Their absolute accuracy is not easy 
to determine because of the difficulty of locating 
the position of a high flying aircraft. However, 
a further series of flights is planned during which 
time theodolite and other photographic methods of 
determining position will be used. 


8. Conclusion 


The general picture of icnospheric effects on VLF 
propagation is that there are no fundamental 
reasons why worldwide navigational aids could not 
be devised in this frequency band. However, there 
is a lack of data to fill in the detail of the picture, 
and only on the basis of such details could one 
assess the ease and convenience of navigation 
using such a system. In this field, theory has out- 
stripped experiment, and the need is for an extensive 
and worldwide monitoring of the phase stability of 
VLF transmissions over long paths during both 
normal and disturbed ionospheric conditions. 
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Figure 8. Diurnal variation of relative phase during a mag- 
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On the Spectrum of Terrestrial Radio Noise at 
Extremely Low Frequencies” 


H. R. Raemer 


Contribution from Applied Research Laboratory, Sylvania Electronic Systems, Waltham, Mass. 
(Received May 30, 1961; revised June 15, 1961) 


A theory of the frequency spectrum of radio noise at extremely low frequency (ELF) is 
presented and the results compared with recent measurements of the first five ““Schumann”’ 
resonant modes (between 8 and 34 c/s) made by Balser and Wagner [1960]. The source of 
this noise is assumed to be return strokes in vertical cloud-ground lightning flashes distributed 
randomly in time, uniformly in angular displacement along the earth relative to the observer 
and with statistics of stroke duration, interstroke intervals and strokes per flash taken from 
studies of thunderstorms reported by J. C. Williams. Thus, the mathematical model for the 
noise sources is an extremely simple one, being analogous to the shot effect in electron devices. 
The electromagnetic model employs the familiar waveguide mode theory, assumes a sharply 
bounded homogeneous ionosphere, and neglects the earth’s magnetic field. 

Agreement between the shape of the theoretical and observed spectrum is good for the 
first three modes and rather poor for the higher modes. 

It is found by matching the theoretical resonant frequencies to the observed resonances 
that the product of effective ionosphere height h and the square root of effective conductivity 
vo; is a decreasing function of frequency. The functional dependence of this quantity on 
frequency is determined and used in the calculation of the mode spectrum. 

Discrepancies between the theory and experimental results are believed to be partially 
due to the artificiality of the sharply bounded homogeneous ionosphere model and to failure 
to give sufficient probability weighting to equatorial regions of abnormally high thunder- 
storm activity. 

These last items are the subjects of continuing work on the extension of the theory. 


1. Introduction 


The theory of extremely low frequency (ELF) propagation, as originally developed by 
Schumann [1957] and extended by Wait [1960a] and others, predicts the existence of a series of 
modes whose resonant frequencies, relative amplitudes, and Q’s can be found by considering the 
earth-air-ionosphere cavity as a linear filter and calculating its frequency response function 
directly from electromagnetic theory. Recently, Balser and Wagner [1960] have observed these 
modes and experimentally determined the detailed shape of the mode spectrum. 

The basis of the cavity analog is the classical Sommerfeld theory of propagation around a 
spherical earth, which leads to the representation of the field between earth and ionosphere as 
a series of “waveguide modes.” If the ionosphere height is extremely small compared 
to wavelength, as is true at ELF, the electric field is essentially radial and a single waveguide 
mode is sufficient to describe the fields at a distance from a localized source. This waveguide 
mode when considered in the time domain can, in turn, be shown (by an appropriate expansion 
of the hypergeometric series into a series of Legendre polynomials) to be a superposition of 


. ; ; : ‘ sn Lc ——— 
waveforms each closely approximating a damped sinusoid of frequency ag gyn ), where c 
oT 


and a are velocity of light and earth radius, respectively and n is a mode index number running 
from zero to infinity. This phenomenon is entirely analogous to the ringing of a resonant 
cavity in response to impulse excitation. 

Experimentally observed VLF noise is thought to consist of a superposition of responses 
from lightning flashes all over the world |Williams, 1959]. Therefore, assuming that this 
hypothesis is correct down to ELF, a theoretical model based on it should yield accurate 


1 Work supported by the Office of Naval Research. 
2 A preliminary note on this work appeared in the J. Geophys. Research, May 1961. 
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calculations of the spectrum of observed ELF noise. It is the purpose of the study reported 
here to construct such a model and use it to calculate the expected ELF noise spectrum. The 
results are compared with the measurements of Balser and Wagner [1960]. 


2. ELF Cavity Mode Theory 


Following Wait [1960a] one can express the significant (radial) component of the 
electric field at an observation point due to the harmonic excitation (at frequency w) of a 
vertical electric dipole of moment Jds on the surface of the earth in the form 


E, (iw) =Ids 33 C,,(0) F,(iw) (1) 


n=0 


where the ionosphere is assumed to be electrically uniform and sharply bounded at height 4h, 
and the following definitions apply: 


C,(0)= eas P,(cos 6)(2n+1) 


éo>= Dielectric constant of free space 
P,,(cos 6)=Legendre polynomial of order n 


6=Polar angle relative to the source 


and 
aay 1 —w?(w) 
P(t) tw [—w?(tw) +? ] 
with 
: C . q 
w?(iw) vt v(v+1)=— (tw)?—a(iw)3/? 
- oe 
R=; n(n+1) 
a 
and 
a= ~< = (1 +z 
helps Mo; Oe 
where 
o,=0, +l, ~o,=complex conductivity of the earth 
and 


ane { Gta) =¢,+ Wwe; ~o,=complex conductivity 


of the ionosphere, where N, e, m, g, are the electron density, charge, mass, and collision frequency 
(Stratton, 1941], and it is assumed that w << g. In the calculation of the appropriate iono- 
spheric electrical parameters the effects of the earth’s magnetic field and the motions of the 
heavy ions are neglected; in addition, the electron collision frequency is taken to be independent 
of the nonrandom electron energy. In general |¢,|>>|¢,|, so in what follows a, which will be 
called the “ionospheric loss parameter,” is approximately a real quantity, equal to 1/(Ayyo,). 

The function E,(iw) is actually the Fourier transform of the field in response to a (vertical) 
current impulse of magnitude J at 6=0 and at t=0, i.e., it can be regarded as the characteristic 
frequency response function of the earth-air-ionosphere cavity. F,,(iw) is the Fourier transform 
of the nth mode excited by an impulse of dipole current of magnitude J. From knowledge 
of all of the F,(tw)’s one can predict the frequency response to any time varying source current 
using standard linear system theory. 
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3. Superposition of Fields From ELF Sources 


The radial field component is expressed here in a spherical coordinate system wherein the 
source is at the pole. Because of azimuthal symmetry about the source, we can without 
changing the form of (1) consider L,(w) as expressed in a system with pole at the observer and 
source at polar angle 6. Suppose that lightning flashes, assumed equivalent to vertical dipole 
sources, are distributed around the earth’s surface. We define £,(w;0,6) as the Fourier trans- 
form of the radial field component at the pole due to a vertical dipole source at the point 6,¢. 
The source current moment waveform at 6,¢ and its Fourier transform are denoted by 7(t;@,¢) 
and J (w;6,6), respectively. 

The voltage induced in a receiver at the observation point at time ¢ due to the presence 
of lightning sources all over the world is proportional to 


t 


N ; 
o(t)=>) dt’ f(t—t’; Oxby) V(t’; Ox,be) (2) 
k=1 


where the index k is used to number the source positions, from 1 through N, and 


o 


f(t; 0,6)=>5 (2n+1)P,(cos 6)f,(t) 


n=0 


where f,(t) is the inverse Fourier transform of F,,(7w), i.e., the impulse response of the nth 
earth-air-ionosphere cavity node. 


4. Calculation of the Power Spectrum 


The evaluation of the power spectrum of v(¢) falls within a familiar class of problems in 
random noise theory, specifically, the spectral analysis of the voltage arising from the “shot 
effect”’ in electron devices. This problem has been treated by Rice [Wax, 1954] and others 
[Davenport and Root, 1958]. The analysis to follow is mathematically equivalent to that of 
the shot effect. 

Consider a voltage, v(t), appearing at the output of a linear filter with impulse response 
F(t) due to a train of input pulses, each of which is denoted by é(t—t,;8). The shape of the 
kth pulse is determined by the components of a vector 8, its magnitude (time integral of ampli- 
tude over duration of pulse) is az, and its time of occurence is & The output voltage is thus 


eo 


r(t)= Sa) f(t—t/)k(t’—te; dt’ (3) 
. -——_ 

We assume a time interval * —7<t<T over which all the action takes place and that there 
are exactly N pulses within that interval. Because of finite filter “memory” and the filter 
realizability requirement,* we can regard as zero the response to pulses occurring at values ¢, 
outside the interval (t— 7) to t. Here J) is the time of occurence of the earliest imput influ- 
encing the filter output at time ¢. The interval (f—7}) is less than T but somewhat greater 
than the filter time constant. Numbering from 1 to N the pulses within the interval — 7 to 7, 
and using (3) we can write 





“~ 
v(H=D) ag (t—t; By) (4) 
k=l = 
where 
t 
g(t—t,; B= f(t—t’ )é(t’ —t,; B,)dt’; —T<t,<T; bent. .i, N, 


t-To 


the response of the filter at time ¢ to a pulse occurring at time ¢;,. 


’ The time T' washes out in the final results; the finite time assumption is strictly for mathematical convenience in the derivation. 
‘i.e., the requirement that the filter can respond only to past (and not to future) inputs. 
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We now assume that: (1) The a,’s are independent and identically distributed random 
variables with common probability density function p(a); (2) for a fixed number of pulses 
within the interval t— 7% to ¢t, the ¢,’s are statistically independent random variables, each 
uniformly distributed throughout the interval — 7 to T; (3) the probability of finding exactly 
N pulses within the interval t— 7) to t, denoted by p(N), is independent of the statistics of the 
a,’s and ;’s; (4) the 6s are independent and identically distributed random variables with 
common probability density function p(8); and (5) the random variables a,;, t,, and 8, are 
all mutually statistically independent. 

By calculating the autocorrelation function of v(t) under the assumptions (1) through (5) 
and taking its Fourier transform, we obtain for the spectrum ° 


A 


S(o)=y73 {(N( N—1))(a)*(g(t) g(t +7))6(w) } + S(w) (5) 
where 
A N)<a*) (7 ' 
Sw) = dBp (B)| G(w; B)? 
= All B-space - - - 
and 
G(w; B)= dte~'*g(t; B) 
By the well-known convolution theorem for Fourier transforms 
G(w; 8) =F (iw)E(@; B) (6) 


where F(iw) is the frequency response function of the filter and £(w;8) is the Fourier transform 
of the pulse. wi 
The first term of (5) is a “DC spike,” and is not required to determine the shape of the 
mode spectrum in the frequency region of interest. We will, therefore, concern ourselves only 
with §(w). 
5. Statistics of Lightning Source Currents 


The result (5) can be adapted to the calculation of the ELF noise spectrum by allowing 
the pulses ¢(t—t,;8,) to represent lightning flashes whose shape parameters, time of occurrence f;, 
and polar angles 6, are assumed to be random variables, statistically independent when 7 #4. 
The angles @, are included among the components of the vector 8;, as are the number of strokes 
in a flash [Williams, chapter VI], and the shaping parameters associated with individual strokes. 

The statistics of the lightning flashes are assumed to obey all of the assumptions (1) 
through (5). It does not necessarily follow that we can with physical justification define a 
spatial probability density function p(@) apart from consideration of other random variables 
associated with flashes. However, this assumption is basic to what follows and we note that 
it is tantamount to the assumption that statistics of lightning flash shape (i.e., those types of 
flash important as sources of VLF noise) are independent of those of terrestrial location of the 
flashes. The essential correctness of the latter assertion is supported by experimental evidence 
[Williams, tables 6-1, 6-2, 6-4, figures 6-1, 6-2, 6-7, 6-11, 6-14]. 

Using (5) and (6), the definition of f(t) given in (2), and the assumptions about lightning 
flashes discussed above, we obtain for the ELF noise spectrum (exclusive of the DC component) 


§(w) =K(E(w,8’)|) 33 Sno) (7) 


m,n=0 
where K is a constant, 8’ represents all lightning flash parameters exclusive of @, and 
Snn(w) =(2m+1)(2n+ 1) F,, (iw) F* (tw) Linn, 


5 Ensemble averaging over Bx is denoted by <_.>, averaging over the time variables tz by a line over the symbol. 
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where 


Dna =20 | dé sin 6P,,,(cos 6)P,,(cos 6) p(@), 


(JE(w; B’) = | dB’ r(a")|Ews B’)|? 
All 8’ space 


p(8’) is the probability density of the set of parameters 8’, and p(@) sin 6d@ is the probability 
that the lightning flash is located between @ and 6+d8. 

We are interested only in relative and not absolute amplitudes; therefore, in numerical 
computations we will fix K by assuming a normalization that equates the average power levels 
in the theoretical and experimental results. 


6. Flash Shape Parameters 


J. C. Williams [ch. V, sec. 5f] has concluded from his studies of thunderstorms that the 
overwhelming cause of VLF noise is the presence of vertical currents in the return strokes 
associated with cloud-to-ground lightning flashes. These flashes consist of from 2 to 12 strokes 
of similar shape. He bas compiled statistics on flash parameters, including the shape of the 
median return stroke [ch. VI, sec. 1, 4, 5, and ch. VII, sec. 1], interstroke interval, and number 
of strokes per flash. He has also made spectral calculations for the VLF noise resulting from 
the flashes, but these are in a higher frequency range and inapplicable to our problem (ch. VII, 
sec. 3). 

The current moment of Williams’ median return stroke is proportional to 


7 4 ‘ ) / (38 )-6) 
(—t)=(| nt’ )at’ ) i(t—t,) = { —16.8(108(e~ F—fe! O-7X18) 4 15,35 (10-3) e~ Fe) XE 
“ 0 c 


Mien, oe 5 -4) | re 3\ ,— (t—t,) / (6.8 —3 f —5.55 3) (¢—t,) >’ 
+103¢ (t—t,.) / (5X10 +.0.45(10°)e (t—t,.) / (6.8X10 }u(t—t,) {1—e 5.955 (10 t—t,) } (8) 


where 2(t) is the velocity of the lightning charge along the cloud-ground channel, the unit step 

is inserted to assure that the stroke has zero value before time ¢, and the lightning current 

i(t—t,), v(t) and t, are measured in amperes, kilometers/second and seconds, respectively. 
The normalized spectrum of (8) is proportional to 


T, ‘> 


nes 4 1 | 
T(w) |2=| e'** I ——- ——____!__—. } | 9 
Fic | >> ie l—%or, 77-7, ( | (9) 


where J, and 7; are respectively the currents and time constants given in (8) and r,=1.80(107*) 
sec. 

It would be desirable in the interest of accuracy to consider the shaping parameters of the 
stroke model as random variables and average the spectrum over these variables. This is a 
formidable mathematical problem in general; and moreover, there is insufficient experimental 
data available to assign probability density functions to these parameters, a measure which 
would be required to evaluate the averages in the general case. The assumption that all shaping 
parameters have small standard deviations relative to their mean values would reduce the 
problem to one requiring only a knowledge of first and second moments of the distribution of 
each parameter. However, experimental knowledge of even these quantities is not available 
in enough detail to allow any meaningful assignment of numbers to the calculation of the 
spectrum of the median return stroke itself. This is tantamount to assuming that all return 
strokes are nearly identical in shape. This is not quite accurate, but the analysis will still 
provide useful insights into the effects of stroke shape on the ELF spectrum. 
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The waveform of a flash consisting of . identical strokes ° is 
£(t)=>3 i(t—t,), (10) 


where 7(t—t,) is given by (8). 
The spectrum of a lightning flash is roughly proportional to the averaged absolute square 


of the Fourier transform of &(t), i.e., 
- es 
EGY =1f@) IF 2a cele) (11) 
[=1 


where the indicated average is over the ensemble of the interstroke time intervals ¢;—*,. 
An obvious classification of terms in the summation of (11) results in 


(E(w) |2) Pe ee 
scuciek'’* Sella hceial aati (12) 


lxk 


where JT is the average interstroke interval and AT, is the deviation of the interval (¢,—t,) 
from (l—k)T. 

To calculate (e'*47%) accurately, we would require a knowledge of the probability density 
function for AT,. An extremely rough estimate of this function is provided by Williams 
ch. VI, figs. 6-7], but it is not specified in sufficient detail to allow a precise representation. A 
rough processing of Williams’ curve shows a Gaussian distribution to be a reasonable approxi- 
mation for this variable, with standard deviation o,, the latter being assumed somewhat smaller 
than 7. The result of a simple integration is 


/ 1wAT jj. 


(e le (13) 


\—=¢ 
In the absence of experimental data on second order statistics of A7;,, it will be assumed 
that the random variable (¢;—(¢;)) is not strongly correlated with (¢,—(¢,)), where k#l. 
It follows that 
07,07 for all 14k (14) 
Invoking (13) and (14) in (12), using the power series for 1/(1—.r) to bring the summation into 
closed form and averaging over all possible values of VM, we arrive at the results. 
p(w)=>) py(w) pW) (15) 
M 


- 
=0 


where p(.M) is the probability of ./ strokes in a flash and 
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par(w)= M(1—e- 7) + ere? (16) 
ee 
—_ = 

Averaging (16) over values of .W/, we obtain 

p(w) (M)(1—e-*?) +e-"? SA, cosn oT (17) 
n=0 


where 
x 
A.= >3 (M—K)p(M). 
M=K 
Actually the strokes are not quite identical, the first being about twice as intense as succeeding ones. However, the approximation will not 


introduce serious errors into an analysis of this type. 
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The parameters 7 and oc, have been roughly estimated from Williams (fig. 6-7) and found 
to be about 33 and 7 msec, respectively. The precision of these estimates is questionable, 
but they are probably correct to within 20 percent. The same remarks apply to similar rough 
estimates made of p(M) from Williams’ cumulative probability distribution curves for the 
number of strokes per flash (fig. 6-1). With these estimates we conclude that the probabilities 
p(1) through p(10) can be roughly approximated by the numbers 0.15, 0.22, 0.23, 0.15, 0.09, 
0.06, 0.04, 0.03, 0.02, 0.01, and 0, that (M) is about 2.5 msec, that the numbers Ay through 
Ag are roughly equal to 3.5, 2.5, 1.6, 1.0, 0.61, 0.36, 0.20, 0.10, 0.04, and 0.01, and that A, 
effectively vanishes for n>9. 


7. Geographical Distribution of Lightning Activity 


The assumption is now introduced that the probability of occurrence of a lightning flash 
within a range of angles 6, 6+d6 during the time interval between t— J) and t is 


or 
sin ap(odo= | do(p(8,¢) sin 6d6)=(po+p(6))dé sin 6; 0<0<r (18) 
0 


where p(6,)dé@ sin 6d¢ is the flash occurrence probability for a region of solid angle dé sin 6d¢ 
at 0,6, po is the part of the 6 distribution that is uniform in the angle 6, and D(@) accounts for 
the existence of regions of abnormally high or low average thunderstorm activity. 

Initially, we will neglect p(@) and assume a completely uniform distribution po. The 
result of applying this assumption to (7) (see appendix) with the aid of the orthonormality 
of the Legendre polynomials is 


S(w) a ae) 
=} (2n+1)|F,, (tw) |? a9) 
Po |E(w, B’)|?) =o 


Considering now the term p(@), it has been established that there are three equatorial 
regions of abnormally high average thunderstorm frequency [Handbook of Geophysics, 1960]. 
These regions are located in (1) South America, (2) Africa, and (3) the Southwest Pacific area. 
We will denote them respectively by R,, R:, and R;. The regions center approximately at 
(50° W, —10°S), (15° E, + 10° N), and (110° E, +10° N) respectively, and during all seasons 
of the year they show an average thunderstorm frequency somewhat higher than is observed 
elsewhere in the world. Speaking more quantitatively, the spatial probability density in 
those regions is from 100 to 1,000 times higher than that throughout most of the earth’s surface. 

There is another near-equatorial region in Southeast Asia (~90° E, +20° N), showing 
substantial year-round thunderstorm activity, but not as high as in R,, R2, and R;. This 
region has not been accounted for in our calculations. Certain nonequatorial land regions 
show intense activity during their summers, e.g., the Southeastern United States (centered 
about (75° W, +30° N)), but this has also been excluded from our analysis. It is desired 
here to concentrate on those features of the ELF noise spectrum that are not season-dependent, 
and this objective is best fulfilled by studying the contributions from the three most prominent 
equatorial regions R,, Rs, and R3. 

In view of the approximate nature of meteorological data, it would not be feasible to attempt 
to construct a precise thunderstorm probability function for the regions under investigation. 
We will, therefore, construct an idealized theoretical model which is convenient for computation 
and which embodies the essential features of the observed distributions. 

We will approximate the probability density functions within the regions ?; by Gaussian 
functions 


dm = 2/92 
, — (0—8;)" 20; 


= (20) 





y 210} 








pO=K; 


where K; is a constant, 6; is the center of the distribution in R; and go; is the standard deviation. 
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Using (7), (18), and (20), we have 


4rPo 5 , N a 
Lon Dan +27 p Kjlinns (21) 
j=1 


2n+1 


where 
eo 0;)? /20? 


bear | d6 sin 6P,,(cos 6)P,,(cos Sh Sameer 
V2 20? 
The integrals Z,,,; have no simple orthogonality property to aid in their evaluation. The 
method to be used here will involve the expansion of the Legendre polynomials P,,(cos @) in 
a Fourier Series of the form [Jahnke-Emde, 1945] 


P,,(cos 0=)=>) dy cos pb (22) 


u=0 
where the coefficients @,, are given in reference 9. From (21) and (22) we have 
m n 
lnaj= =>) > ¥ Anylnv mpunv (23) 
u=0 u=0 
where 
r a= (6—6;)*/20; 
isis ={ dé sin 6 —=— c0s 46 cos vO 
0 \ 2105 
Using trigonometric identities, assuming (justifiably, with the parameter values of regions 
R,, R», R;) that (1) o; is small compared to z and (2) 6; is not near 0 or x (so that the limits of 
the integral are effectively —@ and +), and integrating by completion of the square in 


the exponential, we obtain 
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for k>0 where [5-1 = “| means 5—1 if & is oe” 5 if k is odd, where the subscripts R 


- ~ 


and J refer to “Real part” and “Imaginary part” respectively, and where 


[a,(w) |r=>> (2n+1)[F,, (tw) |2Ony 


n=0 


[a,.(w) >= 2S (2n+1)[F, (to) |r@nu 
n=0 
Normalization of p(@) can be accomplished by integrating (18) subject to the above 
assumptions (1) and (2). The result is 


N —¢?/2 
2po+> 5 K;sindje © =1. (25) 
j=1 


The condition (25) will serve as a constraint on the relative values of p) and the Aj co- 
efficients. 
8. Results and Conclusions 


Note that the theoretical resonant frequencies for a lossless ionosphere are given by the 
De Pees ie : , <<, 
numbers (2m) a n(n+1), which for c/a=46.7 are 10.5, 18.2, 25.6, 33.1, and 40.5. These 
a7 ) 


are sufficiently different from the experimentally observed resonances to warrant the conclu- 
sion that ionospheric losses play a major role in the mode spectrum. The real positive 
quantity a, equal to 1/(hyuo;) is a measure of these losses. To study its effect, the functions 
w,F',(iw) were calculated for various values of a. In the earliest work, the mode series was 
computed for a few values of a (assumed independent of frequency) corresponding to h~80 
km (a compromise between day and night ionospheres) and values of o, between (1075) and 
(10-*)mhos/m. It was found that deviations from observed mode resonant frequencies were 
very substantial (e.g., 20 to 30 percent) and that no constant value of @ would produce near- 
agreement on all five observed modes. The next approach was to allow a to vary with fre- 
quency and use the Balser and Wagner results to infer its frequency variation. This was 
initially done by graphical methods based on plots of real and imaginary parts of w,/’,(iw) 
[Raemer, 1961]. The method has been revised to increase accuracy. The function | w,F, 
a 


(jw)? was computed in terms of normalized variables 7,=/w, and y, 5 (plots are 

Vn 
shown in fig. 1) and its derivative with respect to x, equated to zero to find the maximum. 
Corresponding values of y, were then found for n=1 through 5, resulting in a function a(f) 
required to match theoretical and observed resonant frequencies. This function was then 
fitted to an analytical curve, by the method of least squares and (see fig. 2) a satisfactory fit 


was found to be 
a( f)=7.649—0.1863f +0.01362/2—0.00006602/?— 0.000003388/* (26) 
Note that this differs from the function originally reported by the author [1961], i.e., 


a(f)=6.44 logio f—0.0131(f—S8.0) (27) 
This is due to the difference in computational method, and the function (26) should be regarded 
as the more accurate of the two. 

The fundamental assumption behind the mode-matching technique described above is 
that the observed modes are those arising from the individual nth order frequency response 
functions of the earth-air-ionosphere cavity and are not significantly influenced either by cross- 
terms between modes, overlap between modes, or the details of the lightning flash spectrum. 
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For the moment we confine attention to the uniform angular distribution of sources. 
Using (26) to represent a, the quantity given by (19) has been calculated and is shown plotted 
in figure 3. 

The plot of figure 3 is not appropriate for comparison with the experimental curve, since 
it assumes a white lightning flash spectrum and neglects the fact (see (1)) that the field ampli- 
tude contains a factor proportional to the reciprocal of the effective ionosphere height A. If 
it is assumed that the frequency dependence of a@ is due entirely to the effective ionosphere 
height and not to frequency dependence in the effective conductivity o;, then for consistency 
the power should contain a factor [a(f)]’. Multiplication of the quantity plotted in figure 3 
by this factor leads to degradation rather than enhancement of agreement with experimental 
results, i.e., the envelope increases rapidly with frequency. Assuming our basic model to be 
essentially correct, one of two alternative possibilities follows: (1) The frequency dependence 
in @ is entirely due to h and the spectrum of lightning flashes within the frequency range of 
interest is substantially different from a “white noise” spectrum or (2) the frequency dependence 


in @ is largely due to o;, in which case the lightning spectrum still cannot be “white” noise: 


but need not decay very rapidly with frequency within the range of interest. 

Assuming that the ELF noise is largely caused by return lightning strokes, the suspicion 
that there is a significant decay in the spectrum with increasing frequency is borne out by figure 
4. This figure, showing the estimated spectrum of a lightning flash compounded from (9) and 
(17) shows an average spectral decay of between 3 and 4 db per octave of frequency in the 
region of interest. Inserting this factor into (19) results in the theoretical curve (a) of figure 5. 
Incorporating the frequency dependence of fA into the calculation results in the theoretical 
curve of figure 6. On both figures 5 and 6, the experimental curve is also shown for purposes of 
comparison. In both figures the experimental curve has been normalized in such a way that 
its average power is equal to the average power of the theoretical curve. 
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assumed frequency dependent. 


The minima in the experimental curves are somewhat deeper than those of the theo- 
retical spectra and the envelope of the minima, which is roughly constant on the experi- 
mental curves, decays substantially with frequency in the theoretical results. The rate of 
decay of the envelope of the maxima is’ comparable, but too rapid at the high frequency 
end of the region of interest. Widths of the modes are roughly comparable. Generally the 
agreement of spectral shape is seen to be reasonably good near the lower order modes and 
rather poor near the higher order modes. On the theoretical curves of both figures 5 and 6, 
the fourth mode becomes barely visible and the fifth is completely suppressed. Agreement is 
somewhat better in general in the case where the frequency dependence of 1/h? is not accounted 
for, providing some evidence that the frequency dependence of a@ is predominantly due to the 
conductivity o4. 

The computations on the equatorial region analysis are currently in progress. Preliminary 
results seem to indicate that accounting for these regions has no effect on the depths of the 
modes, but will have some influence on the envelope of the spectrum. 

Other extensions which might reduce the discrepancies between the theory presented here 
and tie Balser and Wagner results are as follows: 

(1) The lightning statistics used here were obtained from data taken at many times of 
year and at many locations. Long term (e.g., year-round) averages of lightning parameters 
are not necessarily the same parameters averaged over a single day. The ELF noise results 
obtained by Balser and Wagner [1960] were obtained on a single day and at a single location. 
Only observed spectra averaged over many locations and times of year would be strictly com- 
parable with the theory discussed here, unless the ELF noise spectrum is found to be essentially 
independent of time and location, a supposition that would require many more experiments 
to ascertain. 
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Figure 6. ELF noise spectrum—effective ionospheric height 








(2) The same remarks apply to the assumption that the lightning which gives rise to 
ELF noise consists entirely of cloud-ground return strokes and that the lightning current is 
entirely vertical. Wait [1960a] has pointed out that there may be significant contributions 
from horizontal current components at the frequencies of interest here. 

(3) Assuming the return lightning stroke moment to be a random function and averaging 
over all random variables would probably make the lightning stroke more “‘noisy,”’ i.e., give it 
a flatter spectrum. This could result in a theoretical spectrum whose envelope bears a closer 
resemblance to that of the observed spectrum. 

(4) The ionosphere model could be improved by accounting for the earth’s magnetic field 
and the effects of heavy ions. 

(5) The probable reason for the fact that our theoretical curves are not as oscillatory as 
the experimental results is that the assumed ionospheric loss parameter is too large, tending 
to “flatten” the modes. This implies that the homogeneous ionsophere model is not sufficiently 
realistic for precise calculation of ELF noise effects. In a model accounting for the variation 
of conductivity with height near the lower edge of the ionosphere, e.g., the exponential model 
of J. R. Wait [1960b, 1960c], the loss parameter would be a complex number.’ This would 
result in a more oscillatory spectral curve, as indicated in recent work by J. Galejs [private 
communication], as vet unpublished, and in some preliminary calculations by the author. 
Work by the author. to infer information about the ionosphere from the Balser and Wagner 
curves with the assumption of a complex loss parameter is currently in progress. 


Appreciation is expressed to Dr. R. Row and Dr. J. Galejs for helpful discussions pertaining 
to the work reported here, and to Mrs. J. Van Horn and Mr. C. Roche for computational 
assistance. 


9. Appendix. Remarks on Analytical Points 


Note that with a uniform distribution, a source could be at the observer’s position. The 
discontinuity of the field at this point is accounted for in the fundamental theory by the fact 
that the hypergeometric series P,(—cos 6) used in the representation of the field has a discon- 
tinuity at @=0. The expansion of this function into a series of Legendre polynomials [Magnus 
and Oberhetinger, 1949], which is basic to our formulation, is valid everywhere but at 6=0. 
Sommerfeld [1949] has shown that P,(—cos @) becomes proportional to log 6 as 6 approaches 
zero. In carrying out the integration indicated in (7) in the vicinity of this point, the product 
(sin 6 log 6) appearing in the integrand is an indeterminate form which can be shown to vanish 
as 6-0. The contribution to the integral at 6=0 therefore vanishes, and no difficulties are 
encountered in extending the expansion to this point. 

Another analytical point worthy of note is that the expansion mentioned above is not 


valid at v=0, +1, +2, +3, ... ete. The point v=0 occurs at zero frequency, which we 
exclude from the analysis. Also, it can be shown that, except at w=0, v always has an imaginary 
part and therefore cannot take on values +1, +2, .. . ete. 
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Reflection and coupling processes exhibited by plane electromagnetic waves propagated 
in an inhomogeneous horizontally-stratified anisotropic ionosphere are associated with dis- 
crete transition points or with continuous coupling regions. These arise when the fourth 
order differential equations are written in first order coupled form, and many terms in these 
equations become infinite at the transition points. This procedure is rendered more precise 
by means of a special linear transformation that reformulates the equations in a new way, 
thereby exhibiting the manner in which local coupling processes are embedded in the more 
general background process of otherwise independently propagated characteristic waves. 
To exhibit the power of the matrix algebra involved, the case of an arbitrary number of 
characteristic waves is considered; moreover, Férsterling-type coupled equations are produced 
in a more generalized form than hitherto considered, and a discussion of the equations govern- 
ing continuous coupling completes the paper. 


1. Introduction 


The equations governing the propagation of plane electromagnetic waves in a horizontally- 
stratified inhomogeneous anisotropic ionosphere may be expressed in various ways, depending 
upon the particular investigation for which the equations are required. If the Oz-axis of a 
right-handed set of Cartesian axes Ox, Oy, Oz is directed vertically upwards, then it often proves 
convenient to use second or fourth order differential equations in the dependent variables E, or 
E,. Heading and Whipple [1952] and Heading [1955] have employed this representation when 
analytical solutions of the differential equations are under consideration. The reason is that 
all known analytical solutions of the field equations may be expressed in terms of generalized 
hypergeometric functions, hypergeometric functions, confluent hypergeometric functions, 
Bessel functions and so on, these functions naturally satisfying second or fourth order differen- 
tial equations. Such functions provide the reflection and conversion coefficients for a particular 
anisotropic model ionosphere, but by inspecting the functional form of the electric field through- 
out the medium no immediate physical interpretation is possible, whereby the actual process 
of wavereflection and wavecoupling may be explicitly exhibited. 

Sets of equations known as coupled equations have therefore been introduced, in order to 
enable approximate solutions for the electromagnetic field to be obtained, and in order to ex- 
hibit directly the important physical processes of wave reflection and wave coupling. Férster- 
ling [1942] considered a coupled system of two second order equations; these refer to an iono- 
spheric model in which the external magnetic field is oblique but in which propagation is vertical. 
Gibbons and Nertney [1951, 1952] have used these coupled equations to obtain approximate 
numerical solutions for various ionospheric models, while Budden [1952] has used them in dis- 
cussing limiting polarization. Heading [1953] and Budden and Clemmow [1957] have given 
similar coupled equations governing the case when incidence is oblique upon an ionosphere in 
which the magnetic field is vertical. 

When both the direction of propagation and of the magnetic field are oblique, four coupled 
equations are used, each of the first order. These were first introduced by Bremmer [1949], 
but more systematically using matric techniques by Clemmow and Heading [1954]. The most 
up-to-date account of these first order coupled equations is found in chapter 18 of the recent 
text by Budden [1961], but there the fundamental simplicity and beauty of the matrix formula- 
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tion is often disguised under a multitude of symbols and approximations. Moreover, certain 
important features of the matrices occurring in the coupled equations together with the possi- 
bility of nonsingular embedding do not as vet appear to have been recognized. 

The coupled equations exhibit those ranges of height in which the propagation of the 
four characteristic waves takes place almost independently. Any point at which many of the 
coefficients in the equations become singular denotes that coupling takes place there between 
the characteristic waves concerned, the other characteristic waves not embraced by the coup- 
ling still being propagated independently. 

The present investigation consists of an examination of the equations near the coupling 
points without the necessity of introducing approximations. It will be proved that the equa- 
tions governing the coupling process at a coupling point are embedded in the equations 
governing the propagation of the waves as a whole, such that all singularities in the coefficients 
of all equations are removed at the coupling point. Moreover, although ionospheric propa- 
gation is represented in terms of four first order equations, vielding four characteristic waves 
in a homogeneous or slowly-varying medium, generalization to wave processes governed by 
n first order simultaneous differential equations yielding n characteristic waves is easily 
possible using matrix notation. The sphere of usefulness of the present investigation will 
thereby be extended if n equations are considered. 


2. Transformation of an rth Order Linear Equation 


The preliminary study of a given linear differential equation rather than a set of first 
order simultaneous linear differential equations is necessary in order to define a principal 
embedded coupled matrix. 

The rth order linear equation in the dependent variable u 


wu =y,_1(z)u7-P t+e,_2(z)u7—-? +... .+01(2)u’ +r 9(2) u, (1) 


where bracketed superscripts denote differentiation with respect to z, and where the 7 
coefficients v%, %1, . . ., Y—-1 are given functions of 2, may be expressed as a linear system using 


matrix notation as follows: 




















ia ‘ m 1 0 0) er 0 0 > wa 4 
u’ 0 0 l 0 0 0 u’ 
d a 9) 
dz a ow eo ee oe ee ee ee ee ee ee a ’ (2) 
gy ?—2 0 0 0 0 eae 0 ] a? 
carrey At V; ls Ve 0 - arn 
or 
u=Tu (3) 


say, where u denotes the column consisting of 7 elements and T denotes the square 7X7 matrix 
shown consisting of the coefficients occurring in (1). 
When the matrix T is diagonalized by means of the transformation 
u=A f, 
the roots of the characteristic equation 
det(T—qI)=0 
are required. This equation has the explicit form 
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1—Y,9q°*—. ..—U1g—M=9, 


q—er-17~ 

possessing the 7 characteristic roots g1, q@2, . . .) Yr. 
The elements of the characteristic vector corresponding to the root g; consist of the 
cofactors of any row of the matrix T—gq;I; in particular, the cofactors of the last row are 1, 








93, V7, » + + Gj’ respectively. Hence the matrix A of the transformation may be taken to 
be 
@ ] ae Ps 
N q2 oe dr 
A= ; 
-_— °° fo 


this is an alternant matrix, whose determinant is given by 


: 
det A= II (q:—q,). (4) 
i, j=1 
i>j 
It follows that A is nonsingular throughout any domain in the complex z-plane in which all 
the r characteristic roots are distinct. 
Equation (3) becomes 


A’f+Af’=TAf, 


or 


f’=A-'TAf—A-'A’T. 


The matrix A7 T A is now a diagonal matrix Q, whose diagonal elements consist of the 
characteristic roots in order; hence 


f’—Qf—A-A'f. (5) 


If solutions of (5) are known, then u=A f are the solutions of (3), and in particular the solution 
of (1) is 


efi tfeb . 2. 4h (6) 


The r eqs (5) constitute a set of first order coupled equations equivalent to the original 
eq (1). The diagonal matrix Q will be called a principal uncoupled matrix, while the terms 
A! A’ f will be called the coupling terms, though strictly speaking coupling arises only through 
the nondiagonal elements of the matrix A~'A’. This product A7'A’, whose explicit elements 
are calculated in section 5, may be called the coupling matriz. It is the presence of this coupling 
matrix that renders the 7 eqs (5) simultaneous and not independent, but the coupling matrix 
vanishes if the coefficients in the original eq (1) are constants, owing to the presence of deriv- 
atives in every element. The 7 transformed equations are then completely independent, and 
their 7 solutions provide 7 characteristic waves independently propagated when a physical 
picture is employed for their interpretation. Then 


Y “4 § ~ 
S;=C; exp (qj2), 
where the C; are arbitrary constants of integration. 


When the coefficients in (1) are not constants, usually some parameter such as the frequency 
or some property of the medium attains values that cause the elements of the coupling matrix 
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to be small in magnitude compared with the elements of the principal uncoupled matrix, 
provided a domain is considered in which A is nonsingular so that its reciprocal exists. Under 
these circumstances, the nondiagonal elements in the coupling matrix are neglected, yielding r 
independent equations, whose solutions represent 7 characteristic waves propagated independ- 
ently through the slowly varying medium. 

In this paper, equations of the form (1) are considered only when all the 7 characteristic 
roots become equal at a particular point 2. To this end, r new functions a, B,... , p are 
introduced, defined by the r linear equations 


Qy=atkh? B49 Pyt+ get +hO-DUHD 5. (7) 
where k=exp (27i/r). The matrix of the coefficients in these equations is nonsingular, 
since it is an alternant whose elements consist of the rth roots of unity 1, k, k®, ... , k’" 
raised to the appropriate powers. The 7 new functions are therefore uniquely defined. 

Now 


Yri=Dy q3 
j=1 
=ratB 3 + y SBOP Et Dy Re Pen, 
j=1 j=1 j=1 


But if p is an integer, 


> me) (1k?) =1—(k?)'=1— (#)?=0, 
j=1 


so if pO, it follows that 


Tr 
Dy PY =, (8) 
j=1 
Hence 
a=v,_1/r. 


For example, when r=2, a and 6 are merely the two terms that arise when the quadratic 
equation is solved using the standard formula, the familiar + signs being associated with the 8. 
When r=3, 6 and ¥ arise when the reduced cubic is solved by Cardan’s method. 

To lay the basis for the following sections, v,_; is now chosen to vanish identically, indi- 
cating that a=0. It will be seen in the general theory, however, that no loss of generality 
occurs on account of this assumption. Coupling is then studied by postulating that the r 
characteristic roots become equal at z=z. Under these circumstances, when @ is omitted 
from (7) and when the r q’s are given equal non-zero values, the rank of the matrix of the coeffi- 
cients of the right-hand side is r—1, since the 7 linear expressions are dependent, while the 
rank of the augmented matrix is7. The equations are therefore inconsistent. If on the other 
hand the r equal values of the q’s are zero, the first y— 1 equations (7) yield B=y= . . . =p=0. 
since the matrix of the coefficients of these first y—1 equations is nonsingular. 

Hence, at 2=2Zp, it follows that 


UA=Qa= ... =g,=90; 


moreover, the coefficients of (1) in which v,_,;=0 consist of sums of products of these r roots, 


so 


U( Zo) =) (Zo) = . . . =V--2(20)=0. 
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This is too general for useful progress to be made, although the theory to be presented is 
equally valid without any further simplifying assumptions. The special choice is made that 


v,(Z) =02(z) =... =v,-2(z) =0, 
leaving a differential equation 
lu 
——-=09( 2). (9) 
dz 


When r=2, it should be noticed that no assumption is in fact made, and it is just this case 
that usually governs wave reflection and wave coupling. The r values of the characteristic 
roots are now given by 


g=E-"(04)"" kp, 


Another reason why these simplifying assumptions are advantageous concerns the possi- 
bility of writing down approximate solutions of (9), since (9) will form what may be termed 
a comparison equation embedded in an equation of higher order. Approximate solutions of 
this latter general equation are given in terms of approximate solutions of the comparison 
equation. When r=2, (9) can be solved approximately in terms of the Airy integral, while 
if r > 2, solutions may be obtained in terms of functions investigated by Heading in a series 
of papers [1957 a and b; 1960], though such approximations have not as yet been published. 
But if (1) is considered with v,_,=0 but with no further assumptions, useful approximate ana- 
lytical solutions are almost impossible to obtain near a point z=2) at which r characteristic 
roots become equal. 

Under these assumptions, matrix T becomes 




















(0 1 Se a a ee 
0 O 1 0 0 O 
pry Set aes etter Soaks PAu rants 2 ee ee (10) 
YF & © 0 1 
we + ee uw eB 
where vp is now written simply as v. 
The characteristic equation g’—v=0 has the r roots 
qy=ki- to" 
and the alternant matrix A factorizes thus: 
fl 0 0 oar “0 \ fl 1 1 rata l 7 
0 vw 0 eapoue 1 k k? si eke | i. 
A=|0 0 oe Ss, OU 1 k? ks ee  aitlias (11) 
L0 0 0 EE cia i de ne a 
=VK 


say, where V and K are defined by these two square matrices. 
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The matrix K being constant, it follows that 


A~'A’=K~'V-'V’K 




















@ 0 0 oem ” 0 0 <n 
0 pi" 0 0 vy’ /rv 0 
—K-! K 
0 0 yp 2/t 0 0 Qy?!"y’ /re 
CRESS SO St or dee Sa ee are J Beles Be cats hs ae ae Fae ia abirin ee) ges eta oe ate J 
i 0 0 ‘ 
P 0 1 0 
= = K 
ne 0 : had 
ree j 


Clearly coupling only becomes large near a point at which v’/v is singular. 
The explicit elements of this product have been calculated by Heading [1960]. They are 


(AA) y=, (12) 


27 


Se ee . 
(A A dij vp r(k?-'—1) (13) 


The coupled eqs (5) then take the form 














¢ 1 1 ‘ 
0 : cre 
k—1 P—1 

= aa ( 0) * 

1 1 
0 fey’? 0 ao ae ae ene 0 — se 

i ns ie f. (14) 

0 0 kyl" val ri ; :' 

Se ae ee a ee ee ee y, k-2—] k-2—] 0 

eee... 


The first square matrix on the right hand side of (14) is the principal uncoupled matrix. 
The coupling terms become singular whenever v’/v is singular; such a point is known as a 
reflection or coupling point. Mathematically speaking, such a point is called a transition 
point. The coupling terms are singular at this point, but the original equation 


u’=Tu 
possesses a special matrix (10), associated in this particular case with no additional coupling 


terms. An?r Xr matrix of the form (10) will be called a principal embedded coupled matrix. 


3. General System 


The equations governing ionospheric radio propagation may be expressed as four linear 
first order differential equations with four dependent variables, given for example by 
Clemmow and Heading [1954]. More generally, consider n dependent variables ¢,, @, .. . , 
é, satisfying the n linear first order equations 
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, la 
2 haa T',e; 


Boa 
In maxtrix notation, 
oe’ = Te, (15) 
where the nXn matrix T contains n? elements each of which is a given function of z. In 


the ionospheric case, n=4, and some of the elements of T are zero, while others are complicated 
functions of height, the explicit forms of which would obscure the essential simplicity of the 
present argument. The matrix T is therefore more complicated than that occurring in (3) or 
the principal embedded coupled matrix (10). 

If the dependent column e is now transformed to the column f by the transformation 


e=Sf, (16) 
it follows that 


f’=S-'Tsf—S-'8’f. 


If the n characteristic roots of T are denoted by q, qs, . - ., Gn, let 8; be the column matrix 
formed from the cofactors of T—q,I taken along any suitable row. Then the transformation 
matrix 

S=(Sis .. . 8,) 


diagonalizes the matrix T; that is, S"'T S=Q consists of the n roots arranged in order down 
the Jeading diagonal. The coupled equations become 


f’=Qf—(s\s....8,)7(siS....8,)’f 
=Qf—S"'S’f, (17) 
valid at all points at which S is nonsingular. 
If desired, a change of independent variable may easily be effected. If ¢=¢(2), then 


a @., 
dz d¢ vs 


vielding 
1 op _o- S 
Py Qf—s dé f. 


df 
dé 
Similarly, a change of dependent variables may be made by placing 
S5=O 2G 
where the n functions ¢;(z) denotes n given functions. In matrix form, let 
f=g, 


where @ is a diagonal matrix consisting of the n functions $,;(Z) arranged in order down the 


leading diagonal. Then 
bg’ + &’g—Qég—S—'S’ Og, 
or 
g’ = b-'Qbg— (#-'b’ + b-'S-18’5)g 
=Qg—(S#)~'(S#)’g, 
implying that these equations would have been produced had S# been used for the transforma- 
tion matrix (16) rather than S. The characteristic vector s; is merely multiplied throughout 
by ?;. 
602247615 601 








Now each cofactor of T—gq,l is a polynomial of degree n—1 at most in qj, the coefficients 
being independent of the suffix 7. Generally, 


(polynomial 1 of degree n—1 at most in q;) 
polynomial 2 of degree n—1 at most in q; 





(polynomial n of degree n—1 at most in q, ) 


(the n coefficients of polynomial 1) (1 > 


the n coefficients of polynomial 2 qj 


a= 


(the n coefficients of polynomial nJ (gj~") 














a. 
qi 
on P 
LG") 








sav, where P denotes the nn matrix formed from the coefficients of the n polynomials. The 
last column of P consists of n—1 zero elements, while the remaining elements is unity (when 
n is even), since only one polynomial can be of degree n—1, the others being of degree n—2. 


Then 


] 1 1 
q 2 ‘ ae 

S=(s\S2 see s,)=P i] 3 “A =PA, 
qi q2" q3" irae 


where A denotes the alternant matrix consisting of the roots and their powers arranged in order. 


It follows that 
S-'S’=A7~'P-'(PA’+P’A) 


==A-*A’+A~*P-"'P’ A. 
The coupled equations (17) take the form 
f’=Qf—A-'A’f—A-'P-'P’ Af, 


similar to the set (5) apart from the additional coupling terms A7!P-!P’Af. As before, Q 
denotes the principal uncoupled matrix, A~! A’ is the primary coupling matrix and A~ P=! P’ A 
is the secondary coupling matrix. 


4. Integral Equation 


Let the matrix A7'A’ be separated into the sum of two parts, D and (A™'A’—D) say, 
where D is an arbitrary diagonal matrix, though this may preferably consist of certain diagonal 
terms chosen from the explicit representation of A7'A’ given in the following section. Then 


f’=(Q—D)f—(A~'A’—D+A“'P7'P’ADf, 
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where the matrix Q—D is diagonal. 
Let f) be a solution of the equation 


f’=(Q—Dff, 


this equation representing in effect n independent linear equations. Let F, denote the nxn 
diagonal matrix containing the elements of f) arranged down the leading diagonal. Consider 
the integral equation 


f—Nf,—F, q F>1(A-'A’—D+A-'P-!P’ A)fdz, (19) 


where N denotes an arbitrary constant diagonal matrix and @ an arbitrary constant. 
This is the required integral equation, for differentiation yields 
f’=Nf,—F, | F,'(A7'A’—D+A7'P'P’ A)fdz—F,Fy1(A7'A’—D+A7'P~'P’ A)f 
a 
=Nf,— FoF '(Nf,—f) —(A-'A’—D+A7'P~'P’ A)f 
from (19). 


But 
Fo= (Q—D)F, 
since all matrices are diagonal, so 
f’=N(Q—D)f,— (Q— D) (Nf,—f)— (A“'A’— D+ A“'P“'P’ A)f 
= (Q—D)f—(A“'A’—D-+ A'P"'P’ADf, 
since N commutes with (Q—D). This demonstrates the equivalence of the differential and 
integral equations. 


5. Explicit Form of the Primary Coupling Matrix 


In order to appreciate more clearly the character of the coupling terms A7'A’f, it is neces- 
sary to evaluate the explicit form of the product A7'A’. 
The reciprocal of the nth order alternant matrix 








al 1 1 or 1 > 
en n Q2 q3 Qn 
— °° £° en: = 


may be calculated to have the form 


























{ (=D"S,-(H=0) (= D"-*S,-2(=0) =S(H=0) So, ’ 
(q1) (q1) ide (q) (q) 
(= 1)" *Sn-1(H=0) (= 1)" "8 n-2(=0) —Si(@=0) So 

(qo) (q2) sine (qo) (q2) q 
(=D""'Sr(G=0) (= 1)"*S 24 =0) =Si(q=0) Sy 

Ss (Qn) (Gn) maa (Gn) (Qn) a 
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Here, So(=1), Si, S2,. . ., S,-1 are the elementary symmetric functions of the n quantities 
M1; Jz, » » +» Un» The symbol S,(q¢;=0) is used to denote the algebraic expression S, when q; is 
replaced by zero. The symbol (q,) is used to denote the product 


(4p) =(Go— D1) Yo— G2) +» (Go Yo-1) (Yo Iv) + «+ (Ip In) = (OF 10) e= 2,5 (20) 
where 


F=det (qI—T). 


It can then be seen that 


(A-') ;=(—1) "IS ,.—)(Q:=0), (qi), 


and 
i-? 


(A’)3=((—1) Gi? q); 
so the general element in the product A~'A’ is given by 


(AA) y=D0 (A aA ay 
k=1 


=>) (—1)"*S,-2(q:=0) (k—1) g* -7q5/( 4). (21) 


k= 


In order to simplify this expression, the polynomial of degree n—1 is formed, whose zeros 


are Gi, Yzy + - +> In but with g; omitted. 
This is 
=q""'—S,(q.=0)g" +. . .+(— 1) "8, (qi=0) "P+ 
q—q; ‘ areee Sp\Fi sas 

n—-1 

=>) (—1)?S,(q:= aie 
p=0 
n 

=>> (—1)"*S,-2(4=0) 7", (with p=n—k) 
=I 


vielding upon differentiation with respect to q 


ra) i n , 
d, \q—g o> (—1)""*S,-1(Q:=0) (k-D gi? (22) 
a. * iy (= 


If +4), this result transforms (21) thus: 


ae fs 
(A A’) = qi) = qI- ==). aj 


. oF oat ‘.. 
_" q—G (9-4)? 


4 


__ G% (4) 9a) 
U—U (Qi) - 


since F' vanishes when q=q;. Out of the n—1 factors in (q;), the factor (q;—q;) obviously 


cancels with the factor (¢i—q;) occurring amongst the n—1 factors of (q;) 
On the other hand, when 7=/, the left-hand side of (22) may be written as 


be ss (¢—p)? 
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vielding 


(res, 
AA =| 2 TI (q— | 
( (qi) 07 p=! (q Up 7 


P ra) n | 
=; — log II (q—q,) 
ti >, v4 |, 
t 


2 , 
— a 

p=1 Yi—Yp 
i 


being the sum of n—1 partial fractions. 
6. Classification of Transition Points 


A specular transition point (a reflection or coupling point in physical terms) exists at a point 
2=2 in the complex plane at which two or more characteristic roots attain equality. 

A transition point of order two exists at a point where g,;=q2; such a point may be denoted 
by z=2.. The first and second diagonal elements each contain one partial fraction with a 
denominator (q:—q2), so this particular partial fraction (but no others down the leading 
digonal) is singular at the transition point. 

Equation (23) shows that every nondiagnonal element is singular at 2,. along rows one and 
two only. But only columns one and two have gj and gj respectively in their numerators; 
these derivatives may be singular at 2,2. There is, however, a difference in the order of magni- 
tude of the various coupling terms near the singularity 2,.. For the coupling coefficients 
(A7'A’) ;; ((=1, 2) are larger in order of magnitude when j= 1, 2 than when7=3, 4,...,. This 
implies that coupling exists mainly only between f, and f2, and that fs, fs, . . . ,f, are relatively 
free from coupling with f, and f,. In fact, the 2X2 matrix containing the first and second 
rows and the first and second columns is the distinct factor in determining the coupled solution 
near 2. 

If gq, and q» are expressed in the form 


“u=at+B, Qg=a—B, 


clearly B=0 at the transition point, while the first part a=(q,+q:) is a more irrelevant term 
in determining the character of the coupling at the transition point. 

Similarly, transition points of order two exist whenever q,—4q,; such points may be repre- 
sented by z=Z,,. Altogether, %n(n—1) sets of second order transition points exist in the 2- 
plane. The solution for f, and f, near 2,» is locally associated with a 22 matrix taken from 
the a- and bth rows and from the a- and 6th columns of A7'A’. 

The determination just considered of the important terms governing the coupling process 
is clearly unsatisfactory, since infinite terms occur throughout. The embedding transformation 
discussed in the next sections removes this uncertainty and yields a clearer picture of the 
coupling process than the matrix A~'A’, which merely shows where the coupling occurs. 

Every nth order equation must of necessity possess these specular transition points of order 
two, since they are merely given by the solutions of the well-defined equations g,=q,. But 
some equations may possess transition points of higher order, not as a matter of necessity but 
exceptionally. If a point exists at which qz=qo=4q:, the point 2,,, is called a third order transi- 
tion point. All elements in rows a, b, ¢ of the matrix A7'A’ are now singular at the point 2,9., 
but only the columns a, 6, c of these rows prove to be of importance, so the 33 matrix taken 
from these three rows and columns is the determining factor for the coupling beteen /,, fo, f.. 

Similar arguments may be used to define an rth (r<n) order transition point. Such a 
point exists if qa=Q=Q= ... at Z=Zan..., Where r of the characteristic roots attain 
equality at this point. The essential coupling between fa, fr, fc, . . . is determined by the 
rXr matrix selected from the appropriate rows and columns of A7!A’. 
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These 7 values of the roots may be expressed according to the scheme (7), but for the 
purpose of this paper it is assumed that gj=a+k?"'8. For the usu:l case when r=2, no assump- 
tion is thereby made. 

It should be noted that of the remaining n—r roots, s of these may also attain equality at 
2=Zane.-.. If their common value y is not equal to @ at 2,,,..., then this set of s variables is 
not coupled to the previous set of 7 variables. 

If various characteristic roots attain equality only as z-> for a certain range of arg 2, the 
coupling is described as continuous coupling rather than specular. 

Finally, irregular coupling may occur at points where the matrix P becomes singular, 
or at points where one or more of the characteristic roots become singular. It is not intended 
to discuss this form of coupling here, since it does not appear to be amenable to systematic 
treatment. 

7. Embedding Transformation 


Specular embedding is defined with respect to a particular transition point 2 of order rSn. 
The first order coupled eqs (18) in the n variables are to be transformed in such a way that 
the r variables associated with the transition point appear in the new equations with a principal 
embedded coupled matrix (similar to (10)), while the remaining n—r variables not associated 
with the transition point appear with a principal uncoupled matrix (a diagonal matrix), in 
such a way that all coupling terms in all equations are nonsingular at 2p. 

The independent propagation of waves not affected by the coupling is thereby exhibited, 
together with the explicit form of the rth order equation satisfied by the 7 coupled waves. This 
coupling is seen to be independent of the remaining n—r waves. 

The general coupled equations 


f’—Qf—A-'A’f—A-'P-'P’ Af 


are assumed to possess a transition point z at which q=q@= ...=q, If qj=at+k? 2, 
where k=e?*"", a and 6 being functions of z such that B=0 at 29, a new column g is defined such 


that 
(I, exp (fadz) oo 
oa OE NE TO Le 


_ eal 
say, where I, denotes the pp unit matrix. Then 


E’g+ Eg’ =QEg—A~'A’Eg— A~'P-'P’ AEg, 








or 
g’=E-'QEg—E~'E’g—E-'A~'A’Eg—E-'A-'P-'P’ AEg. (24) 
Now 
al, 0 
E~'QE—E~'E’=Q— 
0 I, 
e . = ws 6 oon SD 
0 kB oO 0 0 0 
0 0 F 0 0 0 
a 0 0 0 en k’-'80 re 0 
. 6 * Pes 0 ie si 0) 
: “! © 0 O QJ 
BK, 0 
a = 














2 EE TRI nt Tore 








say, where K, is the diagonal matrix with elements 1, k, k?, . . . , k7"! and Q,_, is the diagonal 
matrix with elements ¢,4;, . . - , Qn 
Equation (24) now becomes 


to 
Cr 


BK, 0 
g’ = EAA EAP PAB . (25) 
0 a 


Finally, a new column h is defined, satisfying 


g=R™'h, 
where 
(1 1 1 See | OF 9 
kp kB So chars 0 
B? kB? kp? Leas k27-D Be 0 


B" ~] ta Sama pse-Daer-1 a2. hoe v2gr-1 0 
U0 0 0 ae I,-+J 


B 0) 
=e (26) 
0 I, —T 


say, where B is an alternant matrix similar to (11), 








Equation (25) now becomes 


BK, 0 
aa R~'h—(E7'A7'A’E+E7'A7'P~'P’ AE) Rh, 
0 Q..; 
or 
BK, 0 
“=f ) RW RR-” REE AASB“ A PP AER (27) 
0 Q,,-; 


The original eq (15) is transformed into this form by the transformation 
e=Sf 
=PAf 
—PAEg 
=PAER“'h (28) 


In the next sections, it will be shown that (27) possesses the required embedding properties. 


8. Embedded Equation 


The coupling matrix occurring in (27) will be considered in two stages, firstly 


RR~'’+-RE~'A~'A’ER™ (29) 
and secondly 
RE-'A~'P“'P’ AER“. (30) 
Since R R™'=I, it follows that 
RR-’’=—R’R“. 


Moreover, E and R commute owing to the special forms of E and R; hence (29) becomes 
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—R’R-'+E-'RA~'A’R'E 
—=—E~'R’A~'AR~'E—E~'RA~'’AR“'E 
—=—E~'(R’A~'+RA~-")AR“'E 
=—E-'(RA“)’(AR~)E 
=E-'(AR~')-'(AR~')’E. (31) 


& 
5 


Also matrix (30) takes the form 
(AER~)-'P-'P’(AER~?). (32) 


At all points where coupling may be neglected (and this includes 2) since the next section 
will demonstrate that all coupling matrices are nonsingular at 2), equation (27) becomes 


BK, 0 
WR ( )eon 
0 Q,,-; 
B 0 BK, 0 B"! 0 i 
ast, adi, oP | 
0 = 0  -_ 0 = 


6BK,B™' 0 
3 )a. 
0 ax 


The elements of h are thus separated into two distinct parts. If the first 7 elements are 
denoted by h, and the remaining n—r elements by h,_,, it follows that 


i ~—s8£,.8"4,, (33) 
a hn (34) 
It may immediately be verified that | 
ak. A 4 (8 0 o .. (1 Bo! Bo? oP. 
1 8 kB ke... 0 kp Oo... 1 £"'s"' &-*8- 
8BK,B-'=- 
"re ke ke... Oe SR occ 1 k-2e-! &-4g-2 
ATS Sire om ea eee ‘ kc a ; Spe or ree ae ee or J { 




















io 5 
Oo 0 1 i 
ane (35) 
ae 








Hence the elements of h, are associated with an equation containing a principal embedded 
coupled matrix, while the elements of h,_, are associated with a principal uncoupled matrix. 
In fact, 
h{P=Ph=vh, (36) 
where v vanishes at 2. The separation of the principal terms has therefore been achieved. 
9. Nonsingular Coupling Terms 

Every element in the coupling matrices (31) and (32) must now be shown to be non- 

singular at 2. 
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The elements of A R™ occurring in (31) are considered. Since 











a 1 1 --. ner ince 
api}? & & cae a es 8 
r a Pe a ae Ean ee Ce ar een tea rr Ee 
: l k-Tt sf ie kor? grt 0 0 

= aoe oe e616 SSS, SOS OOS OS ee 
0 0) 0 ? 0 
U0 0 aa 0 Oe as ee 





it follows that the element (AR~!),, when 1<iSn, 15787, is given by 
am), =(e so ( 1 ‘) B75} 

oo) 

(ke-2)3-1 








Dt all ae 


=> i | Seat ile ial 
p=1 


: 
a (a+k?-18)'-1R- (@-D GD B-5+1 

p=1 ; 

tS 1-07 1)(j-1 1 
a, gl at pe Pee st 

p=1 1l=0 


i-—1 r 
==>) ap penne win yr. 


l1=0 Lp=! 


Now if (/—j+1) is a multiple of 7, 
r : 
Sh hP-DU-H4+D mp, 


p=1 


while if (/—j+1) is not a multiple of 7, this sum vanishes by (8). Hence the only terms in 8 
arising in r(AR™'),, are those of the form 6%" where N is an integer. Moreover, the lowest 
value of ((—j+1) is —r+1, so N can only equal 0, 1, 2,.... Hence 


(AR~) y= Darl J iat aid 


where /=Nr+j—1, where N is a non-negative integer chosen so that 
j-1sls<i-1. 
In other words, this element is a polynomial in 6’. 
It follows that (AR~');, contains 8 only through terms such as Nr Bhr-t9°. 


Finally, since (AR~)"!=(adj AR™')/(det AR™!) the elements of adj AR™ are further 


polynomials in 8’, while 
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det AR~'= (det A)/(det R) 


n 
= I (qi—)) Bit2---+-D det K 
i,j=1 


i>j 


from (4) and (26). Now when 7 and 7 are both less than or equal to 7, the relevant factors in 
the numerator are 


(Q:>—Gr-1) - - - (Qr—G1) (Qr-1 — Gr—2) - - » (Gra —- 1) « - » (G2—- G1) 


which is proportional to 


‘Cid «lia 5 oe BB. 
Hence det AR does not contain 8 as a factor, so the matrix AR™ is nonsingular at 2p. 


In the product (AR~)-! (AR™')’, each element contains terms of the form 6*¥7(6%")’, 
where 1420, N2O are integers. This is proportional to 6’—19’, with M20, N>1, the 
lowest power of 8 yielding 6’~'s’. If 8, when expanded in terms of z—Z) commences with a 
term (z—2)?, then 8’~'p’ contains the factor (g—2)”""'. This is nonsingular at 2) if brp—120; 
that is, if 621/r. This implies that v(z) =", when expanded as a power series in terms of 
Z—2, should commence with a term that is linear at least. 

Hence the elements of the primary coupling matrix (31) are nonsingular at z under these 
circumstances. The secondary coupling matrix (32) is also nonsingular at 2, since it has been 
shown that AR~ is nonsingular there. It is assumed of course that P is nonsingular at 2p 
in the definition of specular embedding. 


10. An Explicit Case 
To deal with an explicit case, consider a second order transition point associated with the 


third order differential equation 


, 


wu’? =ru’’ +ru’ + vu. 


u (’ 1 0 x 
a u’ I=] 0 0 1 wu’! 
dz 

u’’ L., v Vo = 


Equation (3) is 


The transformation 


u=| q2 3 f 
if & & 


vields equation (5), which, when expanded by the results of section 5, becomes 


a. By Sy S| 
q1 0 0 mi—fs Gia ie Wm. Aa -@ g—i1 
i Se ae a ee a ee 

‘ t—-h @&-O Gah © oS 2-2 
0 0 @& 








: ee @ , &@ % , & 
Lea ae a—-@' a—-& rege 
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If 2) is a point at which q;=q@, and if ¢=a+ 8, @=a—B8, the coupling matrix has the form 








4 =. bes ul a’ +B" «—f «8 qs qs 
qi— Ys 26 28 I—di Uu—g 28 
a ‘+g’ a’ +6’ a’ —p’ a’ —B’ is 
ee ere = +53 
28 "g—@ 2B 'e—G % a 
a’ +B" a’ +B" a'—B" eae Gh if | 
y @- oa G&—@G Ge —@% G&b—-h Gg—-Qd 
Singularities sonpndas on account of the terms 6’, 1/8, and ’/8. If 6? is approximately linear 
near 2, of the form z—2, B=(z—z)!, so the three singular terms are of the forms (z—2))7}, 


(z—Z)~# and sear e te cian tively. The 6’/8 terms are obviously of a higher order of magni- 
tude than the others; these give rise to the particular matrix 


r =—b 0 
(oi 
28 a LD SORES 
0) 0 0 


but other singalar elements arise in rows 1 and 2 and columns 1 and 2 
The embedding transformation (28), in which P=I, becomes 


u=AER~h 
1 1 1 exp (fadz) 0 0 [’ Bo! 0 
=; “Nn @& 4g 0 exp (fadz) O0];1 —s"! Ojh 
“Gf G&G & 0 0 e 2 
exp (fadz) 0 1 
= aexp (fadz) exp (fadz) q3 | h, 


(oe? +6?) exp (fadz) 2a exp (fadz) Gi 
while the coupled eqs (27) become 
0 1 0 
’=| 8 0 0 | h—E~'(AR~')~'(AR')’Eh 
0 0 4 


using results (31) and (35). The principal matrices yield approximately 


d a: (" (,) 
d2\y,) \e 0] \te 


h3= q3hs. 


exhibiting the embedding property. It should be pointed out that when approximate solutions 
are written down, it may be necessary to retain some of the coupling terms. For example, in 
the last equation, an additional term in A; must be retained on the right hand side if the W.K.B. 
approximation for h3 is to be obtained. 
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11. A Special Property of the Coupling Matrices 
Under the original transformation (16), the coupling matrix is 
A~'A’+A-'P-'P’A, (37) 
while under the embedding transformation (28), the coupling matrix is 
E~'(AR~')-"(AR~')E+ (AER~')~'P~'P’ (AER~'). (38) 


It is supposed that the last coupling matrix occurs when the first 7 elements of f or h are coupled 
together at Zp. 

However, the n—rXn—r square matrix formed by crossing out the first 7 rows and the 
first r columns of both (37) and (38) are identical. This may be simply demonstrated by 
partitioning all the matrices between rows 7 and r+1 and between columns 7 and 7+1, and 
expanding the products. 

Hence the n—r equations (r+1<j<n) in the elements of f 


F=f ;+coupling terms in f,,... , f,+-coupling terms in f,41, .. . fn 


and the n—r equations (r+1<7<n) in the elements of h 


hi=q,h;+coupling terms in h;,.. . , h,+coupling terms in h,41, .. . hy 
are identical as far as the coupling coefficients multiplying f,.;, . . . , f, and h,s:, 


h, respectively, are concerned. 
Hence any coupling diagonal terms not neglected in the equation for fj (r+1<j7Sn) are 
identical with the coupling diagonal terms not neglected in the equation for hj (r+1S5jSn). 


12. Férsterling-Type Coupled Equations 


As far as ionospheric radio propagation is concerned, two differential equations each of the 
second order may be coupled together. The original form given by Forsterling [1942] and that 
considered by Budden and Clemmow [1957] are specially simple, in that the characteristic 
equation for g reduces to a biquadratic. This special feature disguises the more general formu- 
lation that is possible. 

The general formulation of this new type of coupled equation is as follows. Let 2 be a 
transition point of order 7, namely a point at which q=q@= .. . =qg, Then the n coupled 
eqs (27) (with no singularities in the coupling coefficients at 2) may be re-expressed in the form 
of n—r+1 equations in n—r+1 variables, taking the form 


h‘? =6"h,+coupling terms 
h’,_,=Q,,_-h,_,+ coupling terms. (39) 


The coupling terms involve only the n—r+1 variables h, and h,_,. This result may be 
shown merely by eliminating /,, h., . . . , h,, from the last n—r+1 eqs (27), the equation 
involving h; being differentiated ry—1 times in the process. The coupling terms now involve 
derivatives of h,, but the coefficients themselves are also derivatives, which would be small in 
a slowly-varying medium. 

Similarly, if another transition point of order s (2Ss<n—r) occurs such that s of the q’s 
in Q,_, become equal, the number of equations may further be reduced. 

This process is explicitly exhibited for n=4, such that a transition point of order 2 occurs at 
2) at which g;=q2, and another transition point of order 2 at 2, at which q3=4q. 


The given equation 
e’=Te 
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is transformed into 
h’=(0 1 0 0°)h+coupling terms 


2 00 0 


0 0 gq 0 








10 00 q@ 
by the transformation 














e—PAER"h, 
where 
E=(exp ( fadz) 0 00}, B=fl LO Cy: 
0) exp (fadz) 0 0 B —B 0 0 
0 0 j ae 0 0 1 0 
4 0 0 a io + ey 
Similarly, if gz=y+6, gs=v—6, and if 
F=(1 0 0 0 \, S=f1 00 0% 
i 0) 0 Oo 1 0 0 
0 0 exp (Sydz) 0 00 1 1 
0 0 0 exp (fydz)J 10 0 6 —é& 














then the transformation 


e—PAER"FS-j 


(in which E, R7', F, S~! are mutually commutative) transforms the equation into 








fugo 1 6 eo terms (40) 
200 0 
a, 8 8 Ff 
10 0 & OJ 


where the coupling terms are nonsingular at 2 and 2;. 


m= ), n-("") 
J3 J 


Now let 


and two simultaneous matrix equations are formed involving m and n; that is, eq (40) is 


multiplied first by (f = a and then by > “ i) This vields 


0010 00 0 1 
m’=n+Wm- Xn, (41) 
2 0 
a’ = ° ) m+ Y¥m-}Zn 
0 & 
=Am-+Ym-+ Zn, (42) 
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say, where W, X, Y, Z are 22 matrices originating by partitioning the 44 coupling matrix 
in (40). 
Column m is now eliminated by differentiating (42): 
n’’=Am’+A’m+Y’m+Ym/’+Zn’+Z’n 

=(4+Y)m’+(A’+Y’)m+Zn’+Z’n 
=(A+Y)(n+Wm-+ Xn)+(A’+Y’)m+Zn’+Z’'n (from 41) 
=An+(AW+YW-+A’+Y’)m+Zn’+(AX+Y+YX+Z’)n 
=An+(AW+YW-+4’+Y’)(A+ Y)~'(n’—Zn)+ Zn’ + (AX+Y+YX+Z’)n 


from (42). Anis the only term on the right hand side whose coefficient does not involve a 
derivative. Hence A is the principal uncoupled matrix, while all the other terms are coupling 
terms. 

The coefficient of n’ is (AW+YW-+-A’+Y’)(4+Y)7". A subsidiary change of variable 
will eliminate nj from the first equation and nz from the second equation if required. 

The above process may obviously be carried out for any number of variables. 


13. Continuous Coupling 


Let the elements of the original matrix T (see eq (15)) be functions of z that tend to 
constant limits as, say, 2—-—© along the negative real axis. Then the n characteristic roots 
also tend to limiting values as z—>—o. 

If a range of real z existed such that all the elements of T were constants throughout that 
range, then all the coupling coefficients would be zero, in which case the n coupled eqs (18) 
reduce to 

f’—Qf, 
or 
A=, 


implying n-independently propagated waves. 
More generally, if a domain of z is chosen throughout which the coupling coefficients are 
small, the coupled equations reduce to 


S5=GF;—(A'A’ +47! PP’ A) 555;, 


where all nondiagonal coupling terms have been neglected; such approximate equations yield 
the W.K.B.-type solutions with arbitrary constants. These approximate solutions represent 
the n characteristic waves in the medium, and they are independently propagated whenever 
coupling may be neglected. In particular, only one characteristic wave may be considered, 


namely f,, with fA=fs= ... =f,=0. For this particular solution 
e=Sf, 
so 


C1: Cg. . oo 26 ,= Sy? Sn: ... Sn; 


thereby providing in this domain the varying ratios of the values of the n original dependent 
variables associated with this characteristic wave. 

The problem demanding attention concerns what happens to these ratios as z>—. 

The coupling coefficients may not be neglected when: 

i. coupling is specular; the previous embedding theory yields equations dealing with this 
case; 
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li. coupling is irregular; namely, when P~'P’ or when one of the q’s is singular. This 
paper is not concerned with this case; 
ili. coupling is continuous, namely when some of the q’s tend to equality as z>—o. 


To be definite, let q@, q@2,. . . ¢, have a common limit as z——. Then the denominators 
of many of the coupling coefficients considered in section 5 tend to zero. Thus continuous 
coupling between q, g2, . . . , g- is embedded in the background provided by the complete n 


equations. Previously, specular coupling has been explicitly exhibited without approximations. 
Here, however, it proves advantageous to use approximations in the equations. 

Moreover, it is assumed that the coupling coefficients do not become large in the range of 
continuous coupling. Every denominator of the form (q¢,—q;), 1S$7,7<7r, contains a derivative 
in its numerator, and as 2—>— © it is assumed that the derivative is such that the whole ratio 
remains finite at least. 

Above the region of continuous coupling, a definite approximate solution is considered, 
namely /;, with f=f;= . . . f,=0. Within the region, e=S f, with f,4:=f,.= ... =fr=0 
throughout. Then 


r r r 
€):€2. 2. . On, Suhit+2u Safi. .t 2 Snifi 
j= 


I " 
z r s r , ey \ 
=D Si Alf) > D5 Sa(filf) :...: DO Sas Alf) (43) 
j=l j=l j=l 
ratios depending upon the r—1 values f/f, fs/fi, . ~~ , tr/fi. 
These ratios are denoted by 72, 73, . . ., 7, respectively, and their initial values above the 


region of continuous coupling are all zero. 
The coupled equations 
f’=Qf-S"'S’f 


simplify to 


; 
f=ah-D (S'S )ufs (1<i<r) 
j=1 


in which the diagonal elements in S~'S’ may oftentimes be neglected since these would be small 
compared with the q;. 

A new set of r--1 equations may now be derived, with 72,73, . . .,7,as the dependent vari- 
ables. The differentiation of r; (2<i<r) yields 


nt (Lf _L5 
' az hi. 1 fi fi 
1 ! on , 1 " ‘C= , 
=p | ah (S7'S duds Jer 7 | wh (S-*S da | 
fi j=1 fi j=l 
ngs (S'S) yi  o-3 (S'S). 
j=l j=l 


being r—1 simultaneous nonlinear differential equations for 7,73, . . .,7,. These are the appro- 
priate differential equations governing continuous coupling. Integration starts with zero 


values, yielding To, 13, . . +) Tp BS SPO, which yield in turn the required ratios (43 ) of the e’s 
associated with this characteristic wave. In the region of continuous coupling, the original 
ratios S,,:So: . . .:S,, are converted into the final values. 
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Recent studies and measurements have shown that the phase stability of the signals 


in the VLF region is very much higher than in the HF spectrum. 
has caused considerable interest in employing this med- 


excellent coverage characteristics, 


ium for the wide distribution of standard frequencies and time reference. 


This fact, along with its 


Basic limitations 


in stability of the received signals are discussed, including path phase distortion, carrier-to- 
noise and envelope delay variations as related to precise synchronization of clocks, and 


highly accurate frequency calibrations. 


Also included is a discussion of the present services of standard frequency and time sig- 


nal stations throughout the world at HF, 
1. Introduction 


The requirement for better standards for precise 
measurements of frequency and time which are 
readily available has constantly increased with 
advances in various fields of science. The scientific 
unit of time is determined astronomically, and is 
determined to an uncertainty of a few parts in 10° 
in the course of a year. Since frequency is related 
inversely to time interval, this uncertainty, strictly 
speaking, must be transferred to the specification of 
absolute frequency. Frequency generators, how- 
ever, stable to parts in 107° or 101! now exist, and 
form the basis for relative measurements of frequency 
which for many purposes satisfy the increased re- 
quirements mentioned above. 

The history of the development of our present 
system of time is an extremely fascinating subject 
and it is interesting to note Newton’s definition of 
time, “Absolute, true and mathematical time, of 
itself, and by its own nature, flows uniformly on, 
without regard to anything external,” [Mach, 1942]. 
Most of man’s time systems have been based on an 
attempt to find a reference which “flows uniformly 
on.’”’ The very great influence upon our lives and 
surroundings caused by the revolution of the earth 
around the sun and the rotation of the earth on its 
axis has resulted in various definitions * of time in 
terms of years (orbital rotation), days (axial rota- 
tion), and in of these units in hours, min- 
utes, and seconds. Time based on apparent solar 
days (from noon to noon) although apparently 
satisfactory as a time reference for many applica- 
tions was shown by astronomical observations to 
have appreciable nonuniformity. 


! Present address of A. D. Watt and R. W. Plush: DECO Electronics, 
Inc., Boulder Division, 8401 Baseline Road, Boulder, Colo. 

2 An interesting and useful description of time and time signals is contained in 
“United States Naval Observatory Circular No. 49,” published by the U.S. 
Naval Observatory, Washington 25, D.C. (March 8, 1954). See also “Astro- 
nomical Time,’’ G. M. Clemence, Rev. Mod. Phys. 29, 1 (January 1957) and 
“The Way Things Are,” P. W. Bridgman, pp. 135-141 (Harvard University 
Press, Cambridge, Mass., 1959), 





LF, and VLF. 


Since apparent solar days are variable in length 
with season, mean solar time was devised and the 
tabulation of this difference between the apparent 
value and the mean value is called the equation 
of time. This difference has a maximum value a 


little in excess of 16 min. Apparent solar days are 
variable in length with season for two reasons: 
partly because of the variation in the angular 


velocity of the sun along the ecliptic in accordance 
with Kepler’s second law of planetary motion, and 
partly because the inclination of the ecliptic to the 
celestial equator introduces a variation in the rate 
at which the projected coordinate of the apparent 
sun moves along the celestial equator. The mean 
sun is a fictitious body which moves along the 
celestial equator at a uniform rate equal to the 
mean rate of the apparent sun. Because of the 
fact that a solar time reference (date and time of 
day, etc.) will depend upon the observer’s position, 
Universal Time (also known as Gre wich Mean 
time) was established with the mean solar time 
reference based on the prime meridian at Greenwich, 
England. The problem of relative position for civil 
purposes has been solved by dividing the earth into 
various time zones usually differing by one hour 
with lines of demarcation chosen to minimize the 
inconvenience caused by crossing time zones in 
heavily populated areas. 

Three kinds of Universal Times are: UTO, which 
is uncorrected mean solar time; UT1, is Universal 
Time corrected for observed polar motion and repre- 
sents, the true angular rotation of the earth about its 
axis;> UT2, is Universal Time corrected for both 
aed polar motion and seasonal variation in 
speed of rotation of the earth on an extrapolated 
basis. 

The fundamental unit of time was until very 
— the mean solar second which was defined as 


of a mean solar day. This fundamental unit 


86: = 


3 Because of the increased precision with which the earth’s rotation can be 
measured with respect to a reference star rather than the sun, the period of the 
earth’s rotation with respect to the vernal equinox, called a sidereal day, is used 
for precise determination of the earth’s rotation. 
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of time based upon the mean solar day appeared to 
be entirely satisfactory until astronomical observa- 
tions and advances in the fields of communication, 
electronic frequency standards, and precision in- 
strumentation showed that in fact the earth was not 
rotating in a strictly constant manner. Even after 
the corrections of UT1 and UT2, time based on this 
reference did not flow on uniformly as our definition 
of it would require. A more uniform astronomical 
time based on the vearly motion of the earth about 
the sun, called Ephemeris Time, has long been 
known. The present definition of the second is 
] 
31,556 ,925.9747 
1900 at 12 hr Ephemeris Time [Markowitz, 1959]. 
This was adopted by the International Committee 
of Weights and Measures in 1956 and ratified by 
the 11th General Conference on Weights and Meas- 
ures in 1960. Discussion of various astronomical 
and atomic times has been given by Markowitz 
[1959]. 

Although quartz crystal clocks were instrumental 
in revealing many of the variations in rotation of 
the earth, it was not until atomic standards were 
placed in use that a new reference of time became 
conceivable. The use of molecular and atomic 
spectral lines as frequency standards [Lyons, 1952] 
has led to the speculation that perhaps here we have 
a reference which can be made highly independent 
of external variable influences and as a result may 
produce a time base of greatly improved uniformity 
and convenience. The difficulty of keeping molec- 
ular and atomic clocks in operation for long periods 
of time has up until recently prevented their use 
as a basic time reference. A determination of the 
atomic transition frequency of cesium in terms of 
Ephemeris Time by Markowitz, Hall, Essen, and 
Parry [1958] indicated a value of 9,192,631,770 
+20 c/s which is the best value available at present. 
The deceleration in the rotation of the earth about 
its axis as determined by the cesium standard is 
in the order of 1.1610~° radians/vear’, i.e., 0.16 
seconds/year?*, which is in good agreement with 
values obtained by the moon camera of 0.17 seconds/ 
year? for that period (June 1955 to June 1958). 
It is possible [Bullard, 1955] that Ephemeris Time 
and Atomic Time may not have the same rates; 
however, further careful experiments will be re- 
quired to determine this. The stability of the 
best long time interval bases, which at present are 
quartz clocks steered by atomic standards is limited 
to values in the order of one part in 10. Frequency 
comparisons quoted to higher precisions are relative, 
and may be made in two ways. One is to intercom- 
pare two oscillators directly for the short term, 
and the other is to compare one with the mean of 
a group of three or more of the same quality and 
assume that their variations are independent, so 
obtaining a statistical increase in precision for the 
group. The resonance of cesium is known in terms 
of the Ephemeris Second to an accuracy of only 
about 2 parts in 10°, the limitation being set by 


~ 4 This means that if an atomic clock and “earth” clock were started together, 
at the end of a year the ‘‘earth’”’ clock would be 0.08 sec behind the atomic clock. 


of the tropical vear for January 0, 


| 





the uncertainty in the determination of the Ephe- 
meris Second. Thus, high quality clocks closely 
steered by atomic standards will not realize the 
presently defined unit of time with any greater 


; accuracy than this because of the nature of the 


definition itself. Presently, maximum available 
relative frequency stabilities for periods of several 
hours are in the order of: 1 part in 10" for Crystal 
standards and 1 part in 10” for the Ammonia Maser. 


2. Basic Limitations of Precision in the Fre- 
quency and Time Standards Via Radio 
Signals 


With the availability of very stable time references 
and frequency standards, it is obvious that con- 
siderable care must be exercised if these standards 
are to be distributed over large areas of the world 
without appreciable deterioration. If the terrestrial 
propagation of radio wave energy was at a constant 
velocity equal to the velocity of light in a vacuum 
and with a noise free background, it would be a 
relatively simple matter to distribute frequency 
and time signals with essentially no loss in stability 
and accuracy. Under these conditions, the received 
frequency would be unaffected and the time signals 
would have a given correctable constant delay of 
T’,=d/v) where: d is the distance and v is the 
velocity of light ~3 10° km/s. In practice the 
radio energy is transmitted along a path with an 
effective velocity (v) which varies with the charac- 
teristics of and conditions along the transmission 
path so that the actual delay, 7;, may vary also. 
It is obvious that the time indicated at a given 
receiving site is 7,=7,+7,, where the receivin 
clock can be made to indicate transmitted time 7’, 
by subtracting 7, from the received time 7, during 
clock setting. At the high frequencies, 2.5 to 30 
Me/s, presently employed for the distribution of 
frequency and time signals, multiple propagation 
paths frequently exist with appreciably different 
delays which can introduce large apparent errors 
in the received frequency and time. For any given 
path it may be possible by the proper choice of 
frequencies to have essentially one dominant mode 
of propagation; however, even a single mode may 
have appreciable variations in transmission time. 


2.1. Phase Distortion Limitations 


In figure 1 the straight dashed line represents 
the distortionless transmission in free space where 
the phase delay in radians could be given as 


or =Pald 
0 


(1) 





where: f is the frequency in cycles per second, and 
d is the path length expressed in kilometers. 

The actual propagation medium, assuming one 
mode of multiple path to be dominant, may be 
considered as an electrical network with phase and 
amplitude versus frequency characteristics which 
vary as a function of time. The actual time 
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Figure 1. Transmission path phase characteristics. 


delay is 
Ta=d/v, (2) 


where v, is the average path phase velocity in 
kilometers per second. It is well known in electrical 
circuit theory that various types of network delays 
exist, and if we are to consider the phase delay 
defined in eq (1), v7 must be expressed in terms of 
path phase velocity rather than envelope velocity. 

The actual path phase velocity relative to the 
velocity of light, v,/v, is known to be a function 
of frequency and the path involved. Important 
factors may include ionospheric conditions, ground 
conductivity, and surface roughness. In the VLF 
region, Jean, Taylor, and Wait [1960] have shown 
good agreement between experimentally and theo- 
retically determined values of v, over the 4 to 20 
ke/s frequency range. (See also Wait [1961a] and 
Wait and Spies [1961].) At night they have found 
values of v,/v) to be ~1.03 at 4 ke/s, 1.01 at 10 ke/s, 
and 1.003 at 20 ke/s. Many references pertinent 
to the velocity of propagation are contained in 
the paper by Jean, Taylor, and Wait [1960]. In 
general, v, increases with surface conductivity, 
and decreases with increase in ionospheric height 
and earth’s surface roughness. 

An actual physical path phase characteristic may 
look somewhat like the solid wavy curve of figure 1 
where a frequency f, would be delayed by the phase 
¢;. This phase delay in seconds can be written 
as 





=¢1 ‘ 
Ta Qnf; (3) 
where 7; is the equivalent time delay in seconds. It 
should be pointed out that the particular phase 
versus frequency curve illustrated will only apply 
for a given distance and at a given time and that 
it may vary in its exact position as a function of 
time. The effect of these short term variations 
of phase as a function of time is to limit the precision 
with which a given standard frequency can be 
received. 
To obtain the value of the received frequency 
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Ficure 2. Short term propagation time (phase) stability day- 
time paths, normalized to single reflection sky wave. 
from the phase measurements it is clear from the 
definition of frequency that the average departure 
from the transmitted frequency, Af, observed in 
terms of the final and initial phase delays is 

given by: 
— ¢(T)—¢4(0) 
Af ont (4) 
where T is the time interval, in seconds, over which 
the phase comparison is made. 

The relations governing the transmission path 
phase stability and signal-integration times required 
for specific frequency-comparison precisions have 
been described previously in considerable detail 
[Watt and Plush, 1959]. In general, where possible, 
observations should be made near noon at the 
center of the path. At such time, the phase varia- 
tions, due to diurnal and random short term effects, 
are expected to be minimal. Pierce [1957] has 
experimental evidence over a 5200 km path that 
the total diurnal phase variation of the 16 ke/s 
Rugby signal is surprisingly constant throughout 
the year being in the order of 34 +1 usec. Jean? 
over a 7500 km path has observed that the diurnal 
variation of the phase of the 16 ke/s Rugby signal 
has substantial variations in pattern. The observed 
average phase change of ~42 usec yields about the 
same effective change in ionospheric height; 1.e., 
Ah ~16 km using a method described by Wait 
[1959]. The diurnal patterns observed by Jean 


5 A, G. Jean, private communication. See also C. J. Chilton [1961]. 
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appear to change systematically with seasons and 
appreciably particularly during ionospherically dis- 
turbed conditions. The daytime transmission delay 
appears to be much less affected by these disturbances 
than the nighttime phase. 

The short time variations in propagation time 
resulting from path phase instabilities can, in 
general, be assumed to be_randomly distributed 
about an average time delay 7; with a rate of change 
similar to the fade rate. This means that a single 
measured value of the delay ranges roughly between 
T, + o (T,), where o (T,) is the standard deviation 
of the random time delay variable. In general, 
there will be different values for both 7’, and o (7) 
for day or night paths. Expected standard devia- 
tions of short time propagation time delays are 
shown in figure 2 as a function of frequency for 
a single reflection daytime skywave in the VLF and 
LF regions. It should be emphasized that this 
trend does not continue up into the HF region. 
The data are obtained from table B-1 [Watt and 
Plush, 1959] with the additional point at 10.2 ke/s 
from Tibbals.® 


2.2. Time Reference and Clock Setting 


It is obvious that if a single precisely known 
uniform frequency were available at a receiving 
location, that it could be employed in the control 
of a time standard clock, provided: (a) that the 
clock could be properly set initially, and (b) that 
there were no interruptions in the frequency standard 
or clock. Since neither of these can be guaranteed, 
it becomes apparent that some method of establish- 
ing a time or phase reference at the receiving 
location is necessary, and that the problem of trans- 
mitting a precise time reference is much more 


6 M. L. Tibbals, private communication. 
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Figure 3. Standard frequency and time signal waveforms. 





difficult than a standard frequency transmission. 
At the higher frequencies where bandwidth limita- 
tions are not severe, this is accomplished by trans- 
mitting rather short pulses whose leading edges at 
the transmitter are quite steep. The delay of these 
time pulses is also variable but it is not the same as 
the phase delay described in eqs (2) and (3). Since 
the pulses are applied in the form of modulation on 
the rf carrier. the “envelope delay’ as defined by 
Nyquist and Brand [1930] in terms of the slope of 
the phase characteristic, d¢/df, is the important 
factor. As we will see in the next section, the HF 
paths are very instable in both phase and envelope 
time delay and the timing accuracy obtainable with 
these pulses is considerably less than that available 
at the transmitter. 

In view of the much more stable phase delay 
conditions in the VLF region as described by 
Pierce [1957] and Allan, Crombie, and Penton [1956] 
along with its excellent coverage characteristics, 
considerable interest has been exhibited in employing 
this medium for the distribution of standard fre- 
quencies and time reference. The channel band- 
widths available in this frequency region, along with 
the very narrow bandwidths of the transmitting 
antenna, prevent the use of short-time pulses for 
the establishment of coarse time markers which 
could be used in identifying a particular rf cycle. 
An alternate method employing two closely spaced 
and alternately transmitted related frequencies is 
being considered. A description of this method is 
given by Morgan [1961]. The transmitted and 
received waveforms for such a narrow band time 
signal system are shown in figure 3 where it is 
apparent that the envelope of the voltage difference 
produces time reference markers with ambiguities 
spaced by the reciprocal of the frequency difference. 
Coarse time markers may, if desired, be provided 
by the times of frequency switching. On-off keying 
can also be considered for date time code. In a 
dual frequency system. transmitted waveforms are 
given as: 





€;=Ccos wt (5) 
€2=COS wot (6) 


while the envelope of the phasor sum is given as 





E(t) eavetope= 2+2 cos (w2—« )t (7) 


It can be observed that at the transmitter, zero 
time is our phase reference. At the receiver, the 
phases of e; and e, are delayed and these two signals 
are 

€;,=CO0s [w,t+¢)] (8) 

€2,=CO8 [wot +¢2] (9) 
Since these two signals can be employed to phase 
lock high quality local oscillators, it is possible to 
have two highly constant amplitude sinusoids as 
indicated. The envelope voltage waveform at the 
receiver is then 


|< | | received —* 2 +2 cos [ (w2—w )t + (¢.—¢, ) } ( 10) 


where it is evident that the important phase is now 
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—q, and the resulting time delay is 


g2—1 
RGF) 


where 7; is the envelope delay in seconds. 

If the individual frequencies concerned are sepa- 
rated far enough, the path phase fluctuations 
become independent and the standard deviation of 
T, would be 





(11) 


T,)=- y2o(¢) | 
° aa = $5) 
assuming that o(¢:)=o(¢2) =a(¢). 

It can be seen that the time variations will decrease 

with an increase in frequency separation. It should 
be emphasized that as the frequency separation is 
reduced the relative phases tend to become corre- 
lated and as a result the envelope reference time 
variation does not increase without limit since the 
phase difference tends to approach zero. Measured 
standard deviations of envelope delay variations 
for a 1 ke/s frequency spacing in the VLF band are 
shown by Casselman, Heritage, and Tibbals [1959] 
to be in the order of +20 usec. (See also Stone, 
Markowitz, and Hall, 1960.) 7 

Employing the ionospheric roughness parameter 
described by Watt and Plush [1959] the standard 
deviation of short term phase variation; i.e., o(@), 
of a VLF carrier received over a 4200 km path is 
found to be in the order of 0.09 radians. When this 
value is used in (12) the expected o(7)) expressed 
in microseconds is shown in figure 4. Experimental 
values for o(7,) obtained from Heritage and 
Tibbals § for spacing of 200, 1.000, and 3,000 e/s 
appear to substantiate the general level and trend 
of this curve. 

In order to permit identification of a specific cycle 
of the carrier frequency f;, it is obvious that o(T)) 
must be less than 0.5//;—o(7,) which at 20 ke/s will 
be in the order of +22 usec. Until such time as the 
expected reduction in o(7,) due to integration for 
time intervals long compared to the fade periods is 
obtained, it appears that frequency separations in 
the order of 1 ke/s will be required for carrier cycle 
identification in the 20 ke/s region. If the fre- 
quency separation chosen is rather large, the 
(fo—f;) markers may be difficult to resolve from the 
on-off keying envelope. This is due in part to the 
normal build up time required by the narrow an- 
tenna bandwidths, see typical shapes in figure 5, 
and the variation in position of this envelope as 
received due to path envelope delay. Should the 
on-off keying phase markers received with time 
delay variations o(7.) not be stable enough to re- 
solve (f2.—f,) markers, it may be necessary to employ 
an additional frequency ie an intermediate 
reference, viz (fs—fi)< Gf2—f). 

Since the precision of frequency comparison _ is 
defined as the phase jitter divided by the observing 


a( (12) 


7 Stone, Markowitz, and Hall in a recent paper [1960] have indicated VLF 
time comparisons of 500 wsec at Washington, D.C., from NBA in Panama 
Employing a 6 db transmitting bandwidth of 44c/s on fig. 4 and assuming that the 
envelope matching would approach that of two frequencies spaced by the 6 db 
bandwidth we obtain o(7)™~450 usec for a path of comparable length. Actually 
the spacing chosen should likely be less which would yield better agreement. 

8 Private communication. 
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period 7'(7.3 of Watt and Plush, 1959) we can write 
e=4 2o($)/2nfT. (13) 


The time delay for a given phase delay is 7,=¢/2zf, 
and since o(7',)=o(¢)/2zf, we can write (13) as 
d i . 


o(T,)=eT/y2. (14) 


Figure 6 shows the manner in which this relationship 
can be employed to determine the variations expected 
in a given clock relative to a standard clock T 
seconds after it has been set to coincidence with the 
standard clock. It is interesting to note that if an 
all daylight radio path such as the 16 ke/s trans- 
atlantic path observed by Pierce [1957] is employed, 
the maximum standard deviation of time difference 
is less than 2 usec. 

It should be pointed out that the problem of sub- 


tracting out the mean true transmission time, 7%, 
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Ficure 4. VLF transmission path envelope phase¥stability. 
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Figure 6. Clock-time deviations from original setting as a 
function of observing period. (8 Pierce 1957). 


whether it be a phase delay or envelope delay, still 
exists. lonospkeric disturbances caused by solar 
related proton storms or meteoric showers have the 
largest apparent effect on the nighttime phase delay 
(Chilton, 1961]. As a result, the midday values of 


T, are expected to be relatively constant from day 
to day with the exception that solar flares effect the 
day values. Over some paths, it is possible that 
there may be appreciable seasonal variations in 
even the midday value of 7;. Additional studies 
are required before the errors likely in predicting 
T, can be specified as closely as ‘desired. If two or 
possibly three VLF standard frequency and time 
signal stations are operating simultaneously at 
different rf frequencies but with a common fre- 
quency base as well as a common time reference, it 
is expected that the average delays due to transmis- 
sion time can be obtained with greater accuracy. 
Once the average path delays have been subtracted 
out, a clock at a given receiving site can be set with 
a high degree of precision (probably limited mainly 
by the ability to determine absolute average path 
delay), to the time generated by the network of 
VLF standard frequency and time signal broad- 
casting stations. 


2.3. Carrier to Noise Limitations 


The limitations of noise to precision of frequency 
comparison have been previously described [Watt 
and Plush, 1959] and we shall only calculate the 
effects of rms carrier to rms thermal noise ratio on 
time signal accuracy. 

From equation (C1) of Watt and Plush [1959] 
and (3) 

T/ Y 
o( Te) = (15) 


(noise) 242m; 





where o(7,)(noise) is now the standard deviation 
of the time reference in seconds derived from the 
carrier frequency (f,) phase, and N/C is the rms 
noise to carrier ratio in the receiver. With post 
detection filtering or integration following a linear 
detector, N should be calculated in the effective 
bandwidth of the integrator. If we consider noise 
density or in particular, the noise voltage per 1 ke/s 
bandwidth concept where N,x- has the dimensions 


volts/ ke/s, we obtain 


,)___vB,11,000 
9 ,)-—_>Be11,000 ay 
(noise) 24 2rfi( / N lke 


where B, is the effective receiver bandwidth in 
cycles per second. 

When timing is derived by an envelope beat pat- 
tern between two carriers, the effects of noise are 
given by 


NIC 


— 17° 
o(T) 2rr( fo—fi) ( a 

or 
WT) a3 BL, 000 ie 


(noise) 27 ( f, —fi )¢ YN ke 


where it is obvious that the time stability for these 
conditions is poorer for a given C/N; xe/s ratio than 
obtains for ¢(7,) (noise). 

We can olso determine the stability of a time 
reference obtained from a keyed envelope obtained 
at the one half amplitude point of the leading edge 
with coherent detection. For a single sample this 
is shown in the appendix to be 


i 1 : 
o(T.)=(CIN) OB, 


(noise, sampled) 


(19) 


where: o(7,) is in seconds and B, is the transmitter 
6 db bandwidth in e/s. If the carrier is keyed on 
periodically (for example, once each second) the 
resulting envelope reference can be integrated to 
reduce the effects of noise. The improvement is 
equal to the square root of the number of individual 
pulses, and if the repetition rate is F 


1 


Tj (20) 
(C/N)B ATF 
(noise, integrated) 
or ais 
pV B,/1,0007 F 9 
o(T)="(CIN, 1.) B, (21) 


where: 7 is the integration period, 7'F is the number 
of samples averaged, B, is the effective receiver 
bandwidth in ec/s. It is apparent from (21) that 
the time reference derived from the envelope 
leading edge improves directly as the carrier to 
noise ratio and as the first power of the transmitting 
antenna bandwidth and as the square root of the 
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integration time. In all the derivations involving 
the rms carrier to noise ratios, the significant noise 
background considered is of thermal type. The 
actual background will normally be atmospheric 
noise which in a wide bandwidth has a very different 
statistical nature [Watt and Maxwell, 1957]. The 
effects of atmospheric noise upon radio systems are 
in general quite different from those of thermal 
noise [Watt, Coon, Maxwell, and Plush, 1958]. For 
a pulse sampling system with integration, the re- 
quired rms C/N may be appreciably less than indi- 
cated in eqs. (20) and (21). Hefley® [1960] for 
example has found with Loran-C systems that for 
equivalent performance the carrier to atmospheric 
noise ratio can be much less than the carrier to 
thermal noise ratio; however, it must be remembered 
that this system employs a very wide RF band- 
width. ; 

When carrier phase or frequency stability is con- 
sidered, the narrow bandwidths (0.1 ¢/s or less) 
employed will make the atmospheric noise statistics 
the same as thermal noise [Watt and Maxwell, 1957] 
and in these cases the formulas given (15) through 
(18) should apply if linear receivers are employed. 


3. Present Services 


There are at present quite a large number of 
standard-frequency and time-signal stations in use 
throughout the world. VLF, LF, and HF are em- 
ployed for large area coverage with some higher 
frequencies being used for local distribution. A 
description of the stations operating along with 
their location and characteristics is given in Annex 
I of CCIR Report 66 [1956]. 

Annex II of the same document describes the 
characteristics of projected standard frequency and 
time signal stations at four additional locations. In 
general, the high frequency stations all employ pulse 
modulation of the carrier frequency for time signal 
transmissions and many of them in addition have 
standard audiofrequencies employed as modulation 
at various times. The stability as broadcast is 
usually very good, for example with the National 
Bureau of Standards station WWYV [1960a] the carrier 
and audiotones have stability of one part in 10° at 
all times with normal daily deviations of less than 
two parts in 10°. The time signals have essentially 
the same stability and the maximum deviation 
from UT-2 is about +50 msec. Final corrections to 
the frequency as broadcast versus the U.S. Frequen- 
cy Standard are published monthly by the National 
Bureau of Standards [1960b] in the Proceedings of 
the IRE while final corrections to the time signals 
versus U'T-2 time are published by the U.S. Naval 
Observatory. 

3.1. HF Standard Frequency and Time Services 

Because of the variation in the propagation 
medium at HF, the frequency as received is gener- 
ally much less stable than that transmitted. Ac- 
tually, for high precision measurements it is 
frequently necessary to use long averaging periods 





*Private communication, 








of 10 to 30 days. Ideally, the following conditions 
should prevail at the time measurements are being 
made; (a) all daylight (or darkness) over radio path, 
(b) no ionospheric disturbance in progress, (¢) no 
part of the radio path should pass near either 
auroral zone, (d) single mode propagation should be 
existing. In general, frequency determinations can 
be made either by a direct comparison or by means 
of the time pulses. Some idea of the precision of 
frequency comparison obtainable can be seen from 
figure 7 which shows the precisions available from 
WWV and WWVH over paths ranging from approx- 
imately 2,000 to 8,000 km. Line A shows that for 
short period samples, such as 100 sec, the precision 
available is in the order of 2 parts in 10’. If similar 
100-see samples are made each day, after 30 days 
the precision has improved to one part in 10%. It 
should be emphasized that this type of frequency 
comparison does not make full use of the transmis- 
sions and that the improvement with observing 
period is not as great as would be predicted if all 
the actual information received was fully utilized 
for the total periods shown. Line B shows the 
results obtained at Boulder employing time pulses 
compared for approximately 10 sec near aoon each 
day and averaged by means of a clock. Here the 
daily averaging time is 24 hr. The point shown for 
30 days is arrived at as follows. Daily comparisons 
of the local clock versus the received time of WWV 
are used to obtain running 10-day averages of time 
difference. Each day one new daily value is added 
and the 11-days old value is dropped and the re- 
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Line A=100 sec samples, frequency comparisons 1 sample/day at noon. 
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sultant 10 time differences averaged. The 30-day 
running averages of the above 10-day averaged 
values are obtained in a similar manner. The pre- 
cisions obtained are appreciably better, being in the 
order of +2 parts in 10° for 1 day and approximately 
+2 parts in 10" for 30 days. 


3.2. LF Standard Frequency Transmissions 


At present there are at least six low frequency 
standard transmissions being operated including 
MSF at 60 ke/s in Rugby, England; WWVB at 
60 ke/s in Boulder, Colo.; DCF77 at 77.5 ke/s in 
Mainflingen, Germany; OLP at 50 ke/s in Czecho- 
slovakia; HBJ at 96.04 ke/s in Switzerland; and 
A5XA at 1334 ke/s at Fort Monmouth, N.W., 
U.S.A. Other characteristics of these transmissions 
are described in Annex III of CCIR Report 66 
[1956]. Pierce [1957] has shown the standard devia- 
tion of frequency comparison for MSF transmissions 
on a transatlantic path to be in the order of 1 times 
10-” for a 30-minute observing period. Further 
calculations for LF transmissions are compared with 
experimental results by Watt and Plush [1959] and 
shown here in figure 8. Pierce [1958] has shown 
average frequency comparison precisions from 
WWYVB (60 ke/s) as received at Cruft Laboratory, 
of 1.4 parts in 10" in an observing period of about 
one hr. From equation (14) this would correspond 
to a short term time variation o (7) ~ 0.36 usec. 

Recent plans to distribute time signals via Loran-C 
navigation system, pulse transmissions at 100 ke/s, 
have been described by Doherty [1961]. Very precise 
time signal distribution in the order of 0.2 usec 
is possible within the ground wave coverage range 
of about 1,500 km. By employing sky wave modes, 
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| station or network of 2 or 3 stations whose primary 


timing precision in the order of 1 usec appears likely 
out to about 3,500 km. 

It is obvious from this and the preceding sections 
that in a limited area the LF systems (50 to 100 
ke/s) will provide for more precise comparison of 
frequency and time than the VLF systems (15 to 
20 ke/s). The propagation loss, however, is much 
higher in the 50 to 100 ke/s band which restricts the 
useful range to distances of about 3,000 to 4,000 
km. LF transmissions at 60 ke/s in England [Steele, 
1955] have shown that timing can be obtained to 
+15 psec by pulse envelope methods. 


3.3. VLF Standard Frequency Transmissions 


Annex III of CCIR Report 66 [1956] describes 
the transmissions from GBR at 16 ke/s which has 
been used as a basis for much of the pioneer work 
on standard frequency broadcasting in the very low 
frequency region. The precisions available from 
this transmission, and GBZ on 19.6 ke/s, based 
primarily on observations by Pierce [1957, 1958] and 
Allan, Crombie and Penton [1956], are included in 
sections 2.2 and 4. Stabilized transmissions with 
phase locked master-slave relationship from Hawaii 
and San Diego have also served as a valuable source 
of information for VLF path phase stability as 
described by Casselman, Heritage, and Tibbals 
[1959]. Recently interim station broadcasts of the 
U.S. frequency standard have been initiated from 
WWVL on 20 ke/s from near Boulder, Colo. Time 
and stabilized frequency broadcasts from NBA at 
18 ke/s also recently have been initiated.’ Addi- 
tional VLF transmissions are in the process of being 
stabilized in frequency which should make them 
available for path phase stability measurements. 
In the VLF region both theory and experiment 
show that greater observing times are required for 
a given degree of frequency comparison precision 
than is true in the LF region; however, when long 
range paths are considered the low attenuation rates 
obtaining in the VLF band make emissions in this 
frequency region useful for worldwide coverage. 


4. Characteristics and Expected Coverage 
for a VLF Standard Frequency and Time 
Signal Broadcasting Station 


Based on the previous analysis [Watt and Plush, 
1959] the precision of frequency comparison ex- 
pected as a function of observing period 7’ was 
prepared and is shown as figure 8 [Pierce, Winkler, 
and Corke, 1960]." It is obvious when comparing 
figures 7 and 8 that a VLF broadcast will provide 
one or two orders of magnitude improvement over 
present HF broadcasts in precision of frequency 
comparison for a given observing period. This 
great increase in performance has clearly indicated 
the desirability of designing and constructing a 
high quality VLF standard frequency broadcasting 


ES 
10 Precise time and stabilized frequency broadcasts at 18 kc/s over NBA, Sum- 
mit, Canal Zone, were announced in U.S. Naval Observatory Time Service, 
Notice No. 8, 18 November 1959. See also a recent paper by Stone, Markowitz, 
and Hall [1960]. 
11 A recent paper by Pierce, Winkler, and Corke [1960] describes results at 16 
ke/s over a 5,200 km path that indicate possible precisions of 2X10-", 
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function is to provide extremely precise standard 
frequency and time signals. 

Because of the desirability of having such a station 
near the U.S. standard of frequency located at 
Boulder, Colo., a study has been made to determine 
design parameters and service obtainable at various 
locations. Seventeen different receiving sites were 
chosen and the total number of measurements obtain- 
able each day for 99 percent of all hours were 
calculated based on a radiated power of 100 kw. 
Calculations were made for arange of radiated powers 
from 25 to 100 kw and the resulting number of 
measurements obtainable as a function of the cost 
per measurement was determined. The results of 
this study [Watt and Plush, 1959] indicated a 
preference for a radiated power in the order of 100 kw. 

The results of this study are summarized in table 1, 
where the various assumptions made are indicated 
along with a tabulation of the total number of 
measurements obtainable for 99 percent of all hours 
with 100 kw radiated at 20 ke/s ® for the months 
of June and December assuming a precision of fre- 
quency comparison of 1 part in 10°. The method 
of analysis and some of the assumptions made in 
this analysis are contained in the following sections. 


4.1. Propagation Path Attenuation Rates 

It is relatively well known that the attenuation 
rate (a) to be employed in field strength equations 
such as (2) in the paper by Watt and Plush [1959], 
expressed in decibels per 1,000 km (db/K), is depend- 
ent upon a number of factors including: frequency, 
ionospheric conditions, earth’s surface conditions, 
and the earth’s magnetic field. 

Experimental and theoretical studies indicate that 
there is a broad minimum in the attenuation-versus- 
frequency curve which is centered somewhere in the 
vicinity of 16 to 18 ke/s. The shape of this attenua- 
tion curve is likely to vary between day and night 
conditions, and is also likely to be different for 
east-to-west and west-to-east propagation. It would 
appear from an analysis of results obtained from 
many different sources that the east-to-west (mag- 
netic) attenuation rate in some cases may be as much 
as 1 or 2 db/K greater than the west-to-east attenua- 
tion rates. This effect is expected to be more 
pronounced at night; in addition this difference in 
attenuation rate with direction is expected to be 
greater at the lower frequencies; viz, 10 to 14 ke/s, 
and relatively small in the vicinity of 20 ke/s. These 
effects have been discussed by various authors [H. J. 
Round, T. E. Eckersley, K. Tremellen, and F. C. 
Lunnon, 1925; K. G. Budden, 1951 and 1952; J. R. 
Wait, 1958 and 1961b; and D. D. Crombie, 1958]. 

Approximate values for the average daytime at- 
tenuation rate (a) in the vicinity of 20 ke/s are: 
sea ~2, average land ~4, estimated values for perma- 
frost and icecap are 7 and 18. When mixed surface 
conditions are encountered along the path, as is 


2 It should be noted that 20 ke/s has been employed in these calculations 
because: (1) 20 ke/s is 9n internationally assigned standard frequency, (2) the 
power requirements for worldwide coverage are expected to be near a minimum 
at this frequency, (3) considerable propagation data is available in this region. 
Further research on attenuation and phase stability should be conducted over 
the whole VLF spectrum to determine the suitability of other frequencies for 
Standard frequency and time distribution. 





usually the case in practice, the total average path 
attenuation will not always be exactly equal to the 
sum of the attenuations expected from each individ- 
ual portion of the path. This may be due in some 
cases to a transfer of energy from the dominant to 
higher order modes caused by the discontinuities in 
conductivity at the surface. It is also likely that 
low conductivity at the earth’s surface will produce 
appreciably greater attenuation during the day when 
the ionosphere is low than at night when it is higher. 
After determining the various amounts of surface 
types along the individual paths, expected average 
attenuation rates have been assigned for each of the 
receiving locations in table 1. In some cases recipro- 
cal path values have been employed and, in general, 
values quoted are as high as anticipated from a 
consideration of all possible attenuation mechanisms. 
Field strengths anticipated for these particular 
paths have been calculated employing the relation- 
ships given by Watt and Plush [1959]; the results 
are shown in column 5 of table 1 for an assumed 
radiated power of 100 kw at a frequency of 20 ke/s. 


4.2. Expected Availability and Precision of Fre- 
quency Comparison 

The carrier to noise available at the receiving 
location may in some cases be the limiting factor as 
regards the time required to obtain a frequency 
comparison of, say, 1 part in 10°. Median atmos- 
pheric noise fields expected for the various receiving 
locations have been obtained from Crichlow, Smith, 
Morton, and Corliss [1955] and the results presented 
in columns 6 and 7 of table 1. Combining the pre- 
vious information regarding carrier field intensities 
with the noise field intensity data, we have obtained 
the minimum expected observation times required 
to compare frequency to 1 part in 10° for 100 kilo- 
watts radiated as shown in columns 14 and 15. 

Pierce [1957] has shown that the transmission path 
phase is very stable and that the diurnal variations 
caused by sunrise or sunset along the Rugby to Cam- 
bridge path are surprisingly constant throughout the 
vear. During the all sunlight period the path phase 
appears more stable than at night for frequencies of 
16 to 60 ke/s. We have determined, on the basis 
described by Watt and Plush [1959], the minimum 
observing times required to obtain a precision of 1 
part in 10° based on path phase variations expected 
for both day and night conditions. It should be 
noted that « (¢) has been assumed to increase as the 
square root of distance which is what one would 
expect of sections of variable delay wave guides 
connected in series. These results are presented in 
columns 12 and 13. 

Although it is obvious that in a case where mini- 
mum observing times based on path phase variabil- 
ity and carrier-to-noise ratio are nearly the same, 
the combined effect will require a longer observation 
time than indicated for either limitation. We have 
neglected this effect and have considered only the 
one condition which requires the longer period. 

In order to get some idea of the relative service 
obtainable from the proposed standard frequency 
and time signal broadcast, the total number of 
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measurements with a precision of 1 part in 10° 
obtainable per 24 hr for 99 percent of all hours have 
. been obtained by dividing the number of hours of 
all-day or all-night over the particular path by the 
appropriate minimum observing times and adding, 
as shown in columns 16 and 17 for the months of 
June and December. It can be seen from this tabula- 
tion that the service provided at various locations is 
appreciably different; however, i¢ is interesting to 
observe that in all but one case an appreciable 
number of measurements can be made throughout 
4 each day, and that in the worst case, for transmis- 
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where E is_the maximum or locked keyed carrier 
voltage, and dvjdt is the slope observed at the one- 
half carrier envelope amplitude point. 

When thermal type noise is present, the envelope 
of the keyed carrier can_ be expected to have a 
standard deviation in amplitude of o(a) for C/N>1 
with an envelope detector or for a C/N ratios if a 
| synchronous detector 1s employed. If the point at 
which the noise free keyed carrier crosses the //2 
point is chosen as a reference T,, it can be seen that 
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When the receiver effective bandwidth is B,(c/s) 
and we use the convention 0 C/N, xe, 1-€+ the rms 
carrier to rms noise 1n a one ke/s effective band, we 
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The design of ionosondes and the many electronic techniques used in the panoramic 


type are considered. 


A complete ionosonde which is readily transportable, economic to 


construct, and has the additional facility of electronic scan over any portion of the frequency 
range 1.0 to 20 megacycles per second is described briefly. 


1. Introduction 


It is becoming increasingly obvious that one of the 
most valuable ground based instruments for upper 
atmosphere research is the ionosonde recorder. 
Electronic computers and simplified reduction tech- 
niques such as those of Kelso [1954], or more recently 
Schmerling and Ventrice [1959], enable electron 
density distributions to be readily obtained from the 
corresponding h’f records, and recent efforts to 
tabulate these results on a worldwide basis are 
providing invaluable information on global iono- 
spheric morphology. The need for more frequent 
ionospheric samplings is suggested by the presence 
of traveling ionospheric disturbances described by 
Munro [1950] and Heisler [1958]. This work indi- 
cates that the ionosphere is a dynamic medium, in 
which almost continuous changes in ion density are 
occurring. Furthermore, recent discoveries of a 
traveling disturbance—F, relationship by Heisler 
and Whitehead [1960] and the sporadic nature of 
7, itself, suggest that spaced ionosonde recorders 
are essential to study growth of FE, patches and 
mechanics of ionospheric disturbances. 

Such recorders are, of necessity, very complex, 
and involve many diverse electronic techniques. 
The panoramic type recorder as originally described 
by Sulzer [1949] is now invariably used in modified 
forms such as that due to Heisler [1955], or as in 
existing American type C4 recorders [Brown, 1959]. 
However, the equipment is usually bulky and 
expensive. Also there seems to be lack of detailed 
design information available on this type of recorder. 

This paper considers ionosonde characteristics 
necessary to satisfy recent research requirements. 
Design principles of the panoramic type ionosonde 
are discussed in detail and finally a new recorder 
is described which is readily transportable and 
economic to construct, possessing the additional 
facility of electronic scan over any portion of the 
frequency range 1 to 20 Me/s. 


2. Basic Description 


The basic ionosonde recorder consists of a receiver 
tuned synchronously to a transmitter which scans 
over the frequency range fy to f; Me/s in the scan 
period 7’ sec, at the scan rate of S scans per hour. 
The transmitter is pulsed, the pulse length being 





t usec at a repetition frequency of F c/s, where F 
usually is equal to the local a-c supply frequency. 
The directly received pulse and echoes are displayed 
on one axis of a cathode ray oscillograph display 
calibrated in height, and a calibrated frequency 
scale is provided either by movement of the height 
scale across the axis of movement of the recording 
film or by movement of the film itself. 

In older types of recorders synchronization between 
transmitter and receiver was maintained by the use 
of cams and band switching. In the panoramic 
type recorder this is avoided by covering the fre- 
quency range in a single sweep. This is achieved 
by heterodyning an oscillator continuously variable 
from fof, Me/s to fit+f, Me/s with an f, Me/s 
pulsed oscillator and amplifying the resultant f, to 


ti Me/s signal to provide the transmitter output. 


The same variable oscillator is used as the local 
oscillator for a superheterodyne receiver of inter- 
mediate frequency f, Mc/s. This receiver has an 
untuned input stage coupled directly to aerial, and 
is thus tuned at any instant to the transmitter 
frequency insuring automatic tracking. 


3. Design Requirements 


The various factors which must be considered in 
arriving at particular design features are considered 
below: 

3.1. Frequency Range 


The upper frequency limit is normally imposed by 
the maximum value of F; critical frequency. At 
sunspot maximum this may exceed 15 Me/s but is 
seldom greater than 20 Mc/s. This limit may be 
exceeded by the critical frequency of sporadic E 
but this is exceptional and current research practice 
is to consider percentage of time fp, exceeds 5 Me/s. 

Several factors influence the selection of a lower 
frequency limit. At very low frequencies elaborate 
transmitter and receiver antenna systems are neces- 
sary due to radiation inefficiency, and at frequencies 
below 1.5 Me/s broadcast station interference in 
many areas causes considerable lack of resolution 
on the reproduced ionogram. 

Many ionospheric research projects require de- 
tailed examination of a small portion of the fre- 
quency spectrum. It is desirable therefore, that the 
frequency range be variable within the limits dis- 
cussed above. This also readily enables expansion 
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or contraction of frequency range to accommodate 
sunspot Maximum or sunspot Minimum conditions. 


3.2. Scan Period and Scan Rate 


There is now reason to believe that almost con- 
tinuous changes in ion density are occurring in the 
ionosphere, particularly in winter months when day- 
light records taken at intervals of 1 min often show 
marked differences from each other [Heisler, 1958]. 
In order, therefore, to prevent distortion and con- 
sequent misinterpretation of the record, particularly 
in deducing electron density profiles, scan period 
should be no greater than 15 sec. 

Consequent deterioration of integrated signal 
noise ratio as indicated by Dieminger [1959] does 
occur but should preferably be improved by means 
other than increasing scan period, such as careful 
receiver design or by photographic integration as 
discussed in section 4.4. 

The scan rate should be variable to allow several 
scans per minute in special studies e.g., sporadic-E 
and a low rate for economical normal use. 


3.3. Scan Law 


The increasing use of ionograms in electron 
density profile analysis demands higher scaling 
accuracies than previously, particularly in regions 
where ion density changes slowly with height. Irre- 
spective of accuracy of individual frequency calibra- 
tion points therefore, interpolation between these 
points must be accurate. King [1960] has stressed 
the value of a logarithmic law scan in certain simpli- 
fied N(h) reduction techniques, and this type of law 
is also of value in ionospheric prediction analysis. 

However, in practice the usual logarithmic law on 
many ionograms is only approximate and error can 
occur in interpolating between frequency marks. 
Interpolation is expedited by use of a linear scan law 
and small variations in this will usually introduce 
errors of less magnitude than those due to difficul- 
ties in interpolating an approximate logarithmic 
scale. One feature of log scan law is the open fre- 
quency scale at low frequencies where considerable 
accuracy is often required. However, this is usually 
at the expense of accuracy in reading fF». 


3.4. Transmitter Pulse Length and Receiver 
Bandwidth 


Since minimum detectable signal power is inversely 
proportional to pulse length [Herb and Sinsheimer, 
1947], pulse lengths should be as long as possible 
consistent with required resolution. Price [1960] 
suggests that there is little advantage in pulse length 
t being shorter than 100 uS due to loss of resolution 
by ionospheric dispersion. This argument, however, 
does not apply to sporadic £ ionization where a 
pulse length of 35 wS would be required to resolve 
the 5 km stratification observed by Briggs [1958]. 

Pulse length is also a factor in determining receiver 
bandwidth and the optimum value is 1.2/t Me/s 
[Purcell, 1947]. 








4. Design Principles 


The following design principles apply mainly to 
panoramic ionosondes. For purpose of detailed 
discussion the recorder may be considered to con- 
sist of the following main sections: variable frequency 
oscillator, fixed frequency oscillator, transmitter, 
receiver, and calibration system. 


4.1. Variable Frequency Oscillator (VFO) 


This oscillator provides a substantially constant 
output over the frequency range (fo+f,) Me/s to 
(f:+-f-) Me/s in the sean period 7'sec. Since the wide 
frequency spectrum output supplies energy to three 
individual circuits, namely transmitter mixer, re- 
ceiver mixer, and frequency calibrator it is essential 
that oscillator output impedance be low. High 
voltage output is not necessary as receiver and 
calibrator require only a few volts for satisfactory 
operation, and for reasons discussed later in section 
4.3a high level mixing in the transmitter mixer stage 
is difficult to achieve. 

The rate of change of frequency (f-+-f,) Mc/s of the 
oscillator with time ¢ determines the scan law of the 
recorder, and this law is usually logarithmic or linear 
in form. If a rotating capacitor is used to vary the 
oscillator frequency it is possible to shape the rotor 
plate so as to give the desired rate of charge of 
capacitance with rotation angle. In_ practice, 
however, only an approximation to the required law 
is obtainable due to unpredictable edge effects be- 
tween rotor and stator, and an exact law necessitates 
use of a rotating cam profile with a cam follower 
mounted on the capacitor shaft. 

Any suitable rotating capacitor may be made to 
give a similar approximation to the required law by 
applying three point tracking techniques [Green, 
1943] by addition of a series capacitor. 

Electronic control of scan frequency may be 
obtained over a limited range by use of the recently 
developed voltage variable capacitor or Varicap 
[Hammerslag, 1959] or by application of a ring type 
resistance capacity feedback oscillator developed by 
Cormack [1951]. The frequency of this oscillator 
can be varied over a considerable range by applica- 
tion of a control voltage. 

Electronic sweep control has the additional facility 
that any sweep range f) to f; can be obtained by 
suitable preset biasing arrangements while still 
maintaining scan law. <A similar facility in 
mechanical scan arrangements can only be obtained 
by changing cam profiles. This form of control also 
removes the usual limitations on scan rates placed by 
mechanical sweep arrangements. Scans may be 
repeated as often as desirable, and minimum scan 
periods are limited mainly by possible receiver 
detuning before return of ionospheric echoes. 


4.2. Fixed Frequency Oscillator (FFO) 
This oscillates at a frequency of f, Mc/s and since 
the highly sensitive receiver in the recorder has a 
similar intermediate frequency, it is essential that 
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the oscillator be pulse modulated by the transmitter 
keying pulse to prevent continuous receiver paralysis. 
It is extremely difficult even with most elaborate 
oscillator shielding to prevent leakage of signal into 
the receiver circuits. 

The inherent accuracy of frequency calibration 
depends upon the accuracy of this oscillator. How- 
ever the necessary requirements are not stringent 
and it is shown in section 5.5 that there is little 
advantage in having a fixed frequency oscillator 
accuracy greater than about one part in three 
thousand. This can be readily obtained with 
sufficient long term stability by means of a carefully 
constructed pulsed Hartley oscillator [Gamertsfelder 
and Holdam, 1949]. Pulse modulation of crystal 
oscillators is difficult and it is usually necessary to 
crystal control an oscillator at a submultiple fre- 
quency f,/n Me/s and multiply this by pulsed ampli- 
fiers to give the required f, Mc/s output. Unfortu- 
nately there may be leakage of continuous f, Me/s 
harmonic radiation into receiver circuits. 


4.3. Transmitter 


This consists essentially of a mixer to obtain the 
sweep frequency from the variable frequency oscil- 
lator output and the pulsed output of the fixed 
frequency oscillator, followed by a wide band ampli- 
fier feeding the transmitter aerial. These require 
a balanced output and hence push-pull output is 
usually employed, with the additional advantage of 
providing higher power output particularly in class 
B operation. Adoption of this necessitates some 
form of phase inversion in the amplifier chain to 
provide out of phase driving voltages for the output 
stage. In order to obtain sufficient gain with a 
minimum of tubes, wide band coupling techniques 
must be used throughout mixer and amplifier cir- 
cuits. The various approaches which can be made 
to these problems are considered below. 


a. Resistance Capacitance Coupling 


In this type of coupling the top frequency to be 
amplified /,; is approximately equal to bandwidth 
B,. If the amplifier consists of n stages and if B is 
the bandwidth of each stage 


B 
Vn 


The normal maximum plate resistor F in any stage 
is given by 


B,= (1) 


ed 
2 BC 
where C is the interstage shunt capacity. 
R may be increased by a factor of 1.77 to 2.4 by 
the use of two and four terminal network coupling 
techniques, described by Walker and Wallman 
[1948]. The total interstage capacity C should be 
measured as strays are unknown and _ published 
values of tube capacitances are cold values and have 
been found to deviate considerably from hot tube 
values. 
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Some resistors have been found to change value 
considerably at high frequencies, and /? for prefer- 
ence should be deposited carbon on a ceramic former. 

It is possible using network coupling techniques to 
deliberately peak the high frequency gain of an 
individual stage [Walker and Wallman, 1948]. If 
this is done then bandwidth shrinkage as expressed 
by eq (i) is not so severe, and bandwidth. require- 
ments for an individual stage may be relaxed. 

The limitations imposed by maximum power 
rating of tubes may be overcome by pulse switching 
stages for the duration of the transmitter pulse. 
Screen voltages may then be increased by a factor 
N with consequent increase in mutual conductance 
9m to VNGm approximately. Moreover the effective 
grid base of a tube is lengthened by a factor N 
approximately so that higher permissible grid volt- 
age swings become possible, an important feature in 
high level stages. 

One immediate consequence of pulse switching is 
that a large negative pulse appears across plate 
resistor R due to change in plate potential level 
when the valve becomes operative, and small 
coupling time constants consistent with adequate 
coupling at frequency fo are necessary to prevent 
appearance of this on the grid of the following stage. 

The valve chosen as transmitter mixer must have 
a high value of g, since conversion transconductance 
is approximately g,,/3. This restricts choice of 
mixer valves to those with a relatively short grid 
base. Consequently high voltage inputs are not 
possible and high level mixing is difficult to achieve. 


b. Wide Band Transformer Coupling 


One difficult problem associated with the trans- 
mitter design is to provide sufficient power output 
into a low impedance. In practice the push-pull 
output stage is directly coupled to the aerial circuits 
to develop power in an effective load of 600 ohms 
under class A operating conditions. The type of 
tube chosen as final amplifier must therefore have 
high mutual conductance, low plate resistance and 
large negative grid voltage for plate current cutoff, 
conditions which are most readily met by pulse 
modulator tetrodes specially developed for radar 
purposes. 

Performance can be improved by matching aerial 
impedence to the plate resistance of the final ampli- 
fiers by means of a wide band transformer. This 
enables class B operation to be used with its much 
higher plate efficiency and provided sufficient power 
is available from the previous stage, grids may be 
driven positive, consequently valves need not have 
a high value of grid voltage cutoff. This relaxation 
in required characteristics permits use of normal 
high perveance valves and effects considerable 
economy. 

A complete discussion of very wide band trans- 
former design is given by Maurice and Minns [1947]. 
The upper frequency limit is mainly controlled by 
leakage inductance which is minimized by keeping 
number of turns low and careful attention to physical 
design. At low frequencies number of turns must 





be sufficient to provide required inductance, and is 
reduced by using a core of modern powdered iron 
dust materials possessing high initial permeability. 
However, care must be taken to insure that flux 
saturation does not occur at the power outputs 
considered. If the transformer has balanced wind- 
ings, then a high degree of balance can only be 
obtained by inserting an electrostatic shield between 
windings. 
c. Phase Inversion 

If the transmitter has push-pull output, some 
form of phase inversion is necessary in the amplifier 
chain to provide driving voltages 180° out of phase. 
If the phase splitter is used in the early stages of 
the amplifier then all subsequent stages must be 
‘push pull in form, with consequent duplication of 
valves, increase in power consumption and the risk 
of voltage unbalance occurring at the output stage 
grids. When the phase splitter is placed in late 
stages of the amplifier, it may be necessary to 
employ large pulse tetrodes. 

A convenient means of obtaining the necessary 
inversion is by using an amplifier having unity gain 
as described by Heisler [1955], or the cathode 
coupled unbalance to balance system used in the 
C4 recorder [Brown, 1959]. 

Another method of realizing the necessary anti- 
phase voltages is by use of a broad band transformer 
with a balanced secondary winding. This can be 
used to advantage as a means of coupling the driver 
stage to power output stages operated under class 


B conditions. 
d. Transmitter Pulser 


™ This applies a large pulse of suitable polarity to 
either grid or screen circuits rendering the valve 
operative for a period of ¢t wsec. During this period 
the pulser must present a low impedance to the 
appropriate circuits of the controlled valves. This 
permits loading by complex impedences in pulsed 
grid circuits and enables pulsed screen circuits to 
be operated without RF bypass capacitors, both 
factors impairing pulse shape, and inhibiting feed- 
back between common pulsed stages. 

In practice the bootstrap type pulser [Glasoe, 
1948] pulsing screen grids has been found most 
satisfactory 

e. Transmitter Aerial 


Effective power output of the transmitter depends 
to a large extent on the radiation efficiency of the 
associated antenna system. The broadband re- 
quirements demand design compromise particu- 
larly at the low frequency end, where physical 
dimensions of the antenna have to be at least a 
quarter of a wavelength for efficient radiation. 
Usually a 600-ohm delta type antenna is employed 
similar to that described by Cones et al. [1950]. 
Unfortunately, considerable power is dissipated in 
the terminating resistance and at low frequencies 
as much as 80 percent of the power input to the 
aerial may be lost in this manner [Cones, 1951]. 


Where research is concerned with this end of the 
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spectrum, it may be necessary to provide high power 
input to the aerial system to compensate for this. 

There appears to be considerable need for de- 
velopment of suitable aerial designs for this purpose, 
particularly with a view to maintaining aerial 
efficiency at low frequencies. 

4.4. Receiver 

This consists of an untuned aerial stage and mixer 
followed by sufficient intermediate frequency ampli- 
fication, to provide several volts of signal at the 
detector with high signal noise ratio. Since all 
frequencies arriving at the aerial appear on the 
mixer grid, it is necessary to attenuate frequencies 
outside the signal band fy to f; Me/s to prevent 
serious cross modulation. In practice a high pass 
filter is inserted in aerial circuits to remove broadcast 
station interference below f, while a resonant trap 
tuned to f, Me/s prevents breakthrough at the 
intermediate frequency of the receiver. 

A similar restriction of pass band may be obtained 
by using a well designed wide band transformer 
on aerial input circuits. Aerial input is usually 
balanced in form and use of a balance to unbalance 
transformer greatly simplifies front end design. 

Bandwidth requirements are mentioned in section 
3.4. Intermediate frequency f, is usually greater 
than 20 Me/s, hence gain and selectivity can be 
improved by using double frequency conversion. 
As will be described later the frequency calibration 
system uses a 1 Me/s harmonic generator; the value 
of second intermediate frequency should therefore 
be other than integral to prevent interference from 
this source. 

The amplitude of the received pulses covers a large 
dynamic range and if receiver gain is adjusted to 
record weak echoes satisfactorily, strong echoes 
invariably saturate the receiver with consequent 
loss of complex splitting detail. 

This effect may be minimized by some form of 
instantaneous automatic gain control [Lawson and 
Uhlenbeck, 1948a], or alternatively by some form of 
logarithmic receiver such as that described by 
Chambers and Page [1954], which has a logarithmic 
gain response over a very wide range of amplitude 
of input signals. The latter in particular is most 
suitable for measurement of absolute amplitudes 
in ionospheric absorption research applications. 

Automatic gain control is particularly effective in 
discriminating against unwanted CW signals. Sim- 
lar discrimination by differentiation of the received 
echo before application to display should be used 
with caution or loss of complex splitting detail in 
the received echo may occur. The usual RC filter 
is not satisfactory [Lawson and Uhlenbeck, 1950b] 
and a high pass filter with low frequency cutoff 1/(2¢) 
should be used. The technique is only of real value 
if intermediate frequency amplifiers are not over- 
loaded [Haworth and Jordon, 1947]. 

Receiver output should provide echoes of negative 
polarity for application to the grid of the cathode 
ray display tube. Photographic integration on the 
resultant ionogram has been shown to provide con- 
siderable noise discrimination. 
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4.5. Calibration System 


The final ionosonde presentation must be provided 
with calibration of both frequency and height. Fre- 
quency calibration is usually provided by hetero- 
dyning the variable frequency oscillator output 
fy tf, to fi +f, Me/s with harmonics from a 1 Me/s 
harmonic generator and using the resultant zero 
beats to generate short sharp pulses corresponding 
to each megacycle in the range fy to f;. These pulses 
are then applied to intensifier circuits in such a way 
to provide suitable marks on the visual display. 

The generator is usually crystal controlled and 
accuracy of individual marks depends mainly on 
accuracy of f,. However, since the received echo 
falling on a given frequency mark includes a number 
of frequencies in the optimum receiver bandwidth, 
it is doubtful whether frequencies read to any greater 
accuracy than one half of this have any real signifi- 
eance. This demands a FFO stability of 1. part 
in 3,000. 

Height calibration appears on the final display as 
a series of sharp marks corresponding usually to 
50 km height separation. These are generated by 
a 3,000 evele oscillator gated so that it only operates 
during each height display. 

An upper limit to overall accuracy of height cali- 
bration of at least 1 km is imposed by errors in time 
measurement due to rise time of received echoes 
and variations of receiver delay with signal ampli- 
tude [Chance, 1949]. The maximum height range 
usually employed is 1,000 km, and therefore accuracy 
is limited to one part in 1,000. This order of 
stability is easily attained in a well designed RC 
oscillator. 


5. Recorder Description 


An ionosonde is described which covers the fre- 
quency range 1 to 20 Me/s. This has been con- 
structed in the laboratory of the Radio Research 
Board Sydney, based on many of the design tech- 
niques discussed above. The complete recorder weighs 
approximately 150 Ibs and is contained in a 22 in. 
square cubicle which stands 34 in. high (without 
camera and camera boxes). The ionosonde features 
complete electronic operation, with a nominal scan 
period of 15 sec repeated at a scan rate of 1 per 
minute, with a linear or logarithmic sean law. 
However, scan period can be made as short as 2 sec 
and this can be repeated almost continuously. 
Nominal frequency range is 1 to 20 Me/s but any 
part of this range may be selected by the setting of 
two controls. A manual control permits a manual 
scan to be made, or if necessary the ionosonde may 
be set to operate on any particular fixed frequency. 
Transmitter pulse length is variable between 40 and 
100 ysee and the pulse peak power input to a 600 
ohm balanced antenna averages 3 kw at a pulse 
repetition frequency of 50 c/s. Most of the power 
supplies are electronically voltage regulated. The 
complete recorder can be constructed for less than 
£800 (Australian). 
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Because of space considerations it is obviously 
impossible to publish complete circuit diagrams or 
to give detailed technical information. The following 
description therefore lists briefly the most important 
technical features of the equipment. 


5.1. Variable Frequency Oscillator (VFO) 


This is developed from a circuit due to Cormack 
[1951]. Basically the circuit is a resistance capaci- 
tance ring type oscillator, in which the necessary 
phase shift is provided by a triode amplifier and 
four cathode followers. Phase rotation and hence 
oscillator frequency is controlled by varying the 
mutual conductance of the cathode followers with 
an applied voltage. The voltage developed by 
the oscillator is from 1.5 to 1.0 volts rms over the 
frequency range 29 to 48 Me/s and the almost 
linear frequency range shift for applied control 
obtained is amplified to supply 4.0 to 6.0 volts 
rms to the transmitter and receiver mixers. The 
transmitter mixer is supplied directly from this 
amplifier while the receiver mixer and frequency 
calibrator mixer are coupled to the amplifier by 
cathode followers. 


5.2. Fixed Frequency Oscillator (FFO) 


This is a pulsed Hartley oscillator as described 
by Gamertsfelder and Holdam [1949a] in their 
figure 4.45. The oscillator tube operates over a 
linear region of the tube characteristics and is 
therefore free from frequency deviation effects due 
to stray power amplifier pulses in the grid circuits. 
C class pulsed oscillators are particularly prone to 
this. A buffer amplifier is used after the oscillator 
to drive the transmitter mixer cathode. 


5.3. Transmitter 


A type 12BY7 tube was chosen for the transmitter 
mixer from several types tried. The VFO voltage 
is injected into the control grid and the FFO voltage 
into the cathode circuit. The plate circuit is 
coupled to a 6CK6 amplifier by a four terminal 
network with a cutoff frequency of 22 Me/s, which 
inhibits higher order component frequencies from 
reaching the 6CK6 grid. This amplifier is similarly 
coupled to a type 6CM5 amplifier which is trans- 
former coupled to the power amplifiers. The 
transformer is a wide bandwidth type wound on 
a toroid A4 Ferroxcube core, providing out of phase 
voltages to drive the final push-pull amplifiers. 
These are type 6CM5 tubes coupled to antenna 
circuits by another wide bandwidth transformer. 
All amplifiers are screen grid pulsed by a common 
bootstrap type pulser unit. 


5.4. Receiver 


Type 6BX6 tubes are used as RF amplifiers and 
these are coupled to type 6EJ7 tubes as mixers by 
four terminal coupling networks. The variable 
oscillator voltage of average value 1.8 v rms is 








injected into the cathode circuits of the mixer tubes. 
The mixer stage is followed by a single stage of 
IF amplification at 28 Me/s which is then converted 
to 1.9 Me/s and amplified by two additional stages. 
Dynamic range is improved by use of variable MU 
tubes throughout the amplifier and application of 
automatic gain control to the last stage. This is 
followed by a _ conventional compensated diode 
detector and one stage of video and audio amplifica- 
tion. The receiver has an overall gain of 114 db 
and a bandwidth of 30 ke/s. 


5.5. Calibration System 


This provides suitable range marks and frequency 
marks on the final display. The height calibrator 
uses a gated resistance capacity 15 ke/s oscillator 
and this signal is divided to provide a 3 ke/s signal. 
The 15 ke/s and 3 ke/s signals are clipped, suitably 
differentiated, mixed and applied to the display to 
produce fine 10 km and broad 50 km height marks. 
Either or both calibrations may be selected at will. 

The frequency calibrator consists of a 6BLS8 tube 
as a triode 1 Me/s erystal oscillator and pentode 
harmonic generator. This output is mixed with 
the output of the variable frequency oscillator in 
the pentode section of a 6AJ8 tube and the resulting 


heterodyne beats which correspond to 1 Me/s fre- | 


quency intervals are amplified, squared and applied 
to a multivibrator which when operative, removes 
the time base synchronizing signal and hence the 
cathode ray tube display. Special provision is made 
to insure that frequency intervals always occupy the 
same deflection length on the display irrespective of 
duration of scan period, and a later improvement 
places frequency markers on top and bottom of the 
record only. 


“hae. 





Ficure 1. 
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5.6. Control Unit 


The control unit provides the basic variable d-c 
voltage to vary the frequency of the VFO. The same 
voltage is used to provide display tube deflection 
when a linear scan law is desired. An adjustable 
RC charging circuit with suitable controls to vary 
amplitude and initial starting voltage establishes 
start and finish of scan voltage and hence the par- 
ticular frequency scan. A synchronous motor is 
used to provide regular 1 min or longer intervals 
between scans. 


5.7. Display 


Two displays are provided, one is on the face of a 
5 in. long persistence screen cathode ray tube and 
is used for monitoring, while the second display is 
on a similar tube having a normal persistence screen 
which is photographed with either a 35 mm or 16 mm 
camera. 

Typical positive prints of ionograms produced by 
the recorder are shown in figures 1,2, and 3. Figures 
1 and 2 are both examples of 12 sec sweeps. In 
figure 1 the frequency range covered is from 1 to 14 
Me/s, while in figure 2 the frequency range is 2 to 
20 Me/s. Figure 3 shows a very fast sweep of 2 sec 
over the frequency range 0.1 to 12 Me/s. It will 
be noticed that the sweep is so rapid that individual 
strokes of the time base are obvious, yet in spite of 
this the amount of information available does not 
differ materially from that in figures 1 and 2. The 
examples serve to illustrate the great versatility of 
the equipment both in seen period and frequency 


range. 





Typical ionogram covering the frequency range 1.0 to 14 Mc/s in a scan period of 12 sec. 





be 
fea 
wh 
fac 
pre 
des 











Figure 2. Typical ionogram covering the frequency range 2.0 to 20 Mc/s in a sean period of 12 see. 





Figure 3. Typical tonogram covering the frequer 
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6. Conclusions 


Many facets of panoramic ionosonde design have 
been considered. There are possibly other design 
features, particularly for specific research projects 
which have not been included, but it is felt that the 
factors discussed constitute enough information to 
provide a basic ionosonde design. Some of the main 
design features emphasized are as follows: 

(i) The use of electronic scan control provides 
many desirable ionosonde features eco- 
nomically. Rapid scans over the whole or 
any part of the frequency range are pos- 
sible, and sean laws either logarithmic o1 
linear can be selected at will. 

Gi) Considerable economy can be effected in 
transmitter circuits by the use of wide 
band transformer coupling. 
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wcy range 1.0 to 12 Mc/s in a scan period of 2 sec. 


(ili) There is little advantage in using trans- 
mitter pulse lengths less than 100 usec, 
except for FE, investigations where a pulse 
length of 40 usec may be beneficial. 

Present transmitter antenna systems are 
very inefficient at low frequencies. There 
is need for research into suitable aerials 
for this particular purpose. 

(v) To avoid technical design difficulties, the 
frequency range of the ionosonde should 
not be greater than necessary. A range 
of 1.0 to 20 Me/s is suggested. 

The sean period should not exceed 15 sec. 

The sean rate should be frequent at least one 
frame per minute, although this of neces- 
sity is governed by economy in film con- 
sumption. 

(viii) Differentiation of receiver output is ad- 

vantageous, only if intermediate fre- 
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quency stages are not overloaded and 
care is taken in requisite filter design. 

(ix) Considerable noise discrimination is obtained 
if the signal applied to display cathode 
ray tube grid is of negative polarity. 

(x) There is some doubt as to the advantage of a 
logarithmic scan law, and interpolation 
errors in such a case are likely to be 
greater than those occurring in a similar 
linear law. Electronic sean 
particularly adaptable to generation of a 
linear sean law. 

(xi) For most purposes a receiver bandwidth of 
30 ke/s is more than adecuate even for 
sean periods of the order of a few seconds. 

A complete technical description will be given in a 
later paper. 
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A method, based on the maximum median hourly value of fF, month-by-month as 
observed at Washington, D.C., is given for predicting a smoothed annual sunspot number 


immediately at the close of a given month, centered on the month in question. 
It should be a useful suppiement to the MeNish- 


equations and standard errors are given. 


Regression 


Lincoln method of prediction of sunspot numbers, particularly during the first two years 


of the rising part of a sunspot cycle. 


It is capable of use with the observations from any 


ionosphere sounding station operated in a consistent Manner over a period of years, prefer- 


ably during at least two solar cycles. 
1. Introduction 


One of the widely accepted methods used in the 
prediction of sunspot numbers for use in radio 
frequency predictions is the MeNish-Lincoln method 
[1949]. This method involves the use of smoothed 
annual sunspot numbers corrected by departures one 
vear earlier than the time in question. It predicts 
month by month the ‘‘best”’ prediction for 1 year 
ahead. 

The use of a 12-month running average suffers 
of course from the fact that to obtain the sunspot 
number applicable to a certain month, 6 months in 
the future must be observed. When, however, the 
authors extended the technique to the prediction of 
3-month mean sunspot numbers, lower overall 
reliability resulted. The method outlined here has 
the advantage that, immediately upon the comple- 
tion of a month, a prediction of the 12-month running 
average sunspot number centered on that month can 
be made. Of considerable importance is the fact 
that the rms deviation of prediction would seem to 
be of about the same order of magnitude as those 
given in a second paper [McNish and Lincoln, 1954] 
in which the authors examined the reliability of 
their method of prediction and applied it to cycle 18. 

In the proposed method, the highest monthly 
median nh of fof, for any hour, as observed at 
Washington, D.C., is used, in a set of 12 regression 
equations, one for each month, to obtain an estimate 
of the 12-month running average sunspot number 
for the month that has just ended. It is believed 
that this choice of maximum fo/2, rather than fof» 
for a fixed time of day, is of aid in minimizing 
seasonal effects. 

The idea of using fo¥2 as a measure of solar activity 
is not new, but the methods so far proposed have 
been more elaborate, involving for instance such 
things as 5-hourly averages centered on noon, and 
averaging of results from several stations [Phillips, 
1947; Minnis, 1955]. Also much work has been 








done and is continuing on other possible indices 
[Kundu and Denisse, 1958; Minnis and Bazzard, 
1959]. The fact that a reasonable answer can be 
obtained to the problem of predicting solar activity 
by the simple method outlined here, gives reason 
for hope that the much more elaborate methods 
presently under study will eventually result in a highly 
satisfactory index obtained from some combina- 
tion of physical quantities. At the same time, it is 
distinctly worthwhile to have some way of connect- 
ing future ionospheric indices with the long time- 
series of sunspot numbers now available, on which 
so many relationships have previously been based. 

The method of this paper is intended to show 
what can be done almost immediately using data 
from but one ionosphere station (Washington) which 
has been operated in a consistent fashion over a 
long period of time. The initial impetus for this 
study was the sudden drop in solar activity in 
January 1961 as represented by the provisional 
Ziirich monthly sunspot number of 53.5, together 
with the desire to investigate what this sudden 
drop might mean in terms of future predictions of 
critical frequencies. 


2. Correlation Coefficients 


Various results are given in table 1. It will be 
noted that, in general, the coefficients in column 
A exceed those in column B, and those in column 
E exceed those in column D. (However, these 
observed tendencies by themselves are not statisti- 
‘ally significant.) Column C is filled out only 
for January and May, the months for which the 
highest and lowest coefficients were obtained in 
column A. The omission of numbers for the maxi- 
mum period 1957 to 1959 (col. C) and the re- 
duction of high numbers in this period (col. E) 
gave correlations which, when taken with those in 
columns A and D, form an interesting series of 
coefficients for these months. 
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TaBLE 1. Correlation coefficients, maximum hourly median 
foF2 (Washington) versus Ztirich sunspot number for various 
periods; A, C, D, versus 12-month running average sunspot | 
number centered on the month in question; B, versus monthly | 
sunspot number; E, same as D, but with sunspot numbers | 
> 150 used as 150 


A B Cc D E 


1945-1960 1945-1960 | 1941-1956, 1960 1941-1960 1941-1960 
| 
Jan 0. 989 0. 97 0. 985 0. 989 0.975 | 
Feb : . 98 . 94 . 98 98 | 
Mar 97 97 96 | .97 
\ 93 | - 92 | 94 97 =| 
May 88 | 91 | 91 | 90 93 | 
June. . 92 . 89 93 % | 
July. .92 .93 94 96 | 
Aug . 92 . 89 . 93 95 | 
Sep. 95 95 96 98 | 
Oct. . 96 . 93 .97 a 
Nov...- 98 94 98 8 | 
Dec .97 .94 . 98 . 98 


The use of a value of 150 for all running average | 
sunspot numbers over 150, produced an equal or 
better correlation, except for January. This is in 
agreement with previously noted results [Bennington, 
1958; Ostrow, 1959] and with CRPL predictions 
practice during the recent maximum.' In January, 
at the recent maximum period, increased solar 
activity did seem to have a marked effect in raising 
the criticals at Washington to a very high value 
in the middle of the day, higher than C RPL pre- 
dictions based on sunspot number 150. Also. it 
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Figure 1. Scatter of points and regression line for March 
(tables 1 and 2, col. E). 


will be noted that correlations are somewhat higher 
for the winter months, a not unexpected result, as it 
has long been noted at Washington that the range 
of variation in f,/; over the solar cycle, is much 
greater for winter months than for summer months. 
While it might be argued that the apparent seasonal 
effect was due to the way the sunspot numbers 
happened bes fall, e.g., if there had been four num- 
bers for January around 200 and none for June 
greater than 150, it should be remarked that during 
the recent maximum (1956 to 1959) each month 
had one 12-month running average sunspot number 
> 187, and each month had three >150. Column 
B is included to show that correlations can be very 
high between monthly sunspot numbers and critical 
frequency. Figure 1 illustrates the relationship 
between f and FR leading to the correlation coef- 
ficient 0.97 in column E, table 1 , for March. 


3. Regression Equations 


Table 2 gives the linear regression equations for 
each month of the year, based on the coefficients in 
column E of table 1. In addition, equations are 
given for October through February using the 
coefficients of column D, which are greater than or 
equal to the corresponding coefficients in column E. 


TABLE 2. FPegression equations; R= 12-month running average 
sunspot number, f=maximum hourly median foF (Mc/s); 
columns D, E correspond to D, E in table 1 





D E 
Jan P=83.2+20.5(f-9.6) R=78. 4+17.8( f-9.6) 
Feb R=83.4+21.7(-9.6) 
Mar 
Apr 
May =79. 24+-39.4( f-7.3 
June pave. 5+56.0(/-6.8) 
July R=79.8+53.4(/-6.4) 
Aug R=80.0+44.1(f-6.8) 
Sep Steed P=79.9+26.3(f-8.1) 
Oct R=84.2+22.9( f-10.0) PR=79.8+20.3( f-10.0) 
Nov.... RP =84.2+21.4( f-10.3) PR =79.6+18.8( f-10.3) 
Dec P=84.0+22.9( f-10.0) P=79.4+20.2( f-10.0) 


4. Standard Error of an Estimated Sunspot 
Number 


Table 3 gives the low and high values of f (maxi- 
mum hourly median fo>¥2 in Me/s) month by month 
over the 20-year period 1941 to 1960, as well as the 
means, which, of course, are repeated from table 
2. Table 4 gives the standard error of an estimated 
sunspot number, /?, for each value of f given in table 
3. These are computed from the formula 


SR est.) = Sarl iad 1 Dr [Snedecor, 1946] 


where x is the number of observed pairs (f;, 2,), 


f= Dosiln, sez is the sample standard deviation about 


the regression line, and f is the predicting value of 


fof, as before. They may be compared with the 
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TaBLeE 3. Range of values of f (Mc/s), Washington, 1941-1960 
Jan. Feb. Mar. Apr. May June July Aug. Sep Oct Nov. Dee. 
Low 5.9 5.7 5.4 5.3 5.2 5.3 4.8 5.0 5.2 6.2 6.5 6.2 
Mean : 9.6 9.6 9.2 8.2 7.3 6.8 6.4 6.8 8.1 10.0 10.3 10.0 
High a 15.5 14.0 13.2 11.2 10.0 8.0 reer f 8.6 11.0 14.2 14.8 14.0 
TaBLE 4. Standard error of an estimated R for various values of f (see table 3) 
Jan. Jan. Feb. Feb. | Col. E, table 2 Oct. Oct. Nov. Nov. Dec. Dee. 
Col. D, | Col. E, | Col. D, | Col. E, |__ etal _ Col. D, | Col. E, | Col. D, | Col. E, | Col. D, | Col. E, 
table 2 | table 2 | table 2 | table 2 table 2 table 2 table 2 table 2 table 2 table 2 
Mar.) Apr. | May | June | July | Aug. | Sep. 
For low f 10.9 13.3 13.7 11.1 | 14.3 | 16.2 | 23.5 | 18.1 | 17.8 | 19.3 | 13.0 17.8 16.1 12.8 12.4 14.5 12.2 
For mean /_. 10.6 12.8 13.1 10.6 | 13.8 | 15.5 | 22.1 17.1 16.8 | 18.3 | 12.4 17.0 15.4 12.3 11.9 13.8 11.6 
For high f_-- ef 14.1 13.9 11.3 | 14.4 | 16.2 | 23.8 | 17.7 | 17.5 | 19.3 | 13.0 18.0 16.3 13.0 12.6 14.6 12.3 


rms deviation of prediction given in table 1 of 
MeNish and Lincoln [1954], in which the deviations 
range from 7.0 to 19.2. Although these are of the 
same order of magnitude, it must be remembered 
that the MeNish-Lincoln formulas were obtained 
using data from 10 solar eyeles, whereas the standard 
errors of this paper were calculated using less than 
2 cycles. Essentially the sunspot numbers used in 
the MeNish-Lincoln method form a series of numbers 
spaced one solar cycle or about 11 years apart, 
whereas in this method the sunspot numbers are 
spaced 1 vear apart. 


5. Results and Conclusions 


Since this paper was begun, values of maximum 
hourly median fo/2 at Washington have been received 
for January through May 1961. These are listed in 
table 5, together with the values of 12-month running 
average Zurich sunspot number estimated by the 
method of this paper, and those predicted by 
the MeNish-Lincoln method. It must be remem- 
bered, of course, that the latter predictions were 
made, for the corresponding months, 6 months 
earlier than the predictions using the critical fre- 
quencies at Washington. 

At the time of final version of this paper (June 
5, 1961) the provisional 12-month running average 
Ziirich sunspot number, centered on December 1, 
1960, had dropped to 84.3. At the average rate of 
decline for the past 4 months, the number centered 
on the middle of January 1961 would be 76.7, 
compared with the predicted values 75.0 and 97.4 
of table 5. 

TABLE 6. Corr-lation coefficients, regression equations, 


Corr. coeff. 


Station Geomag. 


lat. 
Jan, July 
Fairbanks_- -- _ = : 64° N 0.98 0.91 R 
Washington ___- 50° N 99 92; R 
Christehurch- -- 48° S O4 2.96 R 
Ituancayo- ----- 1° § 97 .96 | R 


«15 years; data missing for July 1955. 
» Standard error of an estimated R for mean /. 


Values of f at Washington, estimated values of R, 
R by McNish-Lincoln method, Jan.- 


TABLE 5. 
predicted values of 
May 1961 


R (Est.) 


| 1961 f (Me/s) R (Pred., MeNish-Lincoln) 
| 

Jar 9.2 97.4 

Feb 8.6 93.8 

Mar 8.6 90. 2 

Apr 7.4 88.5 

May 6.8 86. 2 





« From table 2, col. D; all others from table ‘ 


A quick method of this type should prove espe- 

| cially valuable during the first 2 years of the rising 

part of a solar cycle, to help give an early indication 
of the level and timing of the coming maximum. 

It would of course be desirable to augment the 
sample size for estimating the regression equations of 
table 2 as more years of data become available, in 
particular the years 1961 and 1962, thus completing 
a 22-vear period. Also it should be noted that a 
similar thing could be done for any ionosphere 
sounding station having a long run of well-taken 
data. Thus delays in waiting for information from 
other parts of the world could be avoided. 

An indication of the degree of accuracy at certain 
other locations may be obtained from table 6. Esti- 
mated values of # for January 1961 from the regres- 
sion equations given in this table, are, for the respec- 

| tive stations, 70.6, 75.7, 81.4, and 84.5. 

It is suggested that adequate results might be ob- 
| tained from two locations, one in each hemisphere, 
| 
| 





at about the geomagnetic latitude of Washington, 
each station to cover its own winter months. For the 


standard errors for four stations, 1945 to 1960, Jan., July 





Regression equations » Standard errors 

Jan. July Jan. July 
=96.0+21.2( f-8.3) R=98.5+64.4( f-5.7) 13.9 27.7 
=96.0+20.3(/- R=98.5+60.7(f-6.7) 9.5 24.3 
=96.0+48.7( sR=102.4+2 22.2 2 18.9 
= 96.0+28.7( f-11.2) R=98.5+35. 15.8 17.3 





639 





equinoxes, possibly the station with the longer time 
series would be used. 


Finally, the presentation of two sets of equations | 


in table 2 (columns D and E) is intended to be illus- 
trative of two ways of obtaining the result. One 
should not necessarily switch from one type of equa- 
tion to another, because in small samples such as 
these the reversal in size of correlation coefficients or 
other criteria could happen frequently and not be 
significant. 
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Experimental measurements using scale model techniques have been carried out to 
determine the effectiveness of a ground system of long-wire radials to obtain low angles of 


departure of transmission. 


Since transmission was to be in one direction only, the ground 


wires were laid out to form a ground-plane sector approximately 18° wide centered in the 


direction of transmission. 


The antenna was a base-driven vertical monopole. 


Measure- 


ments were made of the relative response in decibels for the monopole used as a receiving 


antenna at a frequency of 400 megacycles per second. 
always located at a distance of 200 wavelengths. 


The target transmitter antenna was 
At this separation the ground plane 


sector was in the near field of the target transmitting antenna and appropriate corrections 


must be made. 


The received signal strength improvement due to the presence of the ground sector was 


approximately 14 decibels. 


The measured lobe positions of the first and second beam 


maximums and the first null are in good agreement with theory. 


1. Introduction 


In connection with a HF arctic propagation 
project using a base-driven vertical monopole at a 
frequency of 10 Me/s over a ground of lava rock 
having a poor conductivity, a requirement arose to 
evaluate the effectiveness of a planned long-wire 
ground radial system in improving low-angle 
response. 

A common practice has been to use a ground 
system of wire radials centered at the antenna base 
and extending outward for a distance of approxi- 
mately one-half wavelength [Laport, 1952]. How- 
ever, an improved ground system was desired to 
obtain low angles of departure of transmission. 

One method of obtaining a low angle of departure 
would be to elevate the base of the antenna above 
ground. However, this would result in multiple 
lobes in the vertical plane. It was then believed to 
be advantageous to improve on the effectiveness of 
the ground radial system and maintain the antenna 
base at ground level. 

Since the transmission path was in one direction 
only, a ground system was planned which would 
consist of wire radials 3,000 m long (approximately 
100 wavelengths) covering an approximate segment 
of ground 18° wide in the desired direction of 
transmission. 

Before the start of any actual construction, a 
study of this ground system of long-wire radials, 
using modeling techniques, was desired. Measure- 
ments were made to obtain the signal response over 
the center of the ground radial system and at other 
points over the ground-plane sector compared to the 
signal response over the ground path only. This 
would show the effect of the ground-plane sector 
width, the signal improvement due to the long wire 
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radials, and the improvement in response with 
different numbers of wire radials (density of wires 
in the ground-plane sector). 

This report presents the results of the measure- 
ments which were carried out at a frequency of 400 
Me/s using scaled model techniques. 


2. Method of Measurement 


For the purposes of measurement, the ground 
radial system to be measured was scaled for a 
frequency of 400 Me/s (a scaling factor of 40). 
The measurements were carried out at the NBS 
antenna range on Table Mesa near Boulder, Colo. 
The ground area over which the measurements 
were carried out was first cleared of weeds, brush, 
and large rocks in order to make the ground as 
smooth as possible. The ground conductivity is 
considered to be poor and no attempt was made to 
scale the ground conductivity. 

A quarter wavelength monopole, extended verti- 
cally from the center of a flat aluminum disk having 
a diameter of one wavelength. This disk was to 
represent a conventional ground system of a number 
of wires each one-half wavelength long. 

The long-wire ground system initially consisted of 
97 ground wires of No. 22 AWG Copperweld, each 
75 m long. These wires covered an approximate 
18° sector. The ends of the wires, separated by 
\/3(0.191° angular separation) were secured in 
position and bonded to 20 cm-long steel spikes which 
were driven into the ground until the tops of the 
spikes were flush with the ground surface. The 
other ends of the radials were bonded to the edge of 
the flat aluminum disk. 

The vertical monopole was used as a receiving 
antenna. The target transmitter antenna, a 5- 
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element Yagi oriented for vertical polarization, was 
moved in a 40° azimuthal are at a distance of 75 m 
in front of the wire radial ends. The whole ground 
plane sector was not in the far field of the trans- 
mitting Yagi antenna. The measurements were 
made by recording the relative received response in 
decibels by the monopole. Figure 1 is a plan view 
of the ground radial system. The response in decibels 
versus azimuth position of the transmitting antenna 
was Measured at a constant angle of arrival of 1.1°. 
The response in decibels versus angle of arrival in the 
range of 0° to 6.73° was measured at the sector 
center position. The angle of arrival is the included 
angle at the base of the monopole between the 
horizon and the target transmitter antenna. 

The measurements of response in decibels versus 
azimuth position of the transmitting antenna at a 
constant angle of arrival of 1.1° were also carried 
out for ground systems of 49 wires (0.382° angular 
separation) and 25 wires (0.764° angular separation) 
by removing every other wire from the original 97 
and then remaining 49 wires. 


3. Results 


Figure 2 presents the curves of recorded relative re- 
sponse in decibc's versus azimuth position of the trans- 
mitting antenna for the ground system consisting 

= or bd 9 = Oo 
of 97, 49, and 25 wires, measured over a 40° ground 
sector. 


POINTS DESIGNATE LOCI OF 
5-ELEMENT YAGI TARGET 
TRANSMITTER ANTENNA 





| 
| 
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when measured at azimuth positions over the ground- 


plane sector center as is to be expected. The 
received signal is also highest in the case when the 
ground-plane sector was of maximum density of wire 
radials (97 wires). The fall-off of the received 
signal when measured over a path of ground only 
(no long wire radials) is gradual, and at 20° each 
side of the ground-sector center the long wire radials 
do not appear to have any effect on the received 
signal strength. 

Fizure 3 is a plot of the relative response in decibels 
versus the angle of arrival in the forward direction 
with a ground-plane sector of 97 wires. Measure- 
ments at angies above 6.73° were not possible because 
of the maximum height of mounting for the target 
transmitter antenna and separation distance between 
the receiving monopole and transmitting antenna. 


4. Discussion of Results 


On the assumption that the wire radials form a 
well conducting ground plane, and from the geometry 
of figure 1, computations can be made to compare 
theoretically computed values [Wait and Conda, 
1958a] with the experimental results shown in figure 3. 

The computed position of the first beam maximum 
is 3.3°. The experimental value of this first maxi- 
mum is 2.4°. The position of the second maximum 


The received signal is generally higher | is computed to be 6.2° and the experimental value 


97 WIRES SPACED 4/3 AT 
ENDS COVERING SECTOR 
OF APPROXIMATELY 1|8.35° 


4/4 VERTICAL MONOPOLE 
eas ANTENNA 
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Figure 1. 


FLAT ALUMINUM DISK ONE 
WAVELENGTH DIAMETER 
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Plan view drawing of extended wire radial ground system modeled for measurement at 400 Mc/s. 
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FIGURE 2. 


Curves of relative response, decibels, versus azimuthal angle, degrees, for extended wire radial ground systems of 97, 49, 


and 26 wires. 


The angle of signal arrival is 1.10 in all cases. 
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Relative response of the monopole in the vertical plane in the forward direction over an cxtended wire radial system of 


97 wires. 


is 6.3°. The position of the experimental first null 
is at an angle of 4.3°. This compares to a computed 
first null position of 5.0°. 

With the exception of the first maximum, these 
experimental results agree quite well with a simple 
model of conducting half plane to represent the 
ground plane. The finite conductivity of the ground 
outside the ground plane will tend to lower the first 
lobe maximum [Wait and Conda, 1958a]. 

While the results of the model measurements are 
optimistic as regards low angle radiation, special 





considerations must be made when interpreting the 
results of the model measurements as related to the 
HF case. This is necessary since the angle of the 
lobes (for small departure angles) are approximately 
75 
1 150 where 
the long wire radials and 150 is the horizontal sepa- 
ration between transmitting antenna and the receiv- 
ing monopole (see fig. 1). If the separation between 
the transmitting antenna and receiving monopole 
increases from 150 m to infinity as in the field case, 


proportional to 75 is the length of 
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For this experiment 


the lot lobes are raised in elevation. 


1 ° 
yi- y < } Thus the experimentally measured lobe 


positions should be raised by more than 40 percent 
for optimum siting of a point-to-point HF link. 

Another consideration to be made when interpret- 
ing the scaled model measurements as regards the 
HF case is the influence of the earth curvature. At 
these low angles the pattern is cut back due to dif- 
fraction. For the 10 Me/s HF case this would amount 
to approximately 14 degree [Wait and Conda, 1958b]. 
Thus some error could be expected due to the earth 
curvature, particularly at low angles of the order of 
| degree. 

5. Conclusions 


The results of these measurements show that an 
increase in the density (number of wires in the sec- 
tor) improves the relative level of the received signal. 

Measurements of ground conductivity, effect of 
increased wire density, and effect of even longer wire 
radial ground systems were beyond the scope of 
these measurements. In order to completely in- 
vestigate the advantages of ground systems with 
extended wire radials, it would be desirable to make 
a further study with an additional number of wires 
and wire lengths, and to carry out measurements at 
higher angles of arrival. 


This work was carried owt in 1959 in support of 
the HF Arctic Propagation Project No. 85421 spon- 
sored by the U.S.N. BuShips. 

The measurements were made with the assistance 
of W. L. Martin, G. W. Angus, and R. W. Harsh, 
of the National Bureau of Standards, and W. E. 
McKenney and F. A. Foss of Scatter Communica- 
tion Corporation. 
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7. Appendix! 
7.1. Path Geometry 


A very important consideration for evaluating 
measurements as described in this paper is that the 
target transmitter is at a finite distance whereas 
the actual source of the field for a full size antenna 
is a downcoming plane wave (or at least something 
approaching it). Figure 4 shows the geometry of 
the path. 


1 This appendix is based on a private communication from J. R. Wait. 








~ 





FicurEe 4. Drawing of measurement path geometry. 


Transmitting antenna at point A; receiving monopole at point B. 


On the assumption that PB is a well conducting 
ground plane it can be expected that the antenna pat- 
tern first maximum in the vertical plane occurs when 


A=AP+PB—AB= 


of order unity. 


2A 2 
=m 4 when m? is some number 


By trigonometry the following will be obtained: 


ae l 
= 24. (qd—q)2]!/2+q— es 
[(d tan y)?+(d—a)’]"?+a cos Y 


If y is small this simplifies to 





ay’ 
. etc 
d 


A~- 


The angle y which satisfies A=m ; is 


Xr a 1/2 
yorm {x E -5} 


The value of m can be identified with the value of 
u where the function 


w=" {[5+eu 5+ |b. 


isa maximum. C and S are Fresnel integrals [Wait 
and Conda, 1958a]. Thus m—~1.2. 
For the measurements carried out at a frequency 


of 400 Me/s the first maximum is obtained. a@=75 m 
2a 
thus — tH14 and 
S X ( 
=| 
war (1 5)" radians 


~4m (1 -%)" degrees 


Yor2.8m~3.3° 


~2.2 which gives Y~6.2°. 
~5.0°. 


For the second maximum m 
For the first null m~1.8 which gives Y 





644 























7.2. Earth Curvature Influence 


Another consideration in evaluating the described 
measurements as pertained to an actual operational 
antenna is the influence of the earth’s curvature 
[Wait and Conda, 1958b]. 

At these low angles the pattern is ‘‘cut back’’ due 
to diffraction. The general restriction on neglecting 
this effect is then 








,9\13 


¥>>(z) radians 


where ka is the circumference of the earth in wave- 
lengths. For 10 Me/s this works out to be y>'s 
degree. Thus some error could be expected, par- 
ticularly at low angles of the order of 1 degree. 


(Paper 65D6-167) 
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Some boundary value problems involving plasma media, 
J. R. Wait, J. Research NBS 65B (Math. and Math. Phys.) 
No. 2, 1387-150 (April-June 1961). 

A plasma, consisting of a neutral mixture of electrons, ions, and 
molecules, in the presence of a constant magnetic field Hp, 
possesses a dielectric constant which isin the form of a tensor. 
Exact solutions of boundary value problems involving such 
media are obtained for two-dimensional configurations. 
Explicit results are given for the reflection coefficients of 
stratified plasma in planar and cylindrical geometry. 


A survey of spread-F, F. N. Glover, NBS TN 82 (PB 161583) 
(Nov. 1960) $1.78. 

Examples of spread-F forms occurring at different latitudes 
are presented, illustrating the classification of spread into 
range type and frequency type. The occurrence patterns of 
spread-F at different latitudes are correlated with other 
geophysical phenomena. Magnetic latitude and time withia 
the sunspot cycle appreciably affect the pattern of spread 
occurrence. Instrumental techniques and their advantages 
for spread studies are outlined. The principal theoretical 
explanations of spread-F are summarized. A single mech- 
anism need not be postulated as responsible for all types of 
spread occurrence or at all latitudes. 


A diffraction theory for LF sky-wave propagation, J. R. Wait, 
J. Geophys. Research 66, No. 6, 1713-1724 (June 1961). 

The concept that radio waves propagate from transmitter to 
receiver via discrete ionospheric reflections (i.e., hops) is put 
on a sound theoretical basis. The earth and the ionosphere 
are represented by two (smooth) concentric spherical surfaces. 
The (local) reflecting characteristics of these boundaries are 
assumed to be known. Itisshown that the mth hop sky wave 
may be generally represented by a complex integral. The 
geometrical-optical representations for the field are retrieved 
as a special case. For low radio frequencies, it is shown that 
these optical-type formulas, which are often used, have severe 
limitations. Generally, they are not valid near the caustic 
points. Techniques for evaluating the complex integral for 
the sky waves near and beyond the caustics are then discussed 
in some detail. Where possible, a physical interpretation of 
the results is given. 


A diffraction theory for LF sky-wave propagation—an addi- 
tional note, J. R. Wait and A. M. Conda, J. Geophys. Research 
66, No. 6, 1725-1729 (June 1961). 

The behavior of the first-hop sky wave near its geometrical 
horizon is discussed. Using formulas already developed, 
numerical results in this difficult caustic region are presented 
in graphicalform. Itis shown that, sufficiently far away from 
the caustic, the field approaches its geometrical-optical value 
The results should be useful for sky-wave field computation 
at low frequencies and in the interpretation of measured 
field-strength data. 


Photoionization heating in the F region of the atmosphere, 
D.C. Hunt and T. E. Van Zandt, J. Geophys. Research 66, 
No. 6, 1673-1682 (June 1961). 

The thermal economy in the F region of the atmosphere is 
investigated using the barometric equation and the continuity 
equation for thermal energy. It is concluded that photoioni- 
zation is certainly a major source, and possibly the predomi- 
nant source, of F-region heating. The assumptions and ap- 
proximations are discussed and are shown not to invalidate 
the conclusion. The heat flux due to photoionization in the 
F region at low and middle latitudes at sunspot maximum is 
estimated to lie between 0.5 and 2.0 ergs/em? sec, with a most 
probable value of about 1.0 erg/cm? sec. 





The use of interferometer observations of satellites for 
measurement of irregular ionospheric refraction, R. S. Law- 
rence and J. W. Warwick, Annals of the Intern. Geophys. 
— Pt. 2, pp. 566-569 (Pergamon Press, New York, N.Y., 
1961 

We outline a digital method for determination of satellite 
orbital parameters from a single passage through the lobes of a 
simple interferometer. Differences between the observed and 
computed angle-of-arrival indicate that, at least for low- 
altitude passages, the ionospheric irregularities probably occur 
in mas layer and are not strongly aligned along the magnetic 
field. 


Long-distance one-hop F, propagation through the auroral 
zone, L. H. Tveten, J. Geophys. Research 66, No. 6, 1683- 
1684 (June 1961). 

The existence of one-hop F; propagation through the auroral 
zone over a 4,495-km path is described, and an example of a 
sweep-frequency record illustrating this F; mode is shown. 


Incoherent scattering by free electrons as a technique for 
studying the ionosphere and exosphere: Some observations 
and theoretical considerations, K. L. Bowles, pp. 223-242 of 
AGARDograph 42, The upper atmosphere above F2-maxi- 
mum, North Atlantic Treaty Organization, Advisory Group for 
Aeronautical Research and Development, Avionics Panel. 
Papers presented at the Symp. Ionospheric Research Com- 
mittee, Paris, France, May 1959. (Technical Editing and 
Reproduction Ltd., 95 Great Portland Street, London, W. 1, 
England.) 

Incoherent scattering by the free electrons of the ionosphere 
has been suggested as a technique for measuring the electron 
density profile both below and above the F-region maximum. 
This paper reports observations which confirm the existence 
of the incoherent scatter and show that its intensity is es- 
sentially the predicted value. The observed Doppler broad- 
ening is considerably smaller than predicted. In the second 
part of the paper, an explanation for the reduced Doppler 
broadening is offered. The scatter is explained as arising 
from statistical fluctuations of electron density, the distribu- 
tion of which is controlled by the positive ions. 


A model of the F-region above h,,.:/2, J. W. Wright, 
211-222 of AGAR Dograph 42, The upper atmosphere ie 
F2-maximum, North Atlantic Treaty Organization, Advisory 
Group for Aeronautical Research and Development, Avionics 
Panel. Papers presented at the Symp. Ionospheric Research 
Committee, Paris, France, May 1959. 

A simple Chapman model (seale height 100 km) of the 
ionosphere above the peak of the F-region is found to be in 
good agreement with the few observed profiles of this region. 
The ratio of the above-peak electron content implied by the 
model, to below-peak electron content, is found to vary from 
about 2.8 to 4.0, also in good agreement with observations. 
The model is used in conjunction with electron density pro- 
files observed from the ground to construct meridional cross- 


ro 


sections along the 75° W geographic meridian. 


A comparison between theoretical and experimental data on 
phase velocity of VLF radio waves, J. R. Wait, Proc. IRE 
49, No. 6, 1089-1090 (June 1961). 

It can be concluded that there is fairly good agreement 
between measured phase velocity and the mode theory as 
developed for a curved earth-ionosphere waveguide. Fur- 
thermore, the heights of the equivalent reflecting layer are 
fairly consistent with other results. 
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Journal of Research NBS 65A (Phys. and Chem.) No. 5 
(September—October 1961) 75 cents. 


Calibration of a monitor for use in bremsstrahlung beams. 
G. Fuller and E. Hayward. 

Mass spectrometric study of NF2, NF3, N2F2, and NoF,. 
J. T. Herron and V. H. Dibeler. 

Rate of the reaction NO+N, and some heterogeneous re- 
actions observed in the ion source of a mass spectrometer. 
J. T. Herron. 

Synthesis of the A. Van 
Valkenburg. 

Phase equilibria in systems involving the rare earth oxides. 
Part III. The Eu.0;-In,O3 system. 8S. J. Schneider. 

Heats of hydrolysis and formation of dimethoxychloroborane. 
M. V. Kilday, W. H. Johnson, and E. J. Prosen. 

Tritium-labeled compounds VII. Isotope effects in the oxi- 
dation of p-mannitols-C'* and p-mannitols-t to p-fructoses. 
L. T. Sniegoski, H. L. Frush, and H. S. Isbell. 

Franck- oe factors to high vibrational quantum num- 
bers I: N. and N}. R. W. Nicholls. 





humites nMg.SiO, Mg(F,OH) >. 


Journal of Research NBS 65B (Math. and Math. Phys.) No. 

3 (July-September 1961) 75 cents. 

Theory of an accurate intermediary orbit for satellite as- 
tronomy. J. P. Vinti. 

Note on the “baffled piston”’ 

Some results on non-negative 
Mine, and B. N. Moyls. 

Probability inequalities of the 
Savage. 


problem. F. Oberhettinger. 
matrices. M. Marcus, H. 


Tchebycheff type. I. R. 


Journal of Research, NBS 65C (Eng. and Instr.) No. 3 
(July-September 1961) 75 cents. 


Prediction of symptoms of cavitation. R. B. Jacobs. 

Heating and cooling of air flowing through an underground 
tunnel. B. A. Peavy. 

Stress-corrosion cracking of the 

Logan. 

Coatings formed on steel by cathodic protection and their 
evaluation by polarization measurements. W. J. Schwerdt- 
feger and R. J. Manuele. 

Calibration of inductance standards in the Maxwell-Wien 
bridge circuit. T. L. Zapf. 
Calibration of loop — at VLF. A. G. Jean, H. E. 
Taggart, and J. R. Wait. 

Location of the plane be best average definition with low 
contrast resolution patterns. F. E. Washer and W. P. 
Tayman. 

Influence of temperature and relative humidity on the photo- 
graphic response to Co gamma radiation. M. Ehrlich. 


AZ31P magnesium alloy. 


Ideal gas thermodynamic functions and isotope exchange 
functions for the diatomic hydrides, deuterides, and tri- 
trides, L. Haar, A. 8. Friedman, and C. W. Beckett, NBS 
Mono. 20 (1961) $2.75. 

Bibliography of temperature measurement— 
June 1960, C. Halpern and R. J. Moffat, 
(1961) 15 cents. 

Causes of variation in chemical analyses and physical tests 
of portland cement, B. L. Bean and J. R. Dise, NBS Mono. 
28 (1961) 25 cents. 

Thermal expansion of technical solids at low temperatures. 
A compilation from the literature, R. J. Corruccini and 
J. J. Gniewek, NBS Mono. 29 (1961) 20 cents. 

Corrected optical pyrometer readings, D. E. Poland, J. W. 
Green, and J. L. Margrave, NBS Mono. 30 (1961) 55 cents. 

Publications of the National Bureau of Standards July 1, 
1957, to June 30, 1960 (Includes Titles of Papers Published 
in Outside Journals 1950 to 1959), B. L. Arnold, NBS 
Mise. Publ. 240 (1961) $2.25. 

Mean electron density variations of the quiet ionosphere 
May 3, 1959, J. W. Wright, L. R. Wescott and D. J. 
Brown, NBS TN40-3 (PB151399-3) (1960) $1.50. 
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A survey of computer programs for chemical information 
searching, E. C. Marden and H. R. Koller, NBS TN85 
(PB161586) (1961) $2.25. 

Integrated starlight over the sky F. E. Roach, and L. R. 
Megill. Astrophys. J. 133, No. 1, 228-242 (Jan. 1961). 

The electromagnetic fields of a dipole in the presence of a 
thin plasma sheet, J. R. Wait, Appl. Sci. Research 8 
sec. B, 397-417 (1960). 

The nz ture, cause and effect of the porosity in electrodeposits, 
VII. A microscopic examination of nickel-chromium coat- 
ings after atmospheric corrosion, D. W. Ernst and F. 
Ogburn, Plating 48, No. 4, 491-497 (May 1961). 

Physical quantities proposed for radiation measurements, L. 
S. Taylor, Proce. IXth Intern. Congress of Radioiogy, July 
23-30, 1959 (Munchen, Germany), 1301-1302 (1960). 

Improving rectifier circuits, G. F. Montgomery, Electronics 
p. 86-87 (April 7, 1961). 

Fracture characteristics of notched tensile specimens of 
titanium and a titanium alloy, G. W. Geil and N. L. Car- 
wile, Materials Research and Standards (Am. Soc. Testing 
Materials, Philadelphia, Pa.), 1, No. 1, 16-21 (Jan. 1961). 

The robustness of life testing procedures derived from the 
exponential distribution, M. Zelen and M. C. Dannemiller, 
Technometrics 3, No. 1, 29-49 (Feb. 1961). 

Absence of an isotope effect in the fractional recrystallization 
of alpha-D-glucose-1-t, H. 8. Isbell, H. L. Frush, and N. B. 
Holt, Anal, Chem. 33, 225-226 (Feb. 1961). 

Applications of graphs and Boolean matrices to computer 
programming, R. B. Marimont, SIAM Rev. 2, No. 4, 
259-268 (Oct. 1960). 

Comparison of national standards for roentgen measurement, 
H Wyckoff, IXth Intern. Congress of Radiology, July 
23-30, 1959 (Munchen, Germany), 1315-1318 (1960). 

Environmental factors in a family-size underground shelter, 
R. P. Achenbach, Proc. Meeting on Environmental Engi- 
neering in Protective Shelters, Feb. 8-10, 1960, Natl. Acad. 
Sci.-Natl. Research Council, Wash. D.C., 69-118 (1961). 

How to evaluate accuracy, W. J. Youden, Materials Re- 
search and Standards (Am. Soc. Testing Materials, Phila- 
delphia, Pa.), 1, No. 4, 268-271 (April 1961). 

Improved standard for the calibration of vibration pickups, 
R. R. Bouche, Exp. Mech. 1-6 (April 1961). 

a and seatte ring of photons by holmium and erbium, 
E. Fuller and E. Hayward, Proc. Intern. Conf. Nuclear 
ed coy 763-766 (Kingston, Canada, 1960). 

Introductory remarks, U. Fano, Conf. Coherence Properties 
of Electromagnetic Radiation, AFOSR-583, Report on 
Conf., 31-38 (University of Rochester, The Institute of 
Opties, Rochester, N.Y., April 1961). 

Available heat sinks for protected underground installations, 
B. A. Peavy, Proc. Meeting on Environmental Engineering 
in Protective Shelters, Feb. 8-10, 1960, Natl. Acad. Sci.- 
Natl. Research Council, Wash. D.C., 69-118 (1961). 

Interpretation of the appearance potentials of secondary ions, 
M. B. Wallenstein and M. Krauss, J. Chem. Phys. 34, No. 
3, 929-936 (Mar. 1961). 

Observations on the chemiluminescent reaction of nitric 
oxide with atomic oxygen, H. P. Broida, H. I. Schiff, and 
T. M. Sugden, Trans. Faraday Soc. 57, No. 458, 259-265 
(Feb. 1961). 

Studies of the low-temperature distillation of hydrogen 
isotopes, T. M. Flynn, Cryogenics 1, No. 2, 1-5 (Dee. 
1960); Advances Cryogenic Engr. 6, 236-244 (1961). 

Inert enclosed pump for shaped flow of ultraclean solutions, 

. W. Barnes and R. H. Noyce, Rev. Sci. Instr. 32, No. 
3, 353 (Mar. 1961). 

On some partial differential equations of Brownian motion of 
a free particle, A. Ghaffari, (Abstract) Proce. Intern. Conf. on 
Partial Differential Equations and Continuum Mechanics 
(Univ. of Wisconsin, Madison, Wise., June 7-15, 1960), 
348-350 (1961). 

Radiation protection standards, L. 8. Taylor, Radiology 74, 
824-831 (May 1960). 


Atomic beam frequency standards, R. C. Mockler, R. E. 
Beehler, and C. 8. Snider, IRE Trans. Instrumentation 
I-9, No. 2, 120-132 (Sept. 1960). 


The development of more stable ~ blocks, M. R. Meyerson, 
T. R. Young, and W. R. Ney, ASTM Bull. 1, No. 5, 
368-374 (May 1961) 


648 








The power spectrum and its importance in precise frequency 
measurements, J. A. Barnes and R. C. Mockler, IRE Trans. 
Instrumentation I-9, No. 2, 149-155 (Sept. 1960). 

VLF phase perturbation associated with meteor shower 
ionization, C. J. Chilton, J. Geophys. Research 66, No. 2, 
379-383 (Feb. 1961). 

Departures from the Saha equation for ionized helium. 

Condition of detailed balance in the resonance lines, 
R. N. Thomas and J. B. Zirker, Astrophys. J. 183, No. 2, 
588-595 (Mar. 1961). 

Tests for regression coefficients when errors are correlated, 
M. M. Siddiqui, Annals Math. Stat. 31, No. 4, 929-938 
(Dec. 1960). 

Electron diffraction studies on solid a-nitrogen, E. M. Horl 
and L. Marton, Acta Cryst. 14, Pt. I, 11-19 (Jan. 10, 
1961). 

A comment on the NRL solar Lyman-Alpha results, J. T. 
Jeffries and R. N. Thomas, Astrophys. J. 1338, No. 2, 
606-607 (Mar. 1961). 

The nation’s electronic standards program: Where do we 
now stand? H. W. Lance, IRE Trans. Instrumentation 
I-9, No. 2, 94-100 (Sept. 1960). 

Hard gallium alloys for use as low contact resistance electrodes 
and for bonding thermocouples into samples, G. C. Harman, 
Rev. Sci. Instr. 31, No. 7, 717-720 (July 1960). 

Seattering of photons by deformed nuclei, E. G. Fuller and 
E. Hayward, Proc. Intern, Conf. Nuclear Structure, 763-766 
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Solar flare effects in the F region of the ionosphere (letter 
to the editor) R. W. Knecht and K. Davies, Nature 190, 
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Summary of outstanding solar events in July 1959, A. H. 
Shapley and Dorothy Trotter, Union Géodésique et 
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1960 (L’Institut Geographique National, Paris, 1961), 
pp. 72-76. 

Solar disturbances and radio communication forecasts 
James F. Brockman, Sky and Telescope XXI, No. 6, 
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